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We present a unifying framework to study physical systems wfiicli exfiibit topological quantum 
order (TQO). The major guiding principle behind our approach is that of symmetries and entan- 
glement. To this end, we introduce the concept of low-dimensional Gauge-Like Symmetries (GLSs), 
and the physical conservation laws (including topological terms, fractionalization, and the absence 
of quasi-particle excitations) which emerge from them. We prove then sufficient conditions for TQO 
at both zero and finite temperatures. The physical engine for TQO are topological defects associated 
with the restoration of GLSs. These defects propagate freely through the system and enforce TQO. 
Our results are strongest for gapped systems with continuous GLSs. At zero temperature, selection 
rules associated with the GLSs enable us to systematically construct general states with TQO; these 
selection rules do not rely on the existence of a finite gap between the ground states to all other 
excited states. Indices associated with these symmetries correspond to different topological sectors. 
All currently known examples of TQO display GLSs. Other systems exhibiting such symmetries 
include Hamiltonians depicting orbital-dependent spin-exchange and Jahn- Teller effects in transi- 
tion metal orbital compounds, short-range frustrated Klein spin models, and p-|-ip superconducting 
arrays. The symmetry based framework discussed herein allows us to go beyond standard topo- 
logical field theories and systematically engineer new physical models with finite temperature TQO 
(both Abelian and non-Abelian). Furthermore, we analyze the insufficiency of entanglement entropy 
(we introduce SU{N) Klein models on small world networks to make the argument even sharper), 
spectral structures, maximal string correlators, and fractionalization in establishing TQO. We show 
that Kitaev's Toric code model and Wen's plaquette model are equivalent and reduce, by a duality 
mapping, to an Ising chain, demonstrating that despite the spectral gap in these systems the toric 
operator expectation values may vanish once thermal fluctuations are present. This illustrates the 
fact that the quantum states themselves in a particular (operator language) representation encode 
TQO and that the duality mappings, being non-local in the original representation, disentangle the 
order. We present a general algorithm for the construction of long-range string and brane orders in 
general systems with entangled ground states; this algorithm relies on general ground states selec- 
tion rules and becomes of the broadest applicability in gapped systems in arbitrary dimensions. We 
discuss relations to problems in graph theory. 

PACS numbers: 05.30.-d, 11.15.-q, 71.10.-w, 71.10.Pm 



I. INTRODUCTION 

The emergence of complex phenomena in strongly cou- 
pled quantum matter may be the result of quite simple 
guiding principles yet to be discovered. Nevertheless, we 
expect that the key players of Symmetry and Topology 
are rooted in those principles, and their cross fertiliza- 
tion may help to unveil those guiding principles or even 
lead to new phenomena. That symmetry is a guiding 
principle is demonstrated by the classical Landau the- 
ory of phase transitions [l| . Two key contributions form 
the crux of this theory: one is the association of a local 
order parameter 0(r) to a broken (ordered) symmetry 
state, thus to (the lack of) symmetry, and the other is a 
phenomenological description of the transition in terms 
of an analytic functional of 0{r) (Landau free energy). 
Despite its simplicity and success, can all possible phases 
of matter and their transitions (even approximately) be 
described within Landau's framework? 

A new paradigm, Topological Quantum Order (TQO), 
extends the Landau symmetry-breaking framework. At 
its core, TQO is intuitively associated with insensitiv- 
ity to local perturbations: the order is topological. As 



such, TQO cannot be described by local order parame- 
ters. Some examples of systems displaying TQO also dis- 
play ground state (GS) subspace degeneracy highly sensi- 
tivity to boundary conditions or surface topology. What 
are the physical principles and consequences behind such 
an order? How do we mathematically characterize it? 
Although much has been learned about this topic, much 
still remains shrouded in mystery, in part, because there 
is no universally accepted definition of TQO from which 
many results may then be shown to follow. Further inter- 
est in TQO is catalyzed by the prospect of fault-tolerant 
topological quantum computation 

Notwithstanding the lack of an unambiguous accepted 
mathematical definition, several inter-related concepts 
are typically invoked in connection to TQO: symmetry, 
degeneracy, fractionalization of quantum numbers, max- 
imal string correlations {non-local order), among others. 
A fundamental question is what is needed to have a sys- 
tem exhibiting TQO. The current article aims to show 
relations between these different concepts, by rigorously 
defining and establishing the equivalence between some 
of them and more lax relations amongst others. Our main 
contribution is to provide a unifying symmetry principle 
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to characterize existent physical models and engineer new 
systems displaying TQO. Clearly, this principle leads to 
unique physical consequences that we will expand on be- 
low. Most importantly, we (i) prove that systems harbor- 
ing generalized d-dimensional (with d = 0,1,2) Gauge- 
Like Symmetries (d-GLSs) exhibit TQO; (ii) analyze the 
resulting conservation laws and the emergence of topo- 
logical terms in the action of theories in high space di- 
mensions (dimensional reduction); (iii) affirm that the 
structure of the energy spectrum is irrelevant for the ex- 
istence of TQO (the devil is in the states); (iv) establish 
that, fractionalization, string correlators, and entangle- 
ment entropy are insufficient criteria for TQO; (v) report 
on a general algorithm for the construction of string and 
brane correlators; (vi) suggest links between TOO and 
problems in graph theory. The current work [3| is an 
elaborate extension of Rcf. Q in which our central re- 
sults were stated yet due to the lack of space were not 
proven or explained in detail. 

The present work aims to investigate many recurring 
themes in the study of TQO. We start in Section HIl bv 
classifying the degrees of locality/non-locality of general 
operators acting on a physical Hilbert space, and the pos- 
sible spectral structures. 

In Section IIIIl we define both zero and finite temper- 
ature TQO. Thus far, most work on TQO focused on 
the GS properties of (gapped) quantum systems. In this 
work, we formalize a notion of finite temperature TQO. 
This notion is based on the influence of boundary condi- 
tions on both local and topological observables. 

Then, in Section IIVI we review the existence of an 
SU{Ng) type symmetry (or of subgroups thereof) in a 
general system with an iVg-fold degenerate GS. As well 
known, SU (N) symmetry can lead to a fractionalization 
of the basic charge (e.g., the quark charges which are 
multiple of e/3 in SU{3) or N-ality in the general case). 
Thus degeneracy can lead to fractionalization in much 
the same way as the same groups lead to fractionalization 
elsewhere. A similar link was also looked at anew lately 
by Oshikawa and coworkers 0. 

Next, in Section |Vl we introduce one of our key re- 
sults. First, we introduce the concept of GLSs - symme- 
tries that generally lie midway between local symmetries 
(standard gauge symmetries) and global symmetries. In 
general, these symmetry operators act on a d-dimensional 
region. Here, d = correspond to local gauge symmetries 
while d = D operations with D the dimensionality of the 
system correspond to global symmetry operations. Inter- 
mediate GLSs allow for < d < D. In Section IVIl we 
review known examples of TQO and analyze the GLSs 
which are present in each of these systems. In Section 
IVIIl wc introduce several other systems which exhibit 
GLSs. In Section IVIIIl we make links between these 
symmetries and Aharonov-Bohm-type unitary transfor- 
mations. 

In Section HXl wc study the physical consequences of 
having a system endowed with d-dimensional GLSs. In 
particular, we discuss the appearance of new conserved 



Noether and topological charges in high-dimensional sys- 
tems as a result of dimensional reduction. We also ex- 
plain the conditions to break those symmetries, and show 
that well-defined quasi-particle poles are precluded in 
systems that display low-dimensional GLSs. Finally, we 
comment on the relation between charge fractionaliza- 
tion and GLSs, and the appearance of low-dimensional 
topological terms in higher-dimensional theories. 

In Section |Xl we show that there exists an intimate 
relation between GLSs and TQO: GLSs of sufficiently 
low-dimensionality d mandate the appearance of TQO. 
At finite temperature, the proof relies on a recent exten- 
sion of Elitzur's theorem. The results are strongest for 
gapped systems which harbor continuous d < 2 GLSs 
(Section |XB|; here, the requisite T = TQO can be 
proven very generally. 

As we show in Section IXII in other (also potentially 
non-gapped) systems with either continuous and discrete 
GLSs, at zero temperature, the results are also a conse- 
quence of selection rules. These rely on selection rules 
(such as the Wigner-Eckart theorem) as applied to the 
GLS generators. As a particular corollary, these GLSs 
selection rules imply (for Hamiltonians which contain 
finite-range interactions) a degeneracy of energy levels 
in the thermodynamic limit. These selection rules fur- 
ther enable us to construct special GLSs eigenstates in 
which T — TQO appears. To date, the prominent ex- 
amples of zero temperature TQO have GSs which are of 
that special form. In Section IXI II wc illustrate how these 
special states indeed appear as the GSs of the Kitaev and 
related models. Moreover, we show how to systematically 
engineer models that have T = TQO. 

In Section [XIIIl we expand on the general equivalence 
in the spectral structure between TQO systems on the 
one hand and systems with standard local orders on the 
other [J . This affirms that there is no spectral signature 
of TQO. In this Section, we prove that although Kitaev's 
model exhibits TQO, at any finite temperature the Toric 
code operators have vanishing expectation values. This 
occurs notwithstanding the fact that the Kitaev model 
has a finite gap between the ground and all other ex- 
cited states. A spectral gap is not sufficient to ensure 
that topological quantities do not have a vanishing ex- 
pectation value at finite temperature. The crux of this 
demonstration is that the Kitaev model (being a sum 
of commuting operators) can be mapped onto the one- 
dimensional {D ~ 1) Ising model. Similar considerations 
apply to Wen's plaquette model (which as we show is 
none other than a rotated version of Kitaev's model). 
These results may have important implications for Topo- 
logical Quantum Computation. 

In Section lXIVl we examine the entanglement entropy. 
We will show that the quantitative identification of TQO 
as a deviation from an area law form recently suggested 
by several authors differs from the original definition of 
TQO as a robust topological order. At best, all states in 
which TQO appears in the entanglement entropy should 
be a subset of all robust topological orders which are in- 
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sensitive to all quasi-local perturbations. To make our 
arguments sharper we introduce a new class of SU (N) 
Klein models on a small worlds network. 

In Section IXVl we make analogies between the char- 
acterization of TQO orders and that of identifying the 
topology of a graph. In the Appendix, we carry this 
analogy further and construct wavefunctions on graphs 
for which TQO appears that we dubbed Gauge-Graph 
Wavefunctions. 

In Section [XVIl we discuss one-dimensional string and 
higher-dimensional brane correlations. Such orders ap- 
pear in the 5=1 AKLT chains as well as in other elec- 
tronic systems (e.g., doped Hubbard chains). We show 
a general algorithm for the construction of such correla- 
tions. Our approach is detailed for general gapped sys- 
tems (in any dimension) as well as in general entangled 
systems with GS selection rules. In many cases, this con- 
struct leads to the identification of a gauge- like structure. 
This structure should not be confused with the GLSs. It 
is not a symmetry of the system - rather it corresponds 
to a transformation between the GSs of the systems to 
a uniform factorizable state. We show that systems (e.g. 
the AKLT spin chain) can display string orders yet not 
exhibit TQO. Moreover, we show that the hidden non- 
local string order in the AKLT problem is equal to the 
expectation value of a local nematic-type order. 

In Section IX VI I [ we illustrate that systems can have a 
degeneracy which depends on the topology of the mani- 
fold on which they are embedded yet not be topologically 
ordered in the sense of robustness to local perturbations. 
Historically, such a dependence of the GS degeneracy on 
the topology of the manifold on which the system is em- 
bedded sparked much of the study of TQO. 

We conclude in Section [XVIIIl with a summary of our 
results. Technical but quite important details (some of 
them novel) have been relegated to the various Appen- 
dices. These include Sections on degenerate perturbation 
theory, generalized Elitzur's theorem, the Peierls' prob- 
lem and derivation of connecting symmetry operators, re- 
lation between the quantum dimer and Kitaev's models, 
spin S — I chains and the derivation of string operators 
as partial polarizers. 



II. CLASSIFICATION OF LOCAL AND 
NON-LOCAL OPERATORS 

In the bulk of this work, we will focus on quantum 
lattice systems (and their continuum extension) which 
have Ns = Y[^=i sites, with the number of sites 
along each space direction and D the dimensionality 
of the lattice A. Unless stated otherwise, we will assume 
that = L for all directions /i. Associated to each 
lattice site (or mode, or bond, etc.) i € there is 
a Hilbert space Tii of finite dimension T). The Hilbert 
space is the tensor product of the local state spaces, Ji = 
{^;7ii, in the case of distinguishable subsystems, or a 
proper subspace in the case of indistinguishable ones. On 



a few occasions, we will briefly introduce extensions of 
our analysis to several network and graph systems. 

Statements about local order, TQO, fractionalization, 
entanglement, etc., are made relative to the particular de- 
composition used to describe the physical system. Typ- 
ically, the most natural local language Q is physically 
motivated: This is the essence of local Quantum Field 
Theory. For the present purposes we will classify opera- 
tors according to the following: 

(a) Local: Oi (or finite linear combinations). 

(b) Quasi-local: Hiez^J^s i'^^ finite linear combi- 
nations), where Ms represents a finite (non-extensive) 
integer number (even in the thermodynamic limit). 

(c) Non-local: HigzAfs (•-"" finite linear combina- 
tions), where an extensive number Ms < Ns of lattice 
sites are involved (e.g., Ms = Li in a _D-dimcnsional lat- 
tice). 

(d) Global; '^'^ extensive sum of quasi- 
locals or non-locals. 

(e) Quasi-global: X^iez^^^s '-^i' f-^*^ extensive sum 
of quasi-locals or non-locals. Ms < Ns is an extensive 
integer number. 

This classification characterizes either the product 
(cases (a)-(c)) or the sum of local operators (cases (d)- 
(c)). For physical systems not represented by an Ti, with 
a tensor product structure, there is no natural subsys- 
tem decomposition. In general one needs the concept 
of generalized entanglement Q to understand the nature 
of the correlations of their quantum states. In particu- 
lar, as we will see below, the nature of entanglement in 
systems which harbor TQO is of a highly non-local char- 
acter (with respect to the local language) ^ as opposed to 
systems with local order parameters. 

The dynamics of the physical system will be governed 
by a generic Hamiltonian H whose form is constrained, 
in terms of the local language Q , to linear combinations 
of (polynomial in Ns) quasi-local operators. 

To put the subject of TQO in perspective let us start 
by tabulating all of the possible low-energy spectra and 
discuss their relationship to different states of matter. 
Figure [T] shows a schematics of the possible low-energy 
spectra which is self-explanatory. Basically, the spectrum 
may be gapless or gapped and the GS Hilbert subspace 
may be degenerate or not. The main case of interest for 
the present manuscript is case (6). As we will see, how- 
ever, in order to have a TQO state it is not at all manda- 
tory to have an energy spectrum like (6) - symmetry alone 
can potentially mandate that gapless systems also exhibit 
TQO. The Lieb-Schulz-Mattis theorem pertaining to one 
dimension Q and a recent extension by Hastings Q in 
higher dimensions preclude, in the thermodynamic limit, 
rather general (locally interacting) 5=1/2 systems from 
having energy spectra of type (c) shown in Fig. [TJ 
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FIG. 1: Schematic representation of potential energy spectra. 
Case (a) may represent a Fermi liquid, case (fe) a Fractional 
Quantum Hall liquid, case (c) a band insulator, and case (d). 



III. WHAT IS TOPOLOGICAL QUANTUM 
ORDER (TQO)? 

To determine what is needed for TQO, we start by 
defining it. Various discussions in the Uterature empha- 
size the characteristics between only a pair of states while 
others aim to link TQO to the general form of the GS 
sector. To date, all of the treatments seen in the litera- 
ture focused on zero temperature (T = 0) TQO. We will 
discuss r = TQO below and then define its finitc-T ex- 
tension. In later Sections, we provide results on finite-T 
TQO. 



A. Zero temperature TQO 

To allow for greater flexibility and precision, we now 
define rank-n TQO. Given a set of n orthonormal GSs 
{\ga)}a=i,--- ,n, with 1 < n < A'g where Ng is the total 
number of GSs of a given Hamiltonian, rank-n TQO ex- 
ists iff for any bounded operator V with compact support 
(i.e. any quasi-local operator V), 

{9a\V\g,3) = V Sap + Caf3, (1) 

where w is a constant and Caf3 is a correction that vanishes 
in the thermodynamic limit. The physical motivation for 
this definition lies in the irrelevance of finite order terms 
in degenerate perturbation theory for any quasi-local op- 
erator V. (In Appendix \^ we review the rudiments of 
degenerate perturbation theory.) We will now demon- 
strate a simple result concerning the T = definition of 
Eq. ID) and then turn to the finite-T extension. 



Lemma. Given a set of orthonormal GSs {|5q)} which 
span an n-dimensional space, condition ([T]) holds iff for 
any GS in the same space, \g) = J^a «a|ffa) {J2a = 
1), we have that 

{g\V\g}^v + Cg, (2) 

with Cg a correction that vanishes in the thermodynamic 
limit. 

Proof: Let us define the operator Vq ~ PqVPq where 
Pq is the projection operator onto the GS basis. Condi- 
tion ^ then implies that (fflVblg) = v + Cg. If this holds 
true for any \g), it is also true for an orthonormal set 
{\ga)} that diagonalizcs Vq, i.e. Volga) = {v + Ca)\ga)- In 
this basis Vq = X]a(^' + Ca)\ga){ga\- Then, given another 
set of orthonormal GSs [generally, linked by a unitary 
transformation to the eigenstates of Vq] \ga) = X]/3 '^/sls/s) 

{9a\V\gp) = {galVolgp) = v Sap + c^p, (3) 

with vanishing Cap — J^j^-yi'^")*'^^- Similarly, ^ triv- 
ially implies {g\V\g) J2a,pM* ap{ga\V\gp) = 
V + Cg, with Cg ^J2a.pi"'a)*apCap- ■ 

This simple Lemma naturally lends itself to a finite-T 
generalization which we turn to next. 

In passing, we note that a similar but not as easily 
generalizable relation to finite-T follows from the off- 
diagonal portion of Eq. ([T]): Equation ([1]) is satisfied if 
and only if {ga\V\gp) ~ Cap for any two orthogonal states 
\ga) and \gp) which lie in the GS manifold. 

B. Finite temperature TQO 

Although the definition of TQO is cleanest for the 
T = problem, in all physical problems of relevance, 
we must consider the effect of finite temperatures. In 
we introduced for the first time, a finite-T (T > 0) gener- 
alization of TQO. Extending the general T = condition 
for TQO of Eq. we have that 

{V)a^tT{paV) = V + CaaiL), (4) 

(independent of a) with Caa a correction that tends 
to zero in the thermodynamic limit and pa = 
expl—Ha/ikBT)] a density matrix corresponding to the 
Hamiltonian H endowed with boundary terms which fa- 
vor order in the state \ga) [I^l- These boundary terms 
are exactly of the same form of perturbations that allow 
us to manipulate the system initially and set to be in one 
or the other topological GSs at T = for quantum com- 
puting purposes. In later Sections, we will show that the 
excitation spectrum precisely associated with non-local 
d = 1 discrete symmetry operators or d ~ 1,2 continuous 
symmetries (which is dominated by topological defects) 
enables finite-T TQO to appear. We define finite-T TQO 
to hold iff Eq. Q holds at finite temperatures on those 
systems which also exhibit the T = TQO defined by 
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Eq. ID). As will become apparent in later Sections (see, 
e.g., Section IXIII Ap . the general finite-T robustness con- 
dition of Eq. (j4|) by which we define TQO is different 
from finite expectation values of non-local topological op- 
erators; the latter often capture T = TQO in many 
known instances. 

What is the relation between the definition of TQO 
above and its sensitivity to boundary conditions? To 
avoid the use of spurious perturbations such as singular 
projection operators onto a particular GS, we define the 
perturbations to correspond to boundary conditions as- 
sociated with a particular GS. A system exhibits finite-T 
TQO if it (i) obeys the T = TQO conditions of Eq. 
(H]) and (ii) satisfies Eq. (gl) for T > with the expec- 
tation value subscript a corresponding to the constraint 
of boundary conditions which favor the state a. If the 
expectation value of any quasi-local operator V is inde- 
pendent of the boundary conditions a then the system 
exhibits TQO: the topological information on a encoded 
in the boundary conditions is not accessible (and pro- 
tected from) any local measurement. In that sense, TQO 
systems are similar to disordered systems in which no lo- 
cal order is found due to a condensation of topological 
defects. 

To make the definition of finite-T TQO clear and intu- 
itive, let us consider two well-known systems - the U{1) 
gauge theory which exhibits TQO and the Ising ferro- 
magnet in dimensions D > 2 which does not exhibit 
finite-T TQO. Within the U{1) gauge theory, if we fix 
the gauge fields on the boundary of this gauge system, 
the Aharonov-Bohm phase is well specified. This phase 
(and fields on the boundary) are unaltered by fluctua- 
tions at any place which does not involve gauge fields on 
the boundary of the system. By contrast, for a ferromag- 
net, fluctuations inside the system (local perturbations) 
do lift the degeneracy between different boundary condi- 
tions. Here, the boundary conditions and the local flelds 
are not independent of one another. 

To make this distinction more transparent, we may set 
= Xi in Eq. (j4|). Here, Xi is a pointer to a quasi- 
local configuration about site i. That is, Xi = 1 if t^a-t 
quasi-local configuration appears about site i and Xi = 
otherwise. Here, we will have by Bayes' theorem 

(Xi).-T(iH = ^^|^. (5) 

In Eq. P{^\ct) is the conditional probability of find- 
ing that local configuration about site i given that the 
boundary is of the type a. If various boundary configu- 
rations are equally probable - P{a) = P{f3) then Eq. (|4]) 
implies that, up to exponential corrections, 

T(ina) = T(in/3). (6) 

Equation ^ implies that the correlations P(i n a) be- 
tween the local configuration (i) and the boundary con- 
ditions (a) are the same for all boundary conditions. A 
sufficient condition for Eq. ^ to hold is that all correla- 
tions are sufficiently short ranged. Here, the only way to 



determine anything about the boundary is to perform 
a non-local measurement. Degeneracy between differ- 
ent topological sectors or boundary conditions cannot be 
lifted by local perturbations. The system is robust to lo- 
cal fluctuations. If we specify some topological condition 
corresponding to a large set of possible boundary con- 
ditions then as for every unique boundary condition, we 
have {V)a = {V)i3 (up to exponential corrections). Such 
a topological sector can correspond to the total number of 
vortices in an XY system, Hopf invariants, total domain 
wall parity in an Ising chain, etc. 

In their disordered phase, systems often explicitly ex- 
hibits a condensate of the dual (disordering) fields. The 
dual fields (domain walls/ vortices/ • • • ) are topological in 
character - they are a product of an infinite number of 
local operators (i.e. of 0{L'^) local operators with d > 0) 
which contains fields all of the way up to the boundary 
of the system. Here, the disordering (defect) fields are 
given by operators of the form 

Tdefect = n 

with Or a quasi-local operator and Q a region which 
spans 0{L''-) sites/bonds, etc. 



IV. SYMMETRY, DEGENERACY, AND 
FRACTIONALIZATION 

As is well known, the presence of a symmetry allows 
for degeneracies yet does mandate their presence par- 
tout. For instance, the trivial global Z2 up/down symme- 
try of the Ising magnet is intimately tied to the twofold 
degeneracy of its model. More formally, if the group 
of symmetries of the Hamiltonian H \s Q — {fli}, i.e. 
[H,Qi\ = 0, and |^) is an eigenstate of H with eigen- 
value then Qi]"^) is also an eigenstate with the same 
eigenvalue. Symmetries do not mandate the existence of 
degeneracies at all energies. It may happen (as it often 
does) that a symmetry operator when acting on a certain 
state does not alter it at all: the state transforms as a 
singlet under that symmetry. 

What is just as obvious yet not as often alluded to is 
that degeneracies (of any energy level) always imply the 
existence of symmetries (and indeed not the converse). 
The proof of this assertion is trivial. In what follows, we 
focus on the GS manifold. By a simple relabeling, the 
same argument applied to any other energy level. 

As the lowest energy subspacc of H will play a key 
role in the following note, let us assume that the GS 
subspace, Mo = {|.9q)}, a £ Sq = [l,Ng], is degenerate 
with (galgp) = Sai3 for a,/3 G So- Let us now prove 
that this implies the existence of an, at most, SU{Ng) 
symmetry. Let us arrange the eigenstates of H such that 
the first Ng states are the degenerate states in question 
{Iffa}}- The remaining eigenvectors, A4o = {Iv/s)} with 
(3 £ iSjf = Ng -|- 1, • • • , form the orthogonal complement 
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with {va\vi3) = Sap for a,P E S^. Next, we construct a 
matrix A such that 

{9a\A\gp) = {ga\U^\gf3) (if a,/3£So), 

(ga.\A\vp) = (if a e 5o and /3 e S^), (8) 

with Ui any SU{Ng) operator. It is straightforward to 
prove that ^ is a unitary matrix, AA'' = l. More- 
over, in this eigenbasis, we see that [if, ^] = 0. Thus, 
H' = A'' HA commutes with H and shares the same set 
of eigenvalues. This iUustrates that A reahzes an SU (Ng) 
symmetry [l^ of the GS manifold. 

For pedagogical purposes and clarity, we revisit the 
link between symmetry groups (in our case GS sym- 
metry groups) and fractionalization. Although some of 
these ideas have been around, we have not seen the di- 
rect derivation of our central and general approach to 
fractionalization: deducing its existence by an analysis 
of the unitary group which resides on the GS manifold 
which we just derived. Most treatments to date either 
focused on the unitary group of the fundamental inter- 
actions or on counting arguments of GSs or associated 
zero modes. Here, we will show how the unitary group 
acting within the GS manifold directly unifies the two 
treatments in leads to fractionalization. 

The notion of fractionalization is defined relative to the 
original degrees of freedom of the system. For Ng degen- 
erate GSs, we saw earlier that an SU{Ng) group residing 
within the internal space of the degenerate manifold is a 
symmetry of the Hamiltonian, [H, U] = 0. In many cases, 
the SU{Ng) symmetry may splinter into small subgroups 
which suffice to characterize physical aspects of the GS 
manifold. Fractional charges are associated with the cen- 
ter of such symmetry groups. The center of the group 
is enough to characterize the coupling of the elements 
of the group to the U{1) gauge potential (the electric 
charge). The center of the group SU{Ng) is the mod- 
ular additive group ZjVg which can be associated with 
topological homotopy groups. That Z^r is the center 
is seen by noting that if any clement of the group com- 
mutes with all other elements then it must be a multiple 
of the identity. Any such unitary matrix, having equal 
diagonal only elements must have equal pure phases (z) 
along the diagonal. For the matrices to have unit deter- 
minant, = 1. This, in turn, immediately leads us to 
see that the center of the group is . If we complete a 
2tt revolution, then the Aharonov-Bohm phase for a in- 
teger quantized charge would lead to no change (the ex- 
ponentiated phase factor would be 1). For an SU{2) rep- 
resentation however, for half-integer spins, the elements 
of SU{2) suffer a factor of (—1). This may be seen as 
an Aharonov-Bohm phase of a charge q = 1/2. This 
leads us to see that in the fundamental representation 
S = 1/2 and that a general spin 5* can only be integer 
multiples thereof. Similarly, for the lowest rank faith- 
ful representation of SU{3), the phases are e^'^*/'^. The 
latter may be similarly viewed as the Aharonov-Bohm 



phase of a charge 1/3 object when taken around a quan- 
tized monopole. In quantum chromodynamics (governed 
by SU{3)), the fractional 1/N charge adduced to SU{N) 
goes by the name of N-ality. This is related to the sim- 
pler fractional charge as deduced from the Su-Schrieffer 
counting arguments in polyacetylene [l^ . [Tsj . A related 
approach was recently highlighted by [5| . It is important 
to stress that other (non-fundamental) charges may ap- 
pear: our discussion above only leads to a lower bound for 
the simplest always present SU{Ng) symmetry. Further- 
more, in several non-Abelian systems (various Quantum 
Hall states and vortices in p+ip superconductors) [l^ 
there are symmetry operators which link GSs of systems 
which are related by a gauge symmetry - there the con- 
siderations are richer than in the usual iV-ality counting 
arguments. 

Degeneracies imply the existence of symmetries which 
effect general unitary transformations within the degen- 
erate manifold and act as the identity operator outside 
it. In this way, GS degeneracies act as exact GLSs. The 
existence of such unitary symmetry operators (generally 
a subset of SU{N)) allows for fractional charge (A^-ality 
in the SU (3) terminology of quantum chromodynamics) 
implying that degeneracy allows for fractionalization de- 
fined by the center of the symmetry group. The (m- 
)rized Peierls chains [3j[l3l constitute a typical example 
of a system with universal (m— independent) symmetry 
operators, where fractional charge quantized in units of 
e* = e/m with e the electronic charge is known to occur 
In the Appendix, we will provide explicit expres- 
sions for the symmetry operations which link all of the 
GSs to one another. Different Peierls chain GSs break 
discrete symmetries in this D ^ 1 system: this leads to a 
violation of Eq. ([T]) in this system with fractionalization. 
The Fermi number Nf in the Peierls chain and related 
Dirac-like theories is an integral over spectral functions 
[m ; the fractional portion of Nj stems from soliton con- 
tributions and is invariant under local background defor- 
mations psj . 

V. GAUGE-LIKE SYMMETRY: A DEFINITION 

A d-GLS of a theory given by H (or action S) is a group 
of symmetry transformations Gd such that the minimal 
non-empty set of fields (/){ changed by the group opera- 
tions spans a d-dimensional subset C C A. Here, A de- 
notes the entire lattice (or space) on which the theory is 
defined. These transformations can be expressed as [l7j : 

= n gik, (9) 

ieCi 

where C\ denotes the subregion 1, and A — (JjCi. (The 
extension of this definition to the continuum is straight- 
forward.) Gauge (local) symmetries correspond to d = 0, 
while in global symmetries the region influenced by the 
symmetry operation is d ~ _D-dimensional. These sym- 
metries may be Abelian or non-Abelian. In the next two 
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Sections (Sections I VII and IVII[) . we analyze the d-GLSs 
present in various systems in great detail. These Sec- 
tions are lengthy and aimed at being pedagogical. In 
Section IVIl we review the prominent examples of TQO 
and identify the d-GLSs present in each of these sys- 
tems. In Section fVlIi we discuss other well known phys- 
ical systems which exhibit d-GLSs and arc candidates 
for TQO. After this long list of examples in which d- 
GLSs are present we turn in the following Sections to the 
consequences of d-GLSs. After a list of properties relat- 
ing to low-dimensional behavior (for low d) enforced in 
high-dimensional systems (Section IIX[) . we turn to our 
central result concerning the link between TQO and d- 
GLSs. In particular, in Section [Xj we show how exact 
low-dimensional d-GLSs mandate TQO. 



VI. KNOWN EXAMPLES OF TQO AND THEIR 
GAUGE-LIKE SYMMETRIES 

We now list the most prominent examples of TQO and 
identify the d-GLSs present in each of these systems. 
These symmetries either appear exactly in the system 
or appear in some limit in a system which is adiabati- 
cally connected to the system in question. Symmetries 
do not need to be exact in order to exert their influence: 
so long as points in parameter space may be adiabatically 
linked to each other, they lie in the same phase associated 
with the same set of symmetries. At issue is whether or 
not such a smooth adiabatic connection exists. In what 
follows, we start with a list of systems in which the d- 
GLSs are exact (examples (a)-(c)). We then turn to the 
Fractional Quantum Hall Effect (FQHE) (in (d)) where a 
clever choice of parameters adiabatically links the D ~ 2 
system to a D = 1 in which a discrete d = 1 symmetry is 
present. It is important to stress that in any Quantum 
Hall system, the two-dimensional magnetic translation 
group is an exact symmetry which links all GSs. How- 
ever, a more transparent understanding is attained by 
considering the one-dimensional limit. In (e), we com- 
ment on the similarities between the symmetries found 
in the suggested non-Abelian states at the core of vor- 
tices in (p+ip) superconductors to the symmetries of the 
FQHE problem. In both (d) and (e), the physical pic- 
ture behind a one-dimensional {d ~ 1) symmetry which 
links different GSs (an evolution around a toric cycle) is 
intimately linked to suggestions regarding the use of such 
braiding operations in quantum computing. 

(a) In Kitaev's Toric code model on the square lattice, 
the Hamiltonian 

Hk = -J2As-J2Bp (10) 

s p 

where 

^= n n <- (11) 

ij^5tar{s) ijGboundary(p) 



An illustration is provided in the bottom panel of Fig. [2l 
The product defining As spans all bonds ij which have 
site s as an cndpoint (see the cross-shaped object in Fig. 
[5]). The plaquette product Bp spans all bonds which lie in 
the plaquette p (see the plaquette in the bottom panel of 
Fig. [2]). In the presence of periodic boundary conditions, 
this model possesses two d = 1 Z2 symmetries which are 
given by [2| 

These d ~ 1 symmetries are denoted in the lower panel 
of Fig. [2] by the dark solid and dashed lines. Along the 
solid lines a product of the a^j fields is taken while along 
the dashed lines a product of the afj fields is performed. 
The contour Ci represents the solid horizontal line, C2 
denotes the solid vertical line, C[ is the dashed vertical 
line, and is the dashed horizontal line. That these are, 
in our definition of Eq. d — 1 Z2 symmetries follows 
as the regions {Ci} are lines (d = 1) on which Z2 gener- 
ators {<Ji~^'^) are placed. The topological symmetries of 
Eq. are none other than discrete d = 1 symmetries 
in the general classification of GLSs. 

(b) Gauge theories. The Kitaev model of (a) may be 
viewed as a Z2 gauge theory in which usual local gauge 
symmetries [ISj have been removed by hand, leaving bare 
only the two additional d ~ 1 symmetries that a Z2 gauge 
theory exhibits. On a square lattice, the Z2 gauge theory 
Hamiltonian reads 

i/z, =-ii^a^a|,a^;a£-/.,^a^. (13) 

□ >■] 

Here, aside from the d = 1 symmetries Xi 2 of Eq. (|12p . 
there are also the standard local (d = 0) gauge symme- 
tries {^is}. This set of low-dimensional (d = 0, 1) GLSs 
suffices to link all GSs to each other. As a consequence, 
this gauge theory is an example of TQO. (Indeed, gauge 
theories form the most prominent examples of TQO.) 

The middle panel Fig. [2] provides a cartoon of the in- 
teractions of Eq. dini) - the bright (red) plaquette terms 
and crosses (on site magnetic field) as well as the sym- 
metries of the system (the dark dashed lines representing 
the the d = 1 symmetries Xi^2 and the dark solid cross 
representing the d = symmetries captured by {As} of 

Eq- CD)- . _ 

(c) Wen's plaquette model [l9| is given by 

Hw = -KS^afa^^^ erf. - , - cr^^- , (14) 

i 

with all fields lying on the vertices of a square lattice. 
This system displays d~ \ symmetries whose generators 

Op = n <■ (15) 

iG P 

With the line P chosen to be any horizontal or vertical 
line of sites we find that Op\ — 0. Similar symmetry 



operators can be extracted for the other known exam- 
ples of TQO. An illustration of the interactions and 
the symmetries is provided by the top panel of Fig. [21 
Here, the solid dark lines denote the symmetries of 
Eq. ([15]). The plaquette in Fig. [2| denotes the interaction 
terms in Eq. ([T4|) . As in Kitaev's Toric code model (of 
(a)) where all operators appearing in the summand of 
the Hamiltonian commute, all of the operators in the 
summand of Eq. (jl4p commute with one another. We 
will later show (Section IXIII B[) that Wen's model is 
identically the same as Kitaev's model when written in a 
rotated basis. The gauge theory of (b) in the absence of 
an applied field [Eq. ([T3|) with = 0] also has an energy 
spectrum which is identical to that of the Kitaev's Toric 
model (of (a)) or Wen's plaquette mode (of (c)). As we 
will show in Section IXIIIl all of these systems have an 
energy spectrum which is identical to that of the Ising 
chain. Schematically, we write 

Kitaev's model = Wen's model ^ D=2 Ising gauge. 



(d) Fractional Quantum Hall States 

Fractional quantum Hall (FQH) states have long been 
regarded as the quintessential example of TQO. These 
states are realized in D = 2 electronic systems subjected 
to a strong magnetic field B applied perpendicular to 
the plane where the electrons reside. A minimal model 
Hamiltonian (for electrons in a surface of size [Lx x 
Ly]) that captures the FQH physics is 





_rt o <^__^_<^_ 



(16) 



with canonical momentum Hi ^ pi — ^Ai, Ai = A{ri) 
the vector potential corresponding to a uniform magnetic 
field S = V A A, and V{\ri — fj |) a two-body (Coulomb) 
interaction. 

This system has a magnetic translation group (re- 
viewed in the Appendix) which resides in d = 2 dimen- 
sions. It has long been recognized, however, [2^ that 
the symmetry operators linking the different GSs also 
correspond to operators which move a quasi-particle and 
quasi-hole along a d = 1 trajectory - a symmetry of a 
d = 1 character in the scheme of Eq. ([9]). This symmetry 
relates to the magnetic translation group. Recent work 
[2ll . [2^ allows us to make this intuition more precise 
and cast these symmetries as explicit d = 1 symmetry 
operators which nicely fit in our general scheme of clas- 
sifying the symmetries which enable TQO. As shown in 
[21I . [23 ] , the Quantum Hall problem m. D = 2 dimensions 
is adiabatically connected io & D = \ Peierls-type prob- 
lem where the magnetic translation symmetries which we 
briefly review in the Appendix become explicitly discrete 
symmetries of dimension d= 1 [2l|, . 

This dimensional reduction makes explicit the one- 
dimensional (c? = 1) character of the magnetic translation 




_0 o 0--J--O- 



FIG. 2: Schematics of the interactions and symmetries in- 
volved in three of the known examples which display TQO. 
The lower panel represents Kitaev's Toric-code model, the 
middle one corresponds to the Z2 gauge theory and the upper 
panel corresponds to Wen's plaquette model. Dark lines (solid 
or dashed) represent symmetries while the brighter (red) lines 
or crosses correspond to the interaction terms. The hollow 
circles represent the spin locations. See text. 
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group operator which is responsible not only for topolog- 
ical order but also re-explains in greater depth why the 
degeneracy of a Quantum Hall system scales in the way 
that it does, see e.g. [l^, [2^. This relies on the cor- 
respondence with the discrete d = 1 symmetries of the 
Peierls chain to which we briefly alluded to in the Section 
llVl and which we explore in detail in the Appendix. In an 
Abelian Quantum Hall system of filling fraction v = p/q 
with p and q relatively prime integers which is embedded 
on a manifold of genus g (a manifold which has g handles) 
the GS degeneracy scales as q^ as would be expected for 
a system harboring q independent discrete d = I transla- 
tion symmetries. This is the degeneracy of an open chain 
of degeneracy q when punctured by — 1) holes. 

By contrast, the degeneracy for a system with two in- 
dependent discrete symmetries in two dimensions would 
have a degeneracy which is expected to scale as q^ x q^ . 
Similar considerations, albeit slightly more sophisticated, 
apply in the non- Abelian arena [TgI. [2TI. [2^ . 

(e) Majorana Fermions in Vortex cores of p+ip super- 
conductors 

The d = 1 symmetry operators of an evolution of a 
quasi-particle around a toric cycle is very similar to that 
appearing in the Quantum Hall problem (l6j . The con- 
siderations outlined above can be replicated for these sys- 
tems. These d = I (braiding) symmetries in non-Abelian 
states of FQHE and in (p-|-ip) have received much atten- 
tion because of their viable use for quantum computation. 



VII. OTHER SYSTEMS WITH GAUGE-LIKE 
SYMMETRIES 

We now present several other systems which have low- 
dimensional GLSs. Following the proofs to be provided 
in Section IXll many of these systems exhibit TQO. 

In the systems of (a (i)) and (d), a low energy emer- 
gent lower-dimensional GLS appears which is not suffi- 
cient in leading to finite-T TQO. In all other systems 
introduced in this Section, the d-GLSs which appear are 
exact and hold over the entire spectrum. In what follows, 
we will discuss whether or not these systems may exhibit 
broken symmetries - our conclusions rely on the theorem 
of (l7j which is stated anew in Appendix [B) As we will 
explicitly show in Section [Xl if all GSs can be linked by 
the exclusive use of low d-GLSs {d < 2 continuous sym- 
metries or c? < 1 discrete symmetries) then finite-T TQO 
may follow . 

a) Orbitals - In transition metal (TM) systems on cubic 
lattices, each TM atom is surrounded by an octahedral 
cage of oxygens. Crystal fields lift the degeneracy of the 
five 3d orbitals of the TM to two higher energy Cg lev- 
els {{d^z^-r'^) and \dx2_y2)) and to three lower energy t2g 
levels {\dxy), \dxz), and \dyz)). A super-exchange cal- 
culation leads to the Kugel-Khomskii (KK) Hamiltonian 



MM 

H=Y.Krb(Sr-Sr- + \). (17) 

(r\r') 

Here, Sr denotes the spin of the electron at site r and 

^orb are operators acting on the orbital degrees of free- 
dom. For TM-atoms arranged in a cubic lattice, 

= JiiKK' - - + 1), (18) 

where tt" are orbital pseudo-spins and a ~ x,y, z is the 
direction of the bond (r, r'). 

(i) In the Cg compounds, 

n^-y = i(-< ± a/30, n^. = (19) 

This also defines the orbital only "120°-Hamiltonian" 
which is given by 

Horb^jY. E ^?^r+e^- (20) 
r.r' a—x.y,z 

Jahn- Teller effects in Cg compounds also lead, on their 
own, to orbital interactions of the 120°-type [2y|. The 
"120° model" model of Eqs. and ^ displays 

discrete {d = 2) [22]"^^ GLSs (corresponding to planar 
Rubik's-cube-like reflections about internal spin direc- 
tions - Fig. [T]). Here, there is only an emergent symmetry 
in the GS sector of the classical {S 00) rendition of 
the problem. In this rendition the spin-1/2 operators a 
arc replaced by spin-S generators. These symmetry op- 
erators O" are [13, [H, 

rePc 

These operators correspond to a rotation by 180 degrees 
about the directions (— ^,±^) and (0,1) in the inter- 
nal xz spin plane. Here, a = x,y, z and may denote 
any plane orthogonal to the cubic e-a axis. These are dis- 
crete d — 2 symmetries which emerge only within the low 
energy sector and may be broken at finite temperature. 
Within the large S limit, and at finite temperatures, the 
discrete d = 2 Z2 symmetry is broken by entropic fluctu- 
ations [2l,[2i|. 

(ii) In the t2g compounds (e.g., LaTiOa), we have in 
Horb of Eq. (HOD 

K = l<^r- (22) 

This is called the orbital compass model [H, [l^, H^, [27l . 
[28j . The symmetries of this Hamiltonian are given by 
Eqs. (EI]), and (EH). In the D = 3 model of Eq. 
rotations of individual lower-dimensional planes about an 
axis orthogonal to them leave the system invariant [see 
Fig. [3] . It should be noted that here Berry phase terms 
lift the d = 1 Z2 symmetries present in the (large S) 
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classical orbital compass model (the b spin component 
Sb — *■ (1 — 2Sab)Sb along a ray parallel to the spatial 
a— th axis) to a weaker d = 2 symmetry (whose lattice 
version is given in Eq. (HJ)). Insofar as their algebra 
relations are concerned, in 13 = 2 variants of the orbital 
compass model of Eq. (j20p (in which a ~ x,z), where 
the planes P become one-dimensional lines parallel to the 
coordinate axes, we can regard a.s S — 1/2 operators. 
For example, for a certain given plane P orthogonal to 
the X axis, is set to be and for a plane orthogonal to 
the z— axis, we set to be r^. The operators {tj.,Tz} 
satisfy a spin S* = 1/2 SU{2) algebra. In this regard, 
these operators can be viewed as Abelian anyons. Here, 



J^y—JD J^^-Ji) j^y-^ J?^^^ 



O^Oy = OyO^e' 



(23) 



with 6 = 7r/4 (semions). A correspondence between the 
statistics (9) and toric {d = I) translation operators in, 
for example the Quantum Hall problem is reviewed in 

The set of c? = 1 symmetries discussed above does not 
suffice to link all of the GSs. In the d — 2 orbital com- 
pass model, an additional global discrete symmetry is 
required. The D — 2 system displays an additional in- 
dependent Z2 reflection symmetry {ax CTx, CTz — *■ ctx) - 
a rotation by tt about the symmetric line (the 45 degree 
line in the plane), i.e. 



Reflection 



]Jcxp[i^!^(af 



FIG. 3: From Refs. [23, [2J|. The symmetries of Eq. (|2T)) 
applied on the uniform state (at left). 



of sites i on which af = 1 is even on every row and on 
every column is invariant under the d = \ symmetries of 
the D = 2 orbital compass model. 

To make contact with the previously examined models, 
we remark that upon performing standard D = 1 spin 
duality transformations along the diagonals (see [28j). 
we can see that the transverse field Z2 gauge theory of 
Eq. (jl3p is dual to an orbital compass model in which 
half of the interactions corresponding to one polarization 
direction of the spins arc missing. Alternatively, if we 
define, by reference to Eq. (fH 



with 



(24) 



(27) 



(28) 



This symmetry (ORefiection ) is a manifestation of a self- 
duality present in the model. The 13 = 3 orbital com- 
pass model displays similar reflection symmetries (per- 
mutations). As discrete d = 2 Ising symmetries can be 
broken, the D ~ 2 orbital compass systems display finite- 
T nematic orders [H, [li, [13 , 



y V 



(25) 



The additional reflection symmetry of Eq. ([M)) . is needed 
in addition the d = \ discussed earlier, to link all GSs 
to each other. As we will later discuss, such a high- 
dimensional (d = 2) Ising symmetry can be broken (and 
indeed is) [23, [l^, [23] . Such a broken symmetry docs not 
make this system a candidate for TQO. This symmetry 
is removed in strained variants of the orbital compass 
model of Eq. ([^0]) , which are given by 



E 



(26) 



r.r a—x^y.z 



with all {Ja} different. It is readily verified that no > 
1 symmetries exist for the system of Eq. (PB)) . Orbital 
systems approximated by Eq. (|26p are candidates for 
TQO. It should be noted that these symmetries naturally 
allow for a high degeneracy of most levels, yet they do 
not mandate it. For instance, a state which is an equal 
amplitude superposition of all states in which the number 



a star term of the same geometry as that of Ag in Eq. pT|) 
then 



Spec{iJc} = Spec{Horb[Jx = K, Jy ^ h]}. 



(29) 



In Eq. (|29p . Spec denotes the energy spectrum of the 
theory. 

In the classical (large S) limit of the orbital compass 
model, the symmetries discussed above mandate that 
each state is, at least, 2^ degenera te psi . [2^ . As noted 
by [2^ and nicely expanded on by |30l|. in the quantum 
problem, the non-commuting nature of the symmetry op- 
erators Op^ and for two orthogonal planes Pi and P2 
mandates that the GS of the D = 2-dimensional orbital 
compass model is, at least, two fold degenerate. That 
this lower bound is indeed realized was verified numer- 
ically [2^, [131 ■ We find that this lower bound on the 
degeneracy follows, for odd size lattices, even more gen- 
erally as a direct consequence of Kramers' degeneracy 
for the time reversal invariant Hamiltonian of the orbital 
compass model. We further remark here that on a cubic 
lattice {D = 3) of size L x L x L where L is even, all 
symmetry operators of Eqs. (PT|) . and (|22p commute and 
in principle, a single non-degenerate GS may be realized. 
This state can be separated by exponentially small en- 
ergy gaps from all other states in the same topological 
sector. If, at least one of the sides of the lattice (L^.y^z) 
is odd then the GSs are, at least, two fold degenerate. 
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FIG. 4: Schematics of the interactions and symmetries in- 
volved in the classical rendition of three D = 2 examples. 
The left panel represents (i) the local (d = 0) symmetries of 
the Ising gauge theory [Eq. psp ]. The middle panel repre- 
sents (ii) an orbital compass model [Eqs. (|20|) . and (|22|l ] with 
d — 1 symmetries; here the symmetry operations span lines 
[see Eq. (|21|l . The right panel depicts (iii) an XY model with 
d = 2 symmetries; the symmetry here spans the entire D = 2 
dimensional plane. 



A schematic of the c? = 1 symmetries of the orbital 
compass model and its comparison to local and global 
symmetries is provided in Fig. 3] 

b) Spins in transition metal compounds- Following [3l[ , 
we label the three t2g states |dyz), \dxz}, \ dxy) by jF), 
and \Z). In the t2g compounds, hopping is prohibited 
via intermediate oxygen p orbitals between any two elec- 
tronic states of orbital flavor a {a — X, Y, or Z) along 
the a axis of the cubic lattice (see Fig. [5]). As a conse- 
quence, as noted in [3l|, a uniform rotation of all spins, 
whose electronic orbital state is |a), in any given plane 

(P) orthogonal to the a axis c|^^ = Y^n^^^'n Aai-, with 
cr, 7? the spin directions, leaves Eq. p7|) invariant. The 
net spin of the electrons of orbital flavor |a) in any plane 
orthogonal to the cubic a axis is conserved; this consti- 
tutes the conserved topological charge in this case. Here, 
we have d = 2 SU{2) symmetries which are given by (3l| 



dp,a. = [exp{iSf,-9'^)/h], [H,dp,a.]=0, (30) 



with 



(31) 



the sum of all the spins 5*'" in the orbital state a in any 
plane P orthogonal to the direction a (see Fig. [5|). The 
spin of Eq. (pij) constitutes a topological charge. There 
are (3L) conserved topological charges. There is a spin 
continuity equation in every plane conjugate to each of 
the conserved total planar spin. 

Similar to the orbital compass model, an additional 
reflection symmetry is necessary to link all GSs. The 
discrete high-dimensional reflection symmetry may be 
broken to lead to ncmatic type order. The gauge-like 
invariance (Eqs. and ([50]) ) of nematic type order 

parameters with respect to the d ~ 2 SU (2) symmetries 
allows spin nematic order. If the KK Hamiltonian is the 
most dominant spin exchange process, one may predict 
a spin nematic order which onsets at temperatures much 



higher than the currently measured gauge non-invariant 
magnetization [l3l . In the strained orbital system of Eq. 
(PS)) . the discrete high-dimensional reflection symmetry is 
removed. The remaining d = 2 SU(2) symmetries allow 
for TQO. 




FIG. 5: From Ref. [17t |. The anisotropic hopping amplitudes 
leading to the Kugel-Khomskii (KK) Hamiltonian. Similar to 
[sH , the four lobed states denote the 3d orbitals of a transition 
metal while the intermediate small p orbitals are oxygen or- 
bital through which the super-exchange process occurs. The 
dark and bright shades denote positive and negative regions of 
the orbital wave-function. Due to orthogonality with interme- 
diate oxygen p states, in any orbital state \a) (e.g. \Z) = \dxy) 
above), hopping is forbidden between sites separated along 
the cubic a [z above) axis. The ensuing super-exchange (KK) 
Hamiltonian exhibits a d = 2 SU{2) symmetry corresponding 
to a uniform rotation of all spins whose orbital state is \ct) in 
any plane orthogonal to the cubic direction a. 



AD- 
by Eqs. 



2 rendition of the t2g KK Hamiltonian is defined 



18]), and ((221) the xz plane with only two of 



the three t2g orbital flavors (a ~ x, z). In this system, the 



spin conservation of Eq. (|3ip holds for each individual 
line (a continuous d = 1 symmetry). 

c) Superconducting arrays: A Hamiltonian for super- 
conducting {p + ip) grains (e.g. of Sr2Ru04) on a square 
grid, was recently proposed, [33] 



(32) 



Here, the four spin product is the product of all spins 
common to a given plaquette □. The spins reside on 
the vertices on the plaquette (not on its bonds as gauge 
fields). These systems have {d = 1 Z2) symmetries sim- 
ilar to those of the D = 2 orbital compass model. With 
P any row or column. 



0. 



(33) 



In fact, the model of Eq. (|32|) can be shown to not 
only have the same symmetries but to also have the 
same spectrum as the orbital compass model. This is 
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shown by duality transformations j28| which link the 
two systems, 



(p+ip) model 



D=2 orbital compass model. 



There is a deep link between the system of Eq. 
and the Z2 gauge theory of Eq. ([13]). To see this link, we 
may express Eq. (fT^ in terms of sites which are located 
at the centers of these bonds. Rotating the lattice by 
45 degrees, we immediately obtain the Hamiltonian of 
Eq. ([32I) with one difference: the transformed Eq. (fT3|) 
will not correspond to the sum of the products of af^ 
around all plaquettes. Rather, only half of the plaquettes 
appear in the rotated Eq. (jl3[) . This smaller number 
of terms allows for a higher or (in the worst case) an 
equal number of symmetry operations which leave the 
Hamiltonian of Eq. (|13p invariant. In the aftermath, 
all of the d = 1 symmetries of Eq. ((32)) arc symmetries 
of the Z2 lattice gauge theory of Eq. (fT^ (as they must 
be). The local symmetries of the gauge theory (the d = 
operators {As} of Eq. (fTT|) ) are no longer symmetries of 
the system of the more constrained system of Eq. ([32]) . 
These residual d = 1 symmetries are none other than 
the topological symmetries of the I2 lattice gauge theory 
emphasized by Kitaev and Wen 0, [l^ . 

d) Klein models: Klein spin models [s^ have an ex- 
ceptionally high number of low-dimensional symmetries. 
As an example with the highest degree of symmetry, we 
quote the model on the D = i pyrochlore lattice and 
on its 13 = 2 (checkerboard) incarnation. Here, in both 
cases (the pyrochlore and checkerboard lattices), the cor- 
responding short-range Hamiltonian reads 



H = 



{ij)a 



(34) 



with K ~ 4 J/5 > and the projection operator 

onto the state [in each tctrahcdron/plaquettc a = ijkl for 
the pyrochlore/checkerboard lattices respectively] with 
maximal total spin {S — 2). All GSs in these systems 
are superpositions of singlet coverings [13, 



within,,) - -L(iu)-UT)), 



(35) 



and a^^ arbitrary amplitudes (up to an overall normal- 
ization). Here, the dimer coverings are labeled by W. 
The condition for each dimer covering is that a single 



dimer appears in each tetrahedral unit of the pyrochlore 
lattice. These states are separated from all other excited 
states by a finite gap [s^. Here, within the GS manifold 
the systems maps into the six vertex model (see Fig. [7|) , 
the number of dimer state GSs is bounded from below 
by (3/2)^-/2. (On the D = 2 version of the pyrochlore 
lattice - the checkerboard lattice - the dimer state GS 
degeneracy is exactly equal to (4/3)^^'/^.) As the expo- 
nential in volume degeneracy suggests, these states may 
be linked to each by local (d = 0) symmetry operations. 
An illustration of such a process is shown in Fig. [51 




FIG. 6: From Ref. [36(1. All GSs of the Klein model on the 
pyrochlore and checkerboard lattices are composed of dimer 
coverings. The ovals denote singlet dimer states. The arrows 
denote the representation of these dimer states within the 
six vertex model. On each plaquette (tetrahedron) the dimer 
connects the base of the two incoming arrows. The GSs here 
are highly regular. Many other GSs exist. The GSs exhibit 
an emergent local (d = 0) GLS. Local (d = 0) variations of 
dimer coverings product new GSs. 



X X X X X X 
X X X X X X 



FIG. 7: From Ref. [36|]. Standard representation of the six 
vertex states in terms of lines. Every line is composed of links 
whose arrows flow to the right in the vertex representation. 

We illustrated earlier that the existence of a GS de- 
generacy implies the existence of symmetry operators 
which perform unitary transformations within the GS ba- 
sis and which, otherwise, act trivially. The d — symme- 
tries within the GS manifold are exact symmetries of the 
Hamiltonian. Although (due to the d = symmetries 
that this model displays), the system does not possess 
local order, it does display TQO. Different sets of states 
in the GS manifold are characterized (in the line repre- 
sentation of Fig. [7]) by the net number of vertical lines 
on each horizontal row of lattices (and conversely). This 
number is conserved and corresponds to a topologically 
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FIG. 8: From Ref. Local GLSs in the GS manifold. 

Left: a representation of the the six vertex states in a highly 
regular configuration (a) in terms of lines [see Fig. [7] (c). 
Joining the vertex coverings according to the standard pre- 
scription allows on to label all GSs. Right: an elementary 
local (gauge-like) spin interaction process leading to a new 
dimer configuration in the GS basis (b). This corresponds in 
the line representation (d) to a flip of one of the line corners. 



invariant flux of lines. Here, these d = symmetries 
(and TQO) appears exa ctly at the onset of a dcconfined 
T = critical point ^1^. 

The GSs of the Klein model on a pyrochlore lattice can 
be linked by discrete d = 0,1 symmetry operations. The 
d = 1 symmetries link GSs differing by the net number of 
lines (the flux). We can remove the d = 1 symmetries by 
prescribing boundary conditions which fix the number of 
lines at a given value and leave only the d = symme- 
tries. These lead to TQO at T = 0. The GS of the Klein 
model on the checkerboard lattice map onto the six ver- 
tex model at its ice point. Here, the correlations between 
dimers (and finite sets of dimcrs) fall off algebraically (as 
in a dipole gas). 

e) Other systems: similar symmetries were found in 
frustrated spin systems. Ring exchange Bose metals, in 
the absence of nearest-neighbor boson hopping, exhibit 
d = 1 symmetries [s^ . Continuous sliding symmetries of 
Hamiltonians (actions) invariant under arbitrary defor- 
mations along a transverse direction, 

(j){x,y) ^ (p{x,y) + f{y), (36) 

appear in many systems. Amongst others, such systems 
were discovered in works on Quantum Hall liquid crys- 
talline phases [3^, HO], a number of models of lipid bi- 
laycrs with intercalated DNA strands [4l[, and sliding 
Luttinger liquids [4^ . 



VIII. GAUGE-LIKE SYMMETRIES AND 
AHARONOV-BOHM-TYPE UNITARY 
TRANSFORMATIONS 

In this Section we formally write down our group 
symmetry elements and topological defects in terms of 
Aharonov-Bohm-type unitary operators. In the contin- 
uum limit the group elements of d = 1 GLSs can be 
written as a path-ordered (V) product 

U ^-pe^fcA-ds^ (37) 

where C is a closed path in configuration space and A 
is the corresponding connection. On the other hand, a 
defect creating operator in such a system is 

f+ = Ve^^+^"^\ (38) 

where (7+ is an open contour that only spans a portion 
of the entire closed cycle. On a discrete lattice these 
expressions are replaced by equivalent discrete sums, and 
C (C-I-) represents a closed (open) path on the lattice. 
For instance, for the orbital compass model [Eqs. ((^ . 
and ([211)], the d = 1 group elements can be written in a 
form which formally looks like an Aharonov-Bohm (AB) 
phase ^3 [Eq- (PT|) ]. For the same model, the defect 
depicted in Fig. [S] is generated by 

f+ = e'^ ^^■'^^+ (39) 

The unitary operators of Eqs. (pij) . and ([57)1 belong to 
a unitary group different from U{\). As such, Eqs. (|2ip . 
and ((57|) extend the AB phases associated with U{1) gen- 
erators. We can interpret T+ as the creation of a defect- 
antidefect pair and the displacement of each member of 
the pair to the opposite endpoints of C+. The d — 1 
GLS operators linking different GSs (the generators of 
Eqs. ((2T|) . and (|37)) ) correspond to the displacement of 
a defect-antidefect pair along a toric cycle. Formally, 
this is similar to the quasiparticle-quasihole pair opera- 
tor linking different GSs in the FQHE. Monodromy and 
toric topology are manifest in the d — 1 symmetry op- 
erators. The generator A (i.e. the effective connection 
A that appears in the integral for the d = \ symmetry 
operators) is borne by parallel displacement on spin (or 
other) fiber bundles. More complicated topological prop- 
erties are associated with d = 2 symmetries. 

If A is also a generator of a local symmetry (i.e. a 
c? = GLS) then by local symmetry operations we can 
deform the contour C continuously. Here, the group el- 
ement of Eq. p7)) will only depend on the topology of 
the configuration space and C. This is a particular case 
of a more general set of GLS generators that we discuss 
in the present paper. 
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FIG. 9: A soliton (topological defect) in an otherwise uniform 
state of a general (anisotropic) orbital compass model [see Eq. 
(|26|) ] leads to only a finite energy cost. A single energetic bond 
(dashed line) indicates that penalty. The energy-entropy bal- 
ance associated with such d — 1 Ising type domain walls is the 
same as that in a D = 1 Ising system. At all finite tempera- 
tures, entropic contributions overwhelm energy penalties and 
no local order is possible. Order is only manifest in non-local 
quantities associated with topological defects. Similar results 
occur in other systems with low-dimensional GLSs. 



IX. CONSEQUENCES OF GAUGE-LIKE 
SYMMETRIES 



with 

Wick,) = uiY^ 10^12)2 _ 1 ^^2(|^^|2) (41) 

and /i, 1/ = 1, 2, 3. The Lagrangian density C displays the 
continuous d = 1 symmetries 

0p ^ e'M-.}.^.)^^, (42) 

The conserved d = 1 Noether currents are tensors given 
by 

it^i, =i[(t)*,dy(j)t,- {dy(t)l)(j)^,], (43) 

which satisfy d ~ \ conservation laws 

[^u3^.. + drj^^r] = (44) 

(with no summation over repeated indices implicit). 
What is special about d-GLSs is that there is a conser- 
vation law for each line associated with a fixed value of 
all coordinates x^^^ relating to the d — 1 charge 

Qt,{{Xv^^j}) = / dx^, i^,r{x). (45) 



What are the physical consequences of having a sys- 
tem endowed with a symmetry group Qd! In the follow- 
ing Sections, we outline a few of those. First, we remark 
on the conserved Noether charges which are associated 
with these symmetries. We then give several toy exam- 
ples in which topological indices which are usually as- 
sociated with low-dimensional systems make an appear- 
ance in high-dimensional systems. Next, we discuss a 
property which will be of key interest to our forthcoming 
discussion and proof of finite-T TQO in several systems 
- that regarding the absence of spontaneous symmetry 
breaking (SSB). Then, we discuss an extension of that 
result which shows that well-defined quasiparticlc poles 
are precluded in systems which exhibit GLSs of low di- 
mensionality. We will conclude this Section with a brief 
discussion regarding the appearance of lower-dimensional 
topological terms when these symmetries are present. 



A. Conserved Noether charges 

Symmetries generally imply the existence of conserva- 
tion laws and topological charges with associated conti- 
nuity equations. We find that systems with d-GLSs lead 
to conservation laws within d-dimcnsional regions. To il- 
lustrate, consider the (continuum) Euclidean Lagrangian 
density of a complex field 4>{x) ~ {(f)i{x) , 4>2{x) ^ 4'^{x)) 

{X = {xi,X2,X3,)) : 

I^m'/'mP + ^ E I^-'^mP + (40) 



B. Topological indices 

Noether 's theorem stems from the continuous symme- 
tries of the Lagrangian. Its application to the toy system 
of Eq. (gO]) with the d = I symmetries of Eq. (gj) led 
to the conserved d = 1 charges of Eq. (|l5)) . By con- 
trast, topological numbers in field theories are related to 
boundary conditions. Their presence is dictated by the 
requirement of a finite energy (or action). 

In what follows, we will show how in specific continuum 
field theories, topological numbers which are associated 
with low-dimensional systems (of dimension d) can ap- 
pear in high-dimensional systems (of dimension D > d). 

(i) d = 1 Z2 soliton numbers: 

To illustrate the basic premise, consider the theory of 
Eq. (PP]) on a (£) = 2) square lattice for a system with two 

component fields (f> = {(f)i, 4)2) and in which the potential 
term is given by an extension to 2-|-l of the example in 
Section 2.3 of ^ 

+ \^t + \v<l^l{4>\^l), (46) 

with X,ri,m^ > 0. The system given by Eqs. (j40|l . and 
(Uni) displays d ~ 1 Z2 symmetries: these are given by 
(01 — > —4>i) on each horizontal line and (02 — *■ ^^2) on 
each vertical line. The minima of W occur at 

(01,02) = (±1,0), (47) 
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at which W = 0. For A > (77 — m^)^, the functional 
W is positive elsewhere. The requirement of finite en- 
ergy dictates that at infinity, the field must be of the 
form of Eq. ([T7|) . The form of W enables us to define a 
topological charge as that of soliton number in a = 1- 
dimensional system, 

Qsolitonix2) = [(t>l{xi = 00,X2) - (/ilixi = -OO, ^2 )] .(48) 

By virtue of the gradient terms in Eq. (|40p . for finite 
energy penalties, we must have a constant topological 
charge: Qsoiiton{x2) = QsoUton for aU X2. 

Various contending soliton configurations are given as 
a direct product of two D ~ 1 systems (one for (f>i along 
the Xi axis and one for (/)2 along the X2 axis). A solution 
is given by 

(pi{xi,X2) = tanh[m(a;i — c)]. 



4'2{xi,X2) = ±\ sech[m(x2 - c)], (49) 

V V 

with c an arbitrary real constant. The main feature of 
Eq. (|49p is the existence of an effective d = 1-dimensional 
soliton configuration for the field component Such a 
soliton-likc configuration is expected for D ~ 1 Z2 sys- 
tems. 

(ii) Pontryagin indices for d = 2 SU{2) GLSs: 
Next, we consider another toy model in 3+1 dimen- 
sions for which similar topological indices appear for con- 
tinuous d = 2 5*0(3) GLSs. Consider, at every lattice site 
a field composed of an 0(3) vector triad 

$ix-) = i$,ix),Mx),h{x)). (50) 

Here, each (/f^ {i = 1,2,3) is a vector with three com- 
ponents [and thus (/) contains nine components]. In the 
"hard spin" limit of for each i = 1,2,3, we will de- 
mand the normalization condition 

$,{x:) ■ $,{x) = 1. (51) 

For a general Lagrangian density of the form 

M2=0,l,3 

+ \ E (5^303) •(5'^^03) 
^3=0,1,2 

+W, (52) 

with W a functional of {(piix)} which is invariant under 
all global rotations of each (j>i{x): For example, for a 
soft spin implementation of Eq. ([^ . we may set W = 

J2^=ii^i{x) — 1)^. Here, each field exhibits a d = 2 
S0{3) rotation symmetry in each spatial plane [jk] for 



which i ^ j and i k (j ^ k). The static finite energy 
solutions to the Lagrangian density of Eq. ((5^ . exhibit 
three topological indices - one for each of the three fields 
{4>i\^=i- These are given by the Pontryagin indices 

Qt = J eab$i ■ {da$i X db$i) dxjdxk- (53) 

In Eq. (|53p . a, b span all spatial indices different from i: 
(a, b e {j, k}). These (integer) topological indices Qi can 
be related to the homotopy group 

U2{S^) = Z. (54) 

The integers Qi denote the wrapping numbers - the num- 
ber of times that the spatial two sphere [(the subscript in 
Eq. (j54p defined by the plane [jk] with a unique point at 
infinity] can be wrapped onto the the internal sphere 5*^ 
relevant to the three component fields (pi of unit norm. 
The likes of Eq. ([54]) are seen in numerous o(3) spin 
systems in two spatial or in 1-|-1 space-time dimensions 
yet not seen in three spatial dimensions. What is novel 
about systems with d-GLSs with d < D is that they 
allow for defects of lower-dimensional structure in high 
dimensions. 

Later on, we will return to similar forms when we will 
discuss topological terms which appear in systems with 
d-GLSs. 

In a similar fashion a multitude of other field theories 
can be constructed for systems with d-GLSs which reside 
in D > d dimensions and for which topological defects 
are characterized by topological indices which are usu- 
ally reserved for systems in d dimensions. For quantum 
spin systems. Berry phase terms may be inserted into the 
functional W in Eqs. ([5^ . and (|1D|) to enable to the same 
dimensional reduction. 

As wc discuss next, it is precisely these topological 
defects which eradicate local order for systems endowed 
with low d-GLSs and which consequently may enable 
TQO. 

C. Absence of spontaneous symmetry breaking 

Next, let us ask ourselves whether we can sponta- 
neously break d-GLSs. As it turns out [13, and reviewed 
in the Appendix, the absolute values of quantities not in- 
variant under Qd are bounded from above by the expec- 
tation values that they attain in a d-dimensional Hamil- 
tonian H (or corresponding action S) which is globally 
invariant under Qd and preserves the range of the inter- 
actions of the original systems. [In this proof, we take (as 
we must) the perturbing field favoring order to zero only 
after the thermodynamic limit is taken.] The physical 
origin of this effect lies in the large entropic contributions 
associated with low-dimensional defects [see, e.g., Fig.|9]. 
As the expectation values of local observables vanish in 
low-d systems, this bound strictly forbids SSB of non- 
Qd invariant local quantities in systems with interactions 
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of finite range and strength whenever d = (EHtzur's 
theorem) [4^ , d = 1 for both discrete and continuous C/^, 
and (as a consequence of the Mermin-Wagner-Coleman 
theorem (46j ) whenever d = 2 for continuous symmetries 
[l3|- Discrete d = 2 symmetries may be broken (e.g. the 
finite-T transition of the D ~ 2 Ising model, and the 
d = 2 Ising GLS oi D = 3 orbital compass systems). In 
the presence of a finite gap in a system with continu- 
ous d = 2 symmetries, SSB is forbidden even at T = 
[l3|. The engine behind the absence of SSB of non-GLSs 
invariant quantities in high-dimensional systems are the 
same low-dimensional topological defects which destroy 
order in c? < 2 systems (e.g. domain walls/solitons in 
systems with d = I discrete symmetries, vortices in sys- 
tems with d = 2 U{1) symmetries, hedgehogs for d = 2 
SU(2) symmetries). Transitions and crossovers can only 
be discerned by quasi-local symmetry invariant string- 
like quantities (e.g. Wilson-like loops in pure gauge the- 
ories) or, by probing global topological properties (e.g. 
percolation in more general matter coupled gauge the- 
ories [Hi, [iil). Extending the bound of to 
T = 0, we now find that if T = SSB is precluded in 
the d-dimcnsional system then it will also be precluded 
in the higher-dimensional system for quantities not in- 
variant under exact or T = emergent d-GLSs. Exact 
symmetries refer to [U, H] = 0. In emergent symme- 
tries Q unitary operators tl G tJemergent are not bona fide 
symmetries {[U,H] ^ 0) yet become exact at low ener- 
gies: when applied to any GS, the resultant state must 
also reside in the GS manifold, 



f/|5o) 



/3 



(55) 



These operators thus become symmetries when confined 
to the low-energy sector [50l |. 

Putting all of the pieces together, SSB at T = 0+ (even 
for T ^ A with A the gap size in gapped systems) of 
exact or emergent d < 1 discrete or SSB even at T = 
of 0? < 2 continuous symmetries in systems of finite 
interaction strength and range cannot occur. 



D. Absence of quasiparticle excitations 

When the bound of [l3| is applied anew to correlators 
and spectral functions, it dictates the absence of quasi- 
particle (qp) excitations in many instances [5l|] . Here we 
elaborate on this: The bound of [l3| mandates that the 
absolute values of non-symmetry invariant correlators 



\G\^\J2anAU 



(56) 



with f2j C Cj, and {ooj} c-numbers, are bounded from 
above (and from below for G > (e.g., that correspond- 
ing to (|0(k, Li;)p))) by absolute values of the same cor- 
relators \G\ in a c?-dimensional system defined by Cj. In 



particular, {anj} can be chosen to give the Fourier trans- 
formed pair-correlation functions. This leads to stringent 
bounds on viable qp weights and establishes the absence 
of qp excitations in many cases. In high-dimensional sys- 
tems, retarded correlators G _generally exhibit a resonant 
(qp) contribution (e.g., [s^. [53|V Here, 



with 



G = Gre5(k,a;) + Gnon-res(k, w) 



Gres(k,tj) 



CJ — ej. + iO+ 



(57) 



(58) 



In low-dimensional systems, the qp weight ^ and 
the poles of G are often replaced by weaker branch cut be- 
havior. If the momentum k lies in a lower d-dimensional 
region Cj and if no qp resonant terms appear in the cor- 
responding lower-dimensional spectral functions in the 
presence of non-symmetry breaking fields then the up- 
per bound [l3| on the correlator |G| (and on related 
qp weights given by lim^^^^^{LL! — ek)G(k, lj)) of non- 
symmetry invariant quantities mandates the absence of 
normal qps. If quasi-particle fractionalization occurs in 
the lower-dimensional system then its higher-dimensional 
realization follows. 



E. Charge Fractionalization 



As we noted in Section llVi whenever the center of the 
group linking the various GSs is (or contains a subgroup) 
of the Zat type, then the basic primitive charge can be 
quantized in integer multiples of 1/N. Whenever a dis- 
crete d-GLS symmetry is present for which such a frac- 
tionalization occurs, then fractionalization may occur in 
the D-dimcnsional system. If the d-GLSs have a Z^v 
group as their center then a fractional primitive charge 
of size (1/iV) is possible. The d ~ 1 quasi-particle quasi- 
hole insertion operators in the Quantum Hall problem 
[2ll . [2^ lead to fractionalization in just the same fashion 
as they do in the Peierls chain. d-GLSs are obviously not 
required for the appearance of fractional charge (see our 
discussion of 7V-ality in Section IIV[) . The most promi- 
nent of all examples of charge fractionalization [that of 
quarks {SU{N = 3))] occurs in 3+1 dimensions. What 
we wish to stress is that the presence of d-GLSs makes 
the identification of fractional charge more immediate. 



F. Topological terms 

A related consequence of systems with d-GLSs is that 
topological terms which appear in d + 1-dimensional the- 
ories also appear in higher D + 1-dimensional systems 
{D > d). These can be related to Chern number invari- 
ants in these theories whose form is similar to that of 
general topological indices which we discussed in Section 
IIXBI These topological terms appear in actions S which 
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bound quantities not invariant under the d-GLSs. For 
instance, in the isotropic D = 2 (or 2+l)-dimensional 
orbital-less general spin version of the t2g KK mo del [ See 
(b) of Section IVlII for the complete Hamiltonian.] [25|, of 
exchange constant J > 0, the corresponding continuum 
Euclidean action is of the 1+1 form 



S = — I dxdr 
where 



— (Ixdreau m ■ {dam x d^m) (59) 

OTT J 

with fijiy (z {x, t} and e^^, the rank-two Levi Civita sym- 
bol. Here, as in the non-lincar-cr model of a spin-S" chain, 
rh a normalized slowly varying staggered field, g — 2/S, 
Vs — 2JS, = 2TrS, and Str a transverse-field action 
term which does not act on the spin degrees of freedom 
along a given chain [l7l | . Q is the Pontryagin index corre- 
sponding to the mapping between the (l-l-l)-dimcnsional 
space-time (cc, r) plane and the two-sphere on which rh re- 
sides. This (l-l-l)-dimcnsional topological term appears 
in the (2-1-1 )-dimensional KK system even for arbitrary 
large positive coupling J. This, in turn, places bounds 
on the spin correlations and implies, for instance, that in 
D = 2 integer-spin t2g KK systems, a finite correlation 
length exists. 



X. FINITE TEMPERATURE TQO AND 
GAUGE-LIKE SYMMETRIES 

Our central contention is that in all systems known to 
harbor TQO (and in new examples), d-GLSs are present. 
Old examples include: Quantum Hall systems, I2 lattice 
gauge theories, the Toric-code model [4I and other sys- 
tems. In all cases of TQO, we may cast known topologi- 
cal symmetry operators as general low-dimensional d < 2 
GLSs (e.g. in the Toric code model, there are d — 1 sym- 
metry operators spanning toric cycles). The presence 
of these symmetries allows for the existence of freely- 
propagating decoupled d-dimensional topological defects 
(or instantons in {d -\- 1) dimensions of Euclidean space 
time) which eradicate local order. These defects enforce 
TQO. In what follows, we will first investigate finite-T 
TQO in general systems and then turn to gapped sys- 
tems with continuous symmetries. 



A. Sufficient conditions for finite temperature TQO 

We now state a central result: 

Theorem: A system will display TQO for all temper- 
atures T > if it satisfies the following three conditions: 
(i) It obeys the T ~ TQO conditions of Eq. {!]), (ii) all 
of its interactions are of finite range and strength, and 
(Hi) all of its GSs may be linked by discrete d < 1 or by 
continuous d < 2 GLSs U Cz Gd- 



Proof: 

Recall the T > definition of TQO as given the re- 
quirement that both Eqs. ([T]), and dH) are satisfied. For 
any Qd symmetry, we separate V = Vg^ -\- Vg;± where 
Vg_o is the component of V which transforms under the 
singlet representation of Qd (and for which [Vg^, U] = 0) 
and where Vg-± denotes the components of V which 
do not transform as scalars under Qd (and for which 
/ dU UWg.^U = 0). 

We will now first discuss the symmetry invariant com- 
ponent Vgfi and then turn to the symmetry non-invariant 
component Vg._\_. To prove the finite-T relation of Eq. 
(HI, we write the expectation values over a complete set 
of orthonormal states {|a)} which are eigenstates of the 
Hamiltonian H endowed with boundary terms favoring a 
particular GS a, 



{Vg.o)o. 



Ea(a|^e,o|a)e 



E„(a|C/Ve,oC/|a) 



{Vgfi)p- (60) 



Here, we invoked U\a) = \b), and Ea = Eb (as [U,H] = 
0). The term (f)"^ monitors the effect of the boundary 
conditions favoring the state a. If is the Hamiltonian 
endowed with boundary terms corresponding to the GS 
a then (a|i?"|a) = Ea -\- (f>%- In the above derivation, 
't'a = 'l^ua ^ '^^^c which is evident by the application 
of a simultaneous unitary transformation in going from 



\b) and \go 



Next, we turn to the symmetry non-invariant compo- 
nent Vg;±. For the non-symmetry invariant Vg,±, by the 
theorem of [l3|, {Vg^±)a = 0, i.e. for systems with low- 
dimensional GLSs, symmetry breaking is precluded. [In 
the proof of [l3|, we take the perturbing field favoring 
order is taken, as indeed it must, to zero only after the 
thermodynamic limit is taken.] Equation (|60p is valid 
whenever [[/, 1/] = for any symmetry U. However, 
{Vg^±)a = only if is a low-dimensional GLS. In sys- 
tems in which not all GS pairs can be linked by the exclu- 
sive use of low-dimensional GLSs U € Qd {U\ga) = Igp)), 
SSB may occur. 

Taken together, our results for both Vg^o and Vg-± con- 
clude our proof. It is important to note that our sym- 
metry conditions for T > TQO once T = TQO is 
established do not rely on the existence of a spectral gap. 
d-dimcnsional symmetries alone, irrespective of the exis- 
tence or absence of a spectral gap, mandate the appear- 
ance of finite-T TQO. As we will discuss in the next Sec- 
tion, the existence of a gap facilitates the proof of T = 
TQO in systems harboring a continuous symmetry. 

As stated initially, the physical engine behind TQO in 
the presence of d-GLSs are topological defects: For ex- 
ample, domain walls in d < 1 systems, Goldstone modes 
in continuous d < 2 systems. These topological defects 
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proliferate throughout the lattice to enforce TQO (and 
destroy local orders) . The topological character of the de- 
fects is particularly evident when upon duality, the (dual) 
disordering fields have the form of Eq. ([7]). Here it is seen 
how the defects link the system bulk to its boundaries. 

It is worth noting that the homotopy groups of such 
defects can, of course, be cither Abclian or non-Abclian. 
A non-commutativity of the defects implies a topological 
entanglement in the sense that entangled defects cannot 
be undone. Two defect loops can become entangled if 
their homotopy group representations do not commute. 
Such a non-commutativity of dislocations (defects in elas- 
tic solids) is responsible for the phenomena of work hard- 
ening. General defects may allow for operations which in- 
clude and extend on the particular braiding group (which 
appears for FQHE quasiparticles or Majorana fermions in 
the vortex cores oip + ip superconductors) - the principal 
non- Abclian group to have been discussed to date in the 
literature. This may be extended to general non-Abelian 
d-GLSs. 



B. Gapped systems with continuous 
low-dimensional Gauge-Like Symmetries 

We will now show that in systems with a spectral gap 
between the ground and all excited states, and in which 
continuous low-dimensional (d < 2)-GLSs are present, 
TQO can appear at low temperatures. Unlike the case for 
discrete d = 1 symmetries, the existence of a spectral gap 
for d = 1,2 symmetries ensures the existence of T = 
TQO. This then, in conjunction with the central theorem 
of Section |Xl leads to a very strong result: 

Theorem: When in a gapped system of finite interac- 
tion range and strength, the GSs ( each of which can he 
chosen by the application of an infinitesimal field) may 
he linked hy continuous d < 2 GLSs U € Gd- then the 
system displays both zero and finite-T TQO. 

Proof: 

Here, we rely on the discussion of Section IIX CI - in 
particular Corollary IV of Appendix |B] |17| . As Corol- 
lary IV is a T = extension of Elitzur's theorem, the 
proof of the above theorem for T = TQO follows word 
for word the same proof which we furnished earlier for 
T > TQO in Section IXXl As before, we will now first 
discuss the symmetry invariant component Vq^q and then 
turn to the symmetry non-invariant component Vg-_]_. For 
VQfi. (1,9/3) = t^|,9„)), {gc.\VgM = {9c.\U^VgflU\go.) = 
{9p\ygfi\9i3)- 

Next, we turn to the symmetry non-invariant compo- 
nent Vg;j_. For the non-symmetry invariant Vq,j_, by the 
theorem of p7| . and corollary IV in particular, 

{Vg,^)o. = 0, (61) 

at both zero and all finite (non-zero) temperatures. That 
is, for systems with low-dimensional GLSs, symmetry 
breaking is precluded. Eq. (j60p is valid whenever 
[U,V] = for any symmetry U. However, {Vg^±)a = 



only if [/ is a low-dimensional GLS. In systems in which 
not all GS pairs can be hnked {U\ga) = \gp)) by the ex- 
clusive use of low-dimensional GLSs U £ Gd, SSB may 
occur. Taken together, our results for both Vg^ and Vg;± 
conclude our proof. The existence of T = TQO in con- 
junction with the theorem of Section |Xl then establishes 
the existence of TQO for all T > 0. 



XI. GAUGE-LIKE SYMMETRY GROUP 
SELECTION RULES AND THEIR 
CONSEQUENCES 

It is commonly assumed that the existence of a spec- 
tral gap may render T = TQO (that satisfying Eq. ^) 
stable at small positive temperatures. In Section IX Bl 
we rigorously proved (for the first time) that this is in- 
deed so for continuous low d-GLSs for the appearance of 
TQO. We now turn to the more general analysis of T = 
conditions of Eq. ([1]). We will illustrate that symmetry 
alone (without ever invoking the presence of a spectral 
gap) can guarantee that the conditions of Eq. ([T]) are 
satisfied. Unlike finite-T TQO which follows from low-d 
{d < 1 discrete or c? < 2 continuous) GLSs, the T = 
rules which we will discuss apply for general d > 1 GLSs. 
Our approach highlights how the states themselves en- 
code TQO. The selection rules which we will employ be- 
low are independent of the specific Hamiltonian which 
leads to those GSs. 



A. Selection rules for general Gauge-Like 
Symmetries 

Let us first turn to the off-diagonal portion of Eq. ([T|). 
This relation is often seen to be immediately satisfied in 
the basis of GSs which simultaneously diagonalize not 
only the Hamiltonian but also the d-dimensional sym- 
metry operators. Here, we may invoke symmetry based 
selection rules (and, when applicable, the Wigner-Eckart 
theorem in particular) to show that the matrix element 
of any local operator V between two GSs which are eigen- 
states of the d-dimensional symmetry generators vanishes 
in the thermodynamic limit. This, along with a set of 
equations which enable us to satisfy the diagonal por- 
tion of Eq. Ill) will assert the existence of T = TQO. 
For concreteness, in what follows, the notation will refer 
to d = 2 SU{2) symmetries. Our results hold in their 
exact form for symmetry groups in which the Wigner- 
Eckart theorem holds (e.g. integer spin variants of these 
SU{2) symmetries). When applicable, we will comment 
on the occurrence of similar results in systems in which 
the Wigner-Eckart theorem does not apply (e.g. half- 
integer spin variants of these SU{2) symmetries). The 
d = 2 rotational symmetries, which we will use for illus- 
tration, act on two-dimensional planes {Pi}. A realiza- 
tion of these symmetries is afforded by the t2g KK model 
which wc will discussed at length later on [See (b) of Sec- 
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tion IVIIj |55| . The extension symmetries (of dimension 
d 7^ 2) of groups other than SU (2) is straightforward. 

Let us assume that a given quasi-local operator V has 
its support inside a sphere of radius R. This means that 
the operator V spans no more than (2i?) planes and in 
each plane its spans no more than (4i?^) sites. We can 
expand 1^ as a sum of planar operators, each of which 
has its support on a plane Pj, i.e. 



(64) 



for any local operator Vp. which has its support exclu- 
sively on the plane Pj then any orthonormal GSs of the 
form of Eq. dM]) wiU satisfy Eq. Next, let us ex- 

amine in detail this matrix clement, {(l)p\Vp.\(j)p.) . We 
will employ selection rules on the matrix elements of the 
irreducible tensors {T'^'}. Whenever the Wigner-Eckart 
theorem is applicable. 



{a J m\T^'^\aj'm) = 



1 



V2j + 1 ^ ^ 

x{jm;kq\j'm), 



(65) 



we will have that for any local operator Vp, (henceforth 
denoted by V for brevity), the matrix element between 
two orthogonal eigcnstates of the d > 1-dimensional sym- 
metry operator 



(ajiTOi|y|a'j2'7i2) ^ 



(66) 



in the thermodynamic limit (L oo) if ji ^ j2 or 
iTT-i 7^ The proof of this assertion is given by the 

asymptotic form of the Clebsch-Gordan coefficients for 
SU{2) for systems of finite spin S. [For systems of fi- 
nite spin, products of more than S on-site spin operators 
which appear in V (if they indeed appear) can be reduced 
to lower order terms.] 

A similar extension of the more familiar variant of the 
Wigner-Eckart theorem holds for general groups other 
than SU{2) with similar conclusions. Even when the 
Wigner-Eckart theorem is not applicable - as in e.g. half- 
integer spin representations of SU{2) - similar conclu- 
sions may be seen to hold. The proof of this assertion 
for half-integer spins follows by an explicit construction 



of orthogonal eigenstates of the d-dimensional symmetry 
operator in such cases when a = a' ,j = j' . For example, 
for the spin S — 1/2 d = 2 SU{2) symmetry of pertinence 
to the KK model, we have in the sector of maximal total 
spin of the plane j = N/2 with N the number of sites in 
the plane, the states 



\J 



V= J2 Bp,P,-.P2nVp,Vp,---Vp,,, (62) \j 

with bounded coefRcients Bp-^p^...p^j^ . Let us furthermore 
write the GS wavefunctions in the direct product basis 
of individual planes, 

\gi) = J2 ^Pi-PJ<^'A) ® \^%) ® • • ■ \^'p\)- (63) 

Here, L is the number of planes. Let us denote the num- 
ber of states {li/ip^. )} in the sum by Gj. Each state \(p^p.) 
has its support on an entire plane Pj. Now, it is easily 
verified that if 



N 
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N 
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U---iT---T> 



all other states with u down spins 



(67) 



These orthogonal states are related by the exclusive use 
of the d-dimensional symmetry operators. The matrix 
element of any local operator between such orthogo- 
nal eigenstates of the d-dimensional symmetry operators 
tends to zero in the thermodynamic limit {N oo). For 
instance, the matrix element of the spin lowering opera- 
tor at site i. 



N _ N N 



1 



N N- 



0. (68) 



As in the discussion prior to Eq. ^67\ . N denotes the 
number of sites in each d = 2-dimensional plane. Simi- 
larly, the matrix elements of all quasi-local operators are 
bounded from above by (normalization) factors which 
scale to zero in the thermodynamic limit. The factors 
for the case of j = ^ provided in Eq. (p7)) . Here, 
for large u, the matrix element of a local operator be- 
tween two orthogonal states becomes exponentially small 
(©(A^""/^)) as can be seen by invoking the Stirling ap- 
proximation. We next construct states for which the con- 
ditions of Eq. ID) on all quasi-local operators V are the 
strongest and then regress to the more general case. 

For states in which |m — m'| is larger than r (the range 
of the operator V for an 5 = 1/2 system or in S* > 
1/2 systems the number of independent spin fields which 
appear in the various product terms forming a general 
quasi-local operator V), we will have that 



{ajm\V\ajm') = 0, \m — m'\ > r. 



(69) 



Let us next expand the states {|0p )} in the complete 
orthonormal basis of the d-dimensional symmetry gener- 
ators. In particular, let us focus on those states which 
may be written as 



= (r+l)n 



ajm 



ajm) , 



(70) 



with n an integer. In Eq. ([701) . "^"^ omitted the plane 
index Pa. In order to avoid any complications, we chose 



20 



the states in Eq. ([70]) such that there is no ofF-diagonal 
matrix element {ajm\V\ajm') = for any two states 
\ajm) and \ajm') which appear in the sum of Eq. (|70p . 
This vanishing matrix element follows from Eq. (j69p . 

As a consequence of Eq. (|66| , the matrix element of a 
local operator between two states is 



m={r+l)t 



{af-'"'')*{af,^"^){ajm\V\ajm). (71) 



Next, when allowed (e.g. for integer spin renditions of 
SU{2) systems), we expand V in the irreducible repre- 
sentation of the c?-dimensional symmetry. 



akq 



(72) 



In the case of the d ~ 2 planar SU{2) symmetries of 
the KK model, {T"*-'"?} behave as spherical tensors under 
rotations of all spins in the plane. Here, k denotes the 
total spin and q the magnetic quantum number of the 
spherical tensor T. As 1^ is a local operator, the total 
angular momentum that it can carry is bounded from 
above by a finite number (of the order of its range), k,q = 
0(1). Here, and in what follows wc will make extensive 
use of Eq. ((69)l . We now return to the definition of T = 
TQO given by Eq. ([T|). Obeying Eq. ([51]) . from which 
Eq. ID) is satisfied, is tantamount to having 

i"^'?<5„,+c„, ==^(ajm|T"'^-'^|ajm)(ar^")*aj;^" (73) 



with arbitrary constants {t"'"^} for all 1,1' and for finite 
k, q which appear as components of the local operator V. 
By the Wigner-Eckart theorem, Eq. ([75)1 reduces to 

e'"'Su,+cu' = ^^^^=JY^}2^jm;kq\jm) 



x{ar"ri4n, (74) 



for a finite set of (fc, q) values with {jm; kq\jm) denot- 
ing the Clebsch-Gordan coefficients. The double bar in 
"(aj||T"'^||aj)" is used to emphasize that this coefficient 
is independent of m and m'. Selection rules mandate 
that q ~ 0. For m ~ {r + l)n in the sum of Eq. ([7^ . the 
corrections {cw} may be set to zero. 

We are now faced with the task of satisfying Eq. (fM]) 
for a finite set of (fc, q) values and for arbitrarily chosen 
constants {f"'^'}. In general, Eq. ([7^ - a compact short- 
hand for a system of quadratic forms in {a"j^,a"j™} 
with the subscripts R and / denoting the real and imag- 
inary components of a""'"' and for arbitrary engineered 
constants {t"'"^}- admits multiple solutions. Let us now 
show this. If we have y ~ 0{r) possible {k,q) values 
in each representation a for the decomposition of the fi- 
nite range V in Eq. ([75]) then, recalling the definition 
after Eq. on the left hand side of Eq. ([71]), we 

have 0{yG^) independent entries. [We remind the reader 



that G denotes the number of independent states {|(?i>p^.)} 
on each plane Pj which appear in Eq. ([63]) ]. On the 
right-hand side we have far more entries as the symme- 
try indices m also appear: On its own, for high-symmetry 
states (i.e. j = 0{N)) the magnetic number m can span 
0{N) values. Eq. ([74]) generally admits many solutions 
as — > c». 

Let us next avoid the use of Eq. (|69p and examine 
the most general situation. Repeating the steps that we 
undertook above, we will find that 



(aj||r°fe||«j) 



X ^(jm;fc9|jm')(a"^")*(a" 



ajm 
I' 



(75) 



In Eq. ([75)) . locality is manifest in the allowed values of 
the index k: i.e. k < r ~ 0{1). As a consequence of the 
behavior of the Clebsch-Gordan coefficients, (jm; kq\jm') 
tends to zero as |m — m'| > r. As we showed earlier, as 
a consequence of the asymptotic scaling of these coeffi- 
cients, we find that for local operators, the matrix ele- 
ments between different eigenstates of the d-GLSs tend, 
exponentially, to zero in the thermodynamic limit. If we 
do not restrict ourselves to values of m, m' = n{r -\- 1) 



as we have in Eq. (|74p . we will in general need to insert 
these small matrix elements in the sum of Eq. ([75]) . For 
SU{2), the matrix elements of the quasi-local operators 
(for large j = 0{L<^) and finite k = 0{r)) lead to the 
appearance of the general Clebsch-Gordan coefficients 

{jm;kq\jm) = 



(2j + l)[fc!]2(j -hm')!(.? - m')! 



X J2 (-1)^+'"'+'' 

p=0 



{2j + k + l)!(j - m)!(j 4- m)!(fc - g)!(/c + g)! 

(fc -|- J -f m — p) ! (j — m + p) ! 
p\(k — p)\{j + m' — py.(j> + j — k — m')\ 



in Eq. ([75]) . For large j — 0{L'^), these matrix elements 
lead to the exponentially small terms discussed earlier. 

For groups other than SU (2) we have a similar gen- 
eral relation from the Wigner-Eckart theorem [l3| . This 
concludes our construct. Note that in the T = case, 
we do not require that all GSs be linked to each other 
by the exclusive use of the d-dimensional symmetry (if 
that were the case only the magnetic number m changes 
yet j is held fixed in going from state to state). The 
d-dimensional symmetry simply provides a very conve- 
nient basis for which the construction of TQO is simple. 
Our Wigncr-Eckart-typc construction and related vari- 
ants enable us to construct general states with TQO. We 
are currently working on extensions of this method to 
other groups. For the FQHE, the associated group is 
that of magnetic translations (see Section IVip . 

It is noteworthy that, in this construct, we do not de- 
mand that the d-dimcnsional symmetry generators com- 
mute. In the above construction, wc assumed that we 
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have d-dimcnsional symmetry generators, each generator 
acting within a different plane, and showed how to con- 
struct GSs with TQO. The presence/absence of any ad- 
ditional symmetries which do not commute with these is 
inconsequential. Similarly, the dimension d of the GLSs 
need not be small. Any d > 1 GLSs will allow us to fol- 
low this construct. All that matters is that all GSs can 
be written in the form of Eqs. and (|7n)). 



B. A special set of states 

The discrete cyclic symmetry groups Z„ have the sim- 
plest selection rules amongst all GLSs. These discrete 
symmetry groups, either appearing on their own right or 
as homotopy groups for continuous d-GLSs furnish dis- 
crete topological numbers. 

Nearly all of the prominent examples of TQO to date 
exhibit a discrete Abelian d ~ 1 GLS. There the analysis 
simplifies considerably. In Section IXIH we will show in 
detail how the T ^ conditions are satisfied for gen- 
eralized Kitaev type models. The Kitaev model dis- 
plays a d = 1 Z2 symmetry. In what follows, we note 
how for these and other general systems a special set 
of states may be constructed for which TQO appears. 
These states encompass the GSs of the Kitaev model as 
well as those of the Rokhsar-Kivelson (RK) model [56j . 
In the Kitaev and RK models, the role of the representa- 
tion indices of the d-GLS groups is taken on by a discrete 
parity (Z2) index which captures the d = 1 Z2 symme- 
tries which are present in these systems. Our construct 
below is more general albeit still being quite special. 

The states of general Abelian theories can also be 
viewed as particular substates of a richer non-Abelian 
theory, since non-Abelian symmetries (such as the ones 
discussed in the previous Section) contain certain Abelian 
symmetries as special subgroups. For instance, Z,i C 
[/(I) C SU{2) • • • . A very special member of the general 
states of Eq. ([75)1 which exhibit TQO is that of the states 

|g„)-AA J2 /W|co>: (76) 

|c„): |c„)ea 

in which we sum over all states in the local spin repre- 
sentation which lie in the topological sector a, in which 
the number of states in each topological sector is the 
same. The label of the topological sector is generally none 
other than the representation index of a discrete Z„ d- 
dimensional symmetry. The function /(c) in Eq. (|76p is 
such that there is a one-to-one correspondence between 
individual states in the topological sectors |cq) ^ jc/j) 
with /(cq) = f{cp). The coefficients /(c) need to be 
chosen in order to ensure Eq. ([1]). We will return to this 
form in Section IXIII 

As noted above, an example is afforded by the GSs 
of both the Kitaev (Eqs. ([TU)) . and pT|) ). and Rokhsar- 
Kivelson models which we will discuss in more detail in 
later Sections. In the GSs of these systems, /(c) = 1 



in Eq. (|76|) . and the GSs of the completely symmetric 
form of Eq. (|92p . Here, the sectors a are determined by 
the winding number in a d = 1 surface. This winding 
number - a consequence of the d = 1 symmetries that 
these systems display in their low-energy sector - takes on 
the role of the general group representation index which 
we employed in our analysis of the former Section. 

As we will later return to (Section IXIip . in the Kitaev 
model of Eqs. (fTUl) . and pT|) . there is a d = 1 symmetry 
selection rule 

{9c\X{<\9p)=^. a 7^ A (77) 

aef2 

with fi = x,y, OT z if the volume O is finite (the operator 
V is local). This mandates that the off-diagonal condi- 
tion of Eq. ^ is satisfied. We now prove the diagonal 
portion of Eq. ^ for the states of Eq. ([75)) which consti- 
tute the solution to the GS problem of the Kitaev model. 
Let us denote by T^p the operator that generates the one- 
to-one correspondence between the local spin basis states 
in the topological sectors a and (3, i.e. Tq^|cq) = |c^). 
Let the states {|cq_/3)} be eigenstates of ^i^j^o^i with 
A all lattice sites. In the Kitaev model, there is, up to 
a gauge transformation (an application of As), one such 
operator for a given horizontal/ vertical winding number 
sector such as that of a or (3, 

r=n^s-' (78) 

where we may always deform the horizontal/ vertical con- 
tour P by multiplying it by Ag (as As = 1 in the GSs) 
such that it does not traverse the region of support CI of 
the local operator V. All that matters is that P winds 
once in the horizontal and/or vertical directions. Here, 
trivially, 

[T, V] = 0. (79) 

Consequently, the diagonal matrix element of V 

(ffalV^lffa) = {gc.\TWT\g^) = {gplVlgp). (80) 

When combined with the vanishing off-diagonal matrix 
element of V in Eq. ([77]) , this proves the diagonal portion 
ofEq. jSl). 

When expressed in the original SU (2) spin language, 
the maximally symmetric state of the Kitaev model so- 
lution correspond, in the d = I total spin eigenbasis to 
the 4 states highest spin {Stot = L/2) states 

l^±±)- E Cs.si\Stot^L/2,S^,,S'^). (81) 

Sl'l^=L/2~ even/odd 

The superscript v and h refer to vertical and horizontal 
directions. Here, the symmetry operators along the hori- 
zontal and vertical directions commute with one another 
and we can employ the basis states {\S^^^, S"^] S^g^, S'^)} 
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as we have in Eq. (jSip . In the thermodynamic Hmit, 
where the total spin states do not admit off-diagonal ex- 
pectation values of local operators, the constant C in Eq. 
([5T|) is akin to / in the general Eq. ([75)1 . As we noted 
earlier, the Z2 operators (of Eq. (fTTjl ) which appear in 
the Kitaev problem are a particular subset of the larger 
SUi2) group. 

C. Implications for the energy spectrum - exact 
degeneracy of states in the thermodynamic limit 

A corollary of our construct is that for any Hamiltonian 
which is of finite range, i.e. for any Hamiltonian which 

can be expressed as 

H = J2V., (82) 

i 

with quasi-local operators {Vi} which have their support 
in a domain including the site i, degenerate states can 
be immediately constructed. Applied to Eq. ([82l) . our 
conditions imply that all states \(j)i) (and general combi- 
nations thereof such as those appearing in Eq. ([63]) ) are 
degenerate in energy up to the exponential corrections (if 
all cw in Eqs. ((74)) . and ([TS]) are bounded from above by 
exponentially small quantities). This is so as 

\{gi\H\g,)-{g,\H\g,)\ = \J2ic}i-c],,)\ 

i 

< c^L° exp(-a/L) 0, (83) 

on a hypercubic lattice of side i as L — > 00. In Eq. 

for all states {l^a}} [a = 1,1') and for all quasi- 
local operators Vi which have support at a site i, we 
have {gi\Vi\gv) = ViSu, +c\i,, \c\i\ < |ci| exp(-aL'') with 
a constant a > [for a d > 1-GLS], and = maxi{|ci|}. 

D. The failure of this construct for non-topological 

(local) symmetries 

It is important to emphasize that a crucial ingredient 
in our construct, Eq. ([66]) . relies on the eigenbasis of 
d > 1 symmetry generators. In the presence of local 
{d = 0) symmetries alone [i.e. if no other symmetries are 
present], this relation fails. For d = symmetries, the 
eigenstates of the symmetry generator contain a finite 
number of sites and the expectation values of local oper- 
ators between two orthogonal symmetry eigenstates can 
be finite. An example is afforded, e.g., by the classical 
Z2 lattice gauge theory 

H = -KY,<^l^^-k<yH< (84) 

□ 

[which is Hz^ of Eq. with hx — 0] for which the 

local symmetry operators 

G,=\[al, (85) 



with r all nearest neighbors, the site i link orthogonal 
GSs. In this particular case the local operators {Gi} 
arc none other than {As} of Eq. (fTTj) which defines the 
Kitaev model. The local symmetry operators V = Gi 
clearly violate Eq. ([2]) and the Z2 gauge theory of Eq. 
(jMI does not exhibit T = TQO for all of its GSs. A 
partial set of GSs of Eq. ^ does exhibit TQO. This 
partial set for the D = 2 theory is given by the GSs of 
the Kitaev Toric code model which we discuss next. The 
finitc-T TQO for Eq. ([84|) is immediate and follows from 
Elitzur's theorem. 



XII. ENGINEERING TQO: GENERALIZED 
KITAEV-TYPE MODELS 

In this Section, we will show how many models can 
be systematically engineered to have T = TQO. This 
will extend and complement some of the general results 
introduced in Ref. [2| as well as our results of Sections 
Ixra andim] 

We start with a general Hamiltonian H which has GSs 
which may all be linked to one another by the exclusive 
use of discrete d < 1 symmetries or continuous d < 2 
symmetries. Note that H may in general have d = 
symmetries. Henceforth, we discuss systems which have 
such quasi-local symmetries and we label these symme- 
tries by the unitary operators {Gi}. We (i) discuss cases 
in which all of these symmetries commute with one an- 
other. As [H, Gi] = 0, we can simultaneously diagonalize 
Gi and H. (ii) In what follows, we consider models in 
which the GS sector of H contains, amongst others, states 
with a uniform unit eigenvalue under all {G,}, (iii) If a 
local quantity is invariant under all of the d = symme- 
tries {Gi} then it is also invariant under all remaining 
d > symmetries which are necessary to link orthogonal 
GSs. Finally, (iv) we consider systems in which all ba- 
sic local observablcs transform as a singlet under {Gi}, 
G\v,G, = \v■,^V,. If [V, G,;] ^ then Ay,, 7^ 1. 

Let us start by considering 

H = H -h^G^, (86) 

i 

with h > 0. We claim that under the conditions stated, 
the GSs of H display TQO. The proof of this assertion is 
simple: From conditions (ii) and (iv), it follows that for 
all GSs of H, are the subset of the GSs of H which have 
Gi\ga) = \ga)- On the one hand, the expectation value 

{g^\GlVG,\~g^) = Xv49a\V\ga). (87) 
On the other hand, 

Gilfja) = \ga), 

{g^\GlVG,\~g^) = {~go.\V\~g^). (88) 

If there is at least one i for which [V, Gi] 7^ then it 
follows that 

{~9o.\V\~g^) = 0. (89) 
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Putting all of the pieces together, this guarantees that the 
diagonal portion of Eq. ^ is satisfied: (a) If [V, Gi] = 
for all i then {ga\V\ga) = {gf3\V\gp) as V is invari- 
ant under symmetries linking \ga) and \gp) (assumption 
(iii)). (b) If [V^,Gi] 7^ then, as we showed, for all GSs, 
{~go.\V\gc.) =0. 

It is clear that condition (ii) is a necessary condition 
for TQO: For any local symmetry Gi, there exist many 
states GSs for which Gi\ga) = \ga)- 

As we showed earlier in the proof of our finite-T result, 
condition (ii) suffices to establish that the off-diagonal 
portion of Eq. ^ is obeyed. 

Let us now discuss how the Kitaev model is a spe- 
cial instance of such Hamiltonians of Eq. ([86|) . Here, 
H = —J^n'^ij'^jk'^ki'^ii is Hamiltonian of the Z2 
gauge theory which is the sum over all plaquettcs of the 
operator Bp (Eq. (fTTj) '). Corresponding to H are the 
local symmetry operators 
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These are the As=i symmetries of Eq. (|lip . With h — 
1, the corresponding H of Eq. ((86)) is none other than 
the Kitaev model. The d = symmetries of H are just 
{As, Bp} (see Eq. PT|) ). We now verify that conditions 
(i) - (iv) are satisfied for the Kitaev model: (i) All of 
these symmetries commute with one another. 



[As, As,] = [Bp, Bp,] = [As,Bp]^Q. 



(91) 



Condition (ii) is satisfied for any of the four maximally 
symmetric states 0: 



i^,)=2i-^=/^5:ic), 



(92) 



where the sum is performed over all states |c) in the cr^ 
basis which have Bp = 1 on every plaquette and for which 



the two Z2 operators {^1,^2} of Eqs. (|12p assume one 
of the four specific values q = (±1,±1). Note that [tpq) 
is the sum of all GSs of the classical Z2 model (that with 
the transverse field /i = in Eqs. and ([5^ ) which 

belong to the topological sector q. Eq. ((92)) is a special 
realization of Eq. (|75)) . 

(iii) Here, if any local operator V is invariant under 
all of these d = symmetries then V must be invariant 
under the d = 1 symmetries of Eq. (fT2)) . To see this, 
we note that for any such operator V = f{{As}, {Bp}). 
As a functional of {As} and {Bp}, V is automatically 
invariant under the d = 1 symmetries of Eq. (fT^ . Lastly, 
condition (iv) is straightforwardly satisfied: 



As<JijAs 

^s<^klAs 



'kl 



is^k,l). 



and 



AsCrfjAs 



(93) 



(94) 



[Similarly, BpufjBp = 
plaquette p, Bpa'fjBp 



-afj if the bond (ij) lies in the 
: afj if the bond (ij) does not 



lie in p, and Bpaf^Bp = crf^-.] As a consequence of 
these relations, any multinomial Vi in these fields which 
contains, amongst others, fields at the site i satisfies 
G\ViGi = Xv-iVi with Gi = As- Putting all of the pieces 
together, the reader can convince him/herself that Eq. 
([5^ constitutes the solution to the GS problem of the 
Kitaev model. It is noteworthy that here TQO (and the 
GSs of Eqs. ()92p appears for all h > 0. The generaliza- 
tion to other Hamiltonians H is straightforward. Any GS 
of the Kitaev model can be expressed as a projection of 
a reference state to the GS sector (that in which As — 1 
for all sites s and Bp = 1 for all plaquettcs p) : 



1 + B.„ 



(95) 



Here, A/" is a normalization constant and [(f)) is a reference 
state. All GSs of the Kitaev model can be linked to one 
another by the exclusive use of d = 1 symmetries. Due to 
our theorem [see Section |X] TQO follows for all T > 0. 
This is so as our analysis above shows that TQO appears 
at T = 0. 

It should be noted that other prominent examples of 
TQO (e.g. the Quantum Dimer Model ^) have TQO 
because of the very same selection rules for the d-GLSs 
which connect their GSs. In the Appendix [E) we discuss 
certain technical aspects of the Quantum Dimer model. 



XIII. ENERGY SPECTRA AND THERMAL 
FRAGILITY OF SOME ANYONIC SYSTEMS 

The entire focus of our work is on the non-local sym- 
metries {d > 0) which relate and define different topo- 
logically ordered states; we showed how the excitation 
spectrum associated with the restoration of these sym- 
metries through the appearance of low-dimensional topo- 
logical defects can eradicate local orders and enforce the 
appearance of TQO In most cases, this dimensional 
reduction insofar as SSB is concerned does not lend it- 
self to an exact dimensional reduction in the form of the 
associated free energies. In this Section, we dwell on 
several symmetries harboring low d-GLSs for which an 
exact dimensional reduction (of the partition function) 
also occurs. For example, Kitaev's and Wen's models 
of Eqs. HH]), mi), and dH]) display d = 1 Z2 symme- 
tries and fortuitously these D = 2 systems can indeed 
also be mapped onto nearest-neighbor Ising chains (of 
D = d = 1). This accidental exact reduction of the par- 
tition function for these special models points at much 
broader relations - those of the insufhciency of the spec- 
tra in determining whether or not TQO appears and in 
illustrating the thermal fragility of topological quantities 
in systems such as these. Apart from being mathemat- 
ically interesting, these systems and extensions therein 
might pave the way for TQO computing. Kitaev's Toric 
code model [3 was the first surface code [HHl suggested 
for anyonic quantum computing. 
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As we stated in dualities illustrate the fact that 
the quantum states themselves and not the energy spec- 
trum in a particular (operator language) representation 
encode TQO. The energy spectrum of theories with and 
without TQO may be the same. If two systems have the 
same energy spectrum, then they are related by a unitary 
transformation that links the two sets of eigenstates. If 
one system has TQO while the other does not then this 
transformation is generally non-local. 

That the spectrum, on its own, is insufficient to de- 
termine TQO [1] is estabhshed by counter-examples, e.g. 
that of the D = 3 Z2 lattice gauge theory which is dual to 
a. D = 3 Ising model [G^I . Albeit sharing the same energy 
spectrum, in the gauge theory we may find GSs with, in 
the thermodynamic limit, rank-n = 8 (see our defini- 
tion of rank-71 TQO in Section |TTT| finite-T TQO while 
the Ising model harbors a local order. Similarly, the Z2 
gauge theory on a square lattice [given by Eq. (|13p with 
/i = 0] is equivalent to the classical D = 1 Ising model. 
This last, well known, equivalence can be proven by, for 
example, replicating the analysis which we will perform 
in this Section for both Kitaev's Toric code model and 
Wen's plaquette model. The D = 1 Ising model has no 
TQO while the D = 2 Z2 gauge theory (which is dual 
to it) has rank-71 = 4 TQO. [Four GSs of the D = 2 
Z2 gauge theory which satisfy Eq. ([T]) are given by Eq. 

m-] 

Another example is given by the Kitacv Toric code 
model of Eqs. pU]) and (|lip in a system with open 
boundary conditions. [In our convention, on the square 
lattice with open boundary conditions, the star opera- 
tor appears only if all of the 4 bonds which make it up 
appear in the lattice. We do not allow for star opera- 
tors near the boundary of the system where only three 
or two operators appear in the product of Eq. (fTT|) defin- 
ing the operator Ag .] As all of the Ising operators Ag , Bp 
commute with one another, the spectrum can be easily 
determined. Here, the density of states 

^ M' 

9{E) > SE,M~2m—r7JJ TT, (96) 

■^-^ 771 Ai — 771 

m=0 ^ ' 

which is identical to that of a D = 1 Ising model of M 
sites with open boundary conditions or alternatively of 
M decoupled Ising spins in an external magnetic field of 
unit strength. Here, on the square lattice of Ng sites with 
open boundary conditions, there are a total of 

M ^ 2Ng ^ G^/Ws + 5 (97) 

independent fields {As,Bp}. The specific form of Eq. 
(|97p for a square lattice of open boundary conditions with 
Ns ~ (i + 1)^ sites follows as there are plaquettes 
{Bp} and [L — 1)^ stars {^s}- The sum of the number 
of plaquettes and stars is given by Eq. (|97|) . The origin 
of Eq. (j96p and similar relations is the mapping that we 
can perform between each bond (in a, D = 1 Ising model) 
or each site (for a system of decoupled Ising spins in a 



magnetic field) and the M decoupled commuting Ising 
variables {Ag} and {Bp}. An identical correspondence 
between the two systems [Kitaev's model and the Ising 
chain] in the presence of periodic boundary conditions 
will be detailed in Section IXIII Al 

Any system (whether TQO ordered or not) with (for 
finite L) a finite number of states (P), has a spectrum 
which is identical to that of a spin system with T = 
TQO. The proof of this assertion follows from our 
discussion in Section IXI Al To pictorially summarize 
this state of the affairs, we write 

TQO No TQO 

D=3 Ising gauge ^ £'=3 Ising model 

D=2 Ising gauge ^ D~l Ising model 

Kitaev's model — Wen's model ^ _D=1 Ising model 

A spectral equivalence exists between other pairs 
of TQO and non-TQO systems. Thus, we cannot, 
as coined by Kac, "hear the shape of a drum" . The 
information regarding the existence/absence of TQO is 
in the eigenvectors. We will further show that at any 
finite temperature, some of the toric cycle operators (or 
in general other d-GLSs operators) may have vanishing 
expectation values unless they are held fixed by the 
application of an external field. This occurs notwith- 
standing the existence of spectral gaps. Commonly, the 
existence of a finite gap between the ground and excited 
states is assumed to protect properties associated with 
T = TQO up to a small finite temperature. Our results 
show that this assumption is, in general, incorrect. In 
the Letter of [1], we emphasized these properties. In 
the following Sections, we will provide details for these 
claims. 



A. Dimensional Reduction and Thermal fragility in 
Kitaev's Toric code model 

In this Section, we illustrate that, as a consequence 
of a dimensional reduction to a Z? = 1 Ising model, the 
Kitaev Toric code model exhibits a loss of its toric oper- 
ators expectation values at any finite T (i.e. T > 0, no 
matter how small the temperature is). This loss occurs 
notwithstanding the existence of a spectral gap in this 
system. This does not imply that the Kitaev Toric code 
model does not have finite-T TQO. Indeed, as we showed 
by our theorem, this model does exhibit TQO at all T. 
If the fields on a chosen contour (a boundary) are not 
allowed to fluctuate then the toric operators arc robust 
due to the absence of correlations between the boundary 
and all interior sites. Fixing these boundary fields has no 
influence on the expectation value of any quasi- local oper- 
ator V: Eq. (HI) is satisfied. If, however, we were to allow 
the fields along a toric cycle to fluctuate then although 
at T = 0, the GSs have finite values of the Toric code 
operators, at finite T > 0, there is no SSB which leads to 
finite expectation values of these operators. The Kitaev 
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Toric code model has a thermodynamic phase transition 
at T = 0. We similarly discuss Wen's plaquette model. 
We further illustrate in this exactly solvable system the 
important property of adiabaticity: exact GLSs can be 
lifted while not changing the phase of the system and 
discuss within this particular system in the framework of 
our general Ising chain mapping, the dependence of the 
GS degeneracy on topology. 



1. Mapping of the Kitaev model to Ising chains and its 
consequences 

To prove the assertions above, we compute the exact 
(generating) partition function of the Kitaev Toric code 
model in the presence of applied infinitesimal fields. 



tr 



{<Tf^.}exp( -/3 



H 



K 



hx,lXi — hj:,2X2 — hz iZi — ft-z 9^2 



The operators appearing in the argument of the trace 
arc those defined in Eqs. (dU]), (HI]), and The op- 

erators {^1,2} and {Zi^2} [see Eq. ([T^ ] obey (abelian) 
anyonic statistics. For a square lattice having Ng sites 
with periodic boundary conditions. 



(99) 



The partition function of Eq. 



Z = tr cxp - /3 \Hk 

- hx,iXi - hx,2X2 - h,,iZi - h,,2Z2 ) . (100) 



Equation (jlOOp with the constraint of Eq. ([55)1 is none 
other than the partition function of two decoupled circu- 
lar chains of non- interacting S* = 1/2 spins of length 
each. The partition function is given by 



in terms of which the Hamiltonian of the decoupled Ising 
chains is 



H 



(103) 



In what follows, we drop the polarization index (z). This 
polarization index can be forgotten henceforth until the 
end of this Section when we write the second half of the 
duality relations of Eq. (|102p 
The Ising bond fields {ag, 
satisfy 



bp} on each circular chain 



n as = n 



1 



(104) 



along with the same algebra as the star and plaquette 
terms {Ag} and {Bp} of the Kitaev Toric code model 
(see Eq. dn])), 

[as, Us'] = [bp, bp'] = [as, bp] = 0, 

al^bl^ 1. (105) 

The number of independent fields in each case (the Ki- 
taev Toric code model vis a vis two decoupled circular 
Ising chains of length Ns ) is exactly the same - there are 
2Ns independent spin variables aij versus (for Ns Ising 
variables {(Js} and Ns Ising variables {(Jp}) 2iVs indepen- 
dent spin variables {as, bp}. In both cases, these variables 
satisfy the same set of constraints [Eq. (|104p for the Ising 
chains vs Eq. for the Kitaev model variables] and 

the same number of Ising variables which satisfy the same 
algebra and span the same Hilbert space of size 2^^'". 
Moreover, they have exactly the same Hamiltonian (that 
of Eq. pT7)l ). In the presence of additional fields (see last 
terms in Eq. (jlOO|) ). we may define the two component 
vector 



\hi\ 



(106) 



and set 



[(2cosh/3)^= + (2sinh/?)^=]' 
X [cosh/?/ii] [cosh/3/i2] 



(101) 



with hi = y h'^ i + ^. It is important to emphasize that 

Eq. (jlOip constitutes the partition function for the Toric 
code model for arbitrary size systems Ns on the torus. 
That this is the partition follows from, for example, the 
identification of a similar algebra in the D = 1 Ising 
model. Let us consider two decoupled Ising chains, each 
of length Ns and satisfying periodic boundary conditions 
and examine the algebra of their bonds. To this end, let 
us define the bond variables 



(102) 



Qi = Xin^i + ZiUzi, 
Q2 = X2nx2 + Z2nx2 



(107) 



to be the counterparts of two single Ising spins crj^ and 
Up^ which are located at site numbers c and d of the two 
respective Ising chains (that of the s and that of the p 
varieties) . c and d can be chosen to be any integers such 
that < c,d < {Ns — 1). Here, the algebra and set 
of constraints satisfied by the two chain Ising variables 
{as, bp, al^, cTp^} which span a 2^^= -dimensional space is 
exactly the same as that satisfied by the Kitaev vari- 
ables {As,Bp,Qi,Q2} which span the same space. The 
terms {—hiQi) of Eq. (|100p maps onto magnetic field 
terms acting at a single site (c) of each of the two Ising 
chains. The partition function of a circular Ising chain 
with a magnetic field of strength h acting at only one site 
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c and having no on-site magnetic field on all other sites 

is (^[(2cosh/3)^» + (2sinh/3)^-]cosh/3/iy In our Kitaev 

model mapping, we have two such Ising chains. This 
leads to the partition function of Eq. (jlOip . An equiva- 
lent derivation of this partition function employs a high 
temperature expansion (see e.g. [13] for an application 
of these expansions in plaquette actions similar to those 
of the Kitaev model). As seen by inverse Laplace trans- 
forms of the two partition functions (that of Kitaev's 
model and that of the two decoupled Ising rings), our 
mapping implies that the degeneracy of each level in Ki- 
taev's model is equal to the degeneracy of two decoupled 
circular Ising chains at the same energy level. In par- 
ticular, the four GSs of the two decoupled Ising chains 
(two GSs (all spins up or all spins down) for each of the 
two decoupled chains) imply the existence of four GSs of 
Kitaev 's model. Those four GSs of the Kitaev model are 
indeed realized in Eq. (|92|) . 



multiples of a uniform constant gap. We have that 



As we remarked earlier, for a system with open bound- 
ary conditions, each Ag or Bp may be regarded as a de- 
coupled Ising spin. [On the square lattice with open 
boundary conditions, the star operator appears only if 
all of the 4 bonds which make it up appear in the lattice. 
We do not allow for star operators near the boundary of 
the system where only three or two operators appear in 
the product of Eq. pT|) defining the operator Ag.] The 
spins do not interact with one another and the resulting 
partition function is 



(ZA = lim 



d 



InZ 



= lim — ^ tanh(/3/ii); 

hz i^0+ hi 



(X,) = lim 



d 



■InZ 



= lim -^t&nh{f3hi). (109) 

In the above, we compute the finite-T expectation values 
of the Toric code operators (Xi), and (Zi) by inserting 
Eq. pOTj) into Eq. (fT09l) . We see that for all tempera- 
tures T > 0, 



(Zi) = {Z2) = (Xi) = {X2) = 0. 



(110) 



Thus, the existence of a gap in this system may not pro- 
tect a finite expectation value of the Toric code operators 
X12 or Z12 for any finite temperature. The expectation 
value of any local quantity is independent of the toric 
cycle operators. 

As we see here, by the use of our theorem and the 
results of [l7j , we can suggest possible symmetry allowed 
invariants. However, we cannot prove that these must be 
finite at T > 0. The Ising chain mapping further enables 
us to compute all symmetry invariant correlators which 
contain a finite number of fields 



^open = [cosh/?] 



M 



cosh /3/ii cosh /3/i2 , (108) 



with M given by Eq. ([57]) . The only meaningful ther- 
modynamic properties are given by the thermodynamic 
limit Ns 00 where, as they must, the results of Eqs. 
pOip . and p08p for periodic and open boundary condi- 
tions lead to identical conclusions. 

Alternatively, we briefiy note, that Kitaev's model can 
be identified as that of two decoupled Ising gauge theories 
on a square lattice with h = Q in Eq. . One of these 
gauge theories corresponds to the plaquettes Bp and the 
other to the sum of the star operators Ag. As the parti- 
tion function of the d ~ 2 Ising gauge theory is equivalent 
to the one of an Ising chain, the correspondence between 
Kitaev's Toric code model and the Ising chain immedi- 
ately follows. We now comment on a central observation 
in this Section: As the partition functions Z{f3) for D = 1 
Ising model and of the Kitaev Toric code model are the 
same, the Laplace transforms, the densities of states, are 
also the same in both systems. 

As a particular consequence of the spectral equivalence 
between the D = 1 Ising model and the Kitaev Toric 
code model, all non-degenerate states of the Kitaev Toric 
code model are separated from one another by integer 



= (tanh/3)ll^ll + ll^ll. (Ill) 



Here, [IS*]] is the number of star operators in the product 
of Eq. ()111|) . ]|P]| is the number of plaquette operators 
in the product. For an explicit proof of Eq. (jllip see 
our discussion in (6l| . The factorization manifest in Eq. 
(jllip illustrates that there are no connected correlation 
functions at any finite temperature. Although it is redun- 
dant for this simple case, in other systems with TQO, we 
may similarly employ overlap parameters and more com- 
plicated multi-particle correlators in order to discern a 
transition. 

As the result of Eq. (|110p is counter to some com- 
mon lore, we elaborate on it in a more general context 
and derive the same result from related general consider- 
ations. For any finite size system, the partition function 
is an analytic function (as it is the sum of a finite number 
of analytic functions). As the magnetizations (Xi), (Zi) 
are odd in the external fields yet must be analytic (and 
continuous) in the applied fields we trivially have that 
for any finite size torus at a finite inverse temperature /3, 

limh^_,-,o+{Z,{h,,P)) = liT[lh^,^o+{X,{h^.P)} = 0. As 
always, if SSB happens, it must occur in the thermody- 
namic limit. An alternate derivation of Eq. (|110p is as 
follows. Define on a system of size ^/TT^ x ^/TT^ (with 
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C2 



FIG. 10: The two contours Ci=i,2 which define the toric oper- 
ators of Eq. (|112p . The two contours Ci,2 are separated 
from each other by a distance y. The spins at the center of 
each link {ij) are aij. The shaded area denotes the region R 
[see Eq. (|113|l ]. The system is on a torus of horizontal and 
vertical cycles of size L = y/Wl. On each link ij there is an 
S = 1/2 field. 



Ns — > oo), the two topological operators 



(112) 



Here, we choose the two contours Ci_2 to be horizontal 
lines that encircle the torus (see Fig. [TU]) . The two con- 
tours are parallel to each other and are displaced relative 
to each other by a vertical distance (number of bonds) 
y. Let us denote the ribbon circumscribed between Ci 
and C2 by R (including Ci and C2). As for any bond 



[crfJ^ = 1 we have that 



PGR 

Now, by our Ising mapping [see [6l|], 
peR peR 



(113) 



(114) 



Here, ||-R|| denotes the number of plaquettes in R. Thus, 
(ZciZc,) = [tanh/3]^^. (115) 

Thus, [el 

lim(Zc,^c.) =0. (116) 

y— too 

This, in turn, suggests that 

(Zcu) = 0, (117) 



which is consistent with the result of Eqs. (|109p . and 
(jllOp . Even if these expectation values did not vanish 
(as they do not for general finite hi), for any finite /3, 
there is no ab initio reason for us to expect (^c^) to 
be the same for all contours {Ch} that encircle one of 
the toric cycles (e.g. the horizontal cycle). This is so as 
for non-zero temperatures As and Bp are not identical 
to one. This implies that, for instance, for two contours 
C and C which differ by the addition of a plaquette p, 
the expectation value {Zc') = (ZcBp) which is gener- 
ally different from (Zc)- Thus, at finite temperatures 
and fields, not only the topology of the contours Ci is 
important; their precise shape is also important. 

Insofar as the operators {Xi} and {Zi} are concerned, 
ergodicity remains unbroken at all finite temperatures. 
Generally, we can discern transitions by ergodicity break- 
ing at low temperatures which means that infinitesimal 
external fields lift ergodicity (states with positive magne- 
tization have higher weight even for infinitesimal external 
hi). Phase space is not sampled at low temperatures - 
ergodicity is broken and states with positive magnetiza- 
tion (or negative ones) for infinitesimal positive (nega- 
tive) fields are of larger weight. For the case at hand, 
just as in the D = I Ising model and free spins systems, 
ergodicity remains unbroken at all finite temperatures. 

While, obviously, Eq. (|110p holds for all positive tem- 
peratures, at T = different orders of limits can be taken 
in Eq. (jlOip [63t . The only way to have finite topolog- 
ical (Xi) or (Zi) is to keep them finite by force (apply 
external perturbations hi). In the Toric code, we proved 
that (Xi), (Zi) can remain finite at T > only if they 
arc kept finite by hand by for example coupling to a non- 
local external field. The simplest Toric code devised for 
pedagogical purposes works well at T = 0. It might have 
additional effects at T > 0, no matter how small T is. We 
reiterate our main result here: the gap is inconsequential 
in the context of keeping (Xi) and (Zi) finite. That these 
expectation values differ from their T = value relates 
to the non-analyticity of the free energy at T = 0. That 
is, it relates to the T = transition of the free spin (Ising 
chain) system. 

The physical engine for the absence of order in the 
Ising chains (and thus also in Kitaev's Toric code model) 
are topological defects (solitons) which proliferate and 
propagate freely at any finite temperature. The explicit 
form of these dual operators can be written to augment 
half of the duality transformations written in Eq. (|102p . 
The forms in Eq. (|102p appear as a duality for the spin 
for one component (say S^) while the other spin compo- 
nent (5"^) amounts to a string product (a soliton) on the 
dual lattice. In what follows, we will write the duality 
relations for the second half of these dual relations (that 
for S^) on the two decoupled Ising chains, e.g. [Ill, \Et\ . 



i<.s 

bp ^ IT' 

2<p 



(118) 
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Here, [as,a'^] = [bp,b'p\ = and {as,as} = {bp, bp} = 
0. The duality enables an identification of the d = 1 
symmetries and makes non-local terms appear local and 
vice versa. To make contact with our earlier discussion 
of Section [IIII regarding the unimportance of local events 
on the boundary in triggering the stabilization of TQO, 
we remark that the d = 1 operators of Eq. ([7]) are none 
other than the dual operators in Eq. (jllSp . 



2. Adiabaticity 

Our finite temperature solution allows us to demon- 
strate adiabaticity of the Kitaev model. That is, exact 
GLSs can be lifted while not changing the phase of the 
system. Let us re-examine the consequences of our solu- 
tion of Eq. (fTOO]) . The addition of, e.g., hx,i ^ to the 
Hamiltonian Hk does not lead to any transition within 
the system. Both the Kitaev model and its counterpart 
with any additional h^^i ^ lead to free energies that 
are everywhere analytic apart from T = 0. Nevertheless, 
insofar as GLSs are concerned augmenting Hk by an ad- 
ditional hx,i lifts the d = 1 symmetry Zi. Thus, the 
Kitaev model exhibits finite temperature adiabaticity. 

A similar phenomenon occurs in the adiabatic evolu- 
tion from a pure Ising gauge theory [Eq. (|84p which is 
equivalent to of Eq. with h = 0]] to another 

system. Consider the following change in Hamiltonians 

= H^,{K = 0)-KY,^Ij- (119) 

The 2,2 gauge theoryof Eq. ((84)) maps onto an Ising 
chain. By contrast [43|, iJzsC^a; = 0, /ij) maps onto a 
D = 2 Ising model in a finite magnetic field. Thus, 
the systems described by Hz^{hx ^ hz = 0) as well 
as 77^2 (^a: — 0, /12 ^ 0) exhibit no finite temperature 
singularities. Thus, -ffz2(^a: = hz = 0) and Hz^[hx = 
0, hz ^ 0) describe systems in the same phase. However, 
-Hz2 {hx = hz = 0) exhibits the local {d — 0) symmetries 
exemplified by {A^} of Eq. pT|) . as well as the d = 1 
symmetries {Xi, X2, Z2} of Eq. p^ . On the other 
hand, Hi^{hx = 0,hz ^ 0) only has the d ~ \ sym- 
metries of {^1,^2}. In other words, by setting hz 7^ 0, 
the d~ \ operators {Xx^X'i} arc no longer symmetries 
(and furthermore all local symmetries have been lifted). 
Nevertheless, even if a finite hz lifts d~ \ GLSs, no tran- 
sition occurs. Similar to the Kitaev model, in line with 
the lifting of symmetries, the application of a finite longi- 
tudinal field in a gauge theory \hz in Eq. (|119p ] lifts the 
GS degeneracy and leads to a single GS. Nevertheless, 
at any finite temperature, the systems are adiabatically 
linked to each other. This is analogous to the adiabatic 
continuity seen in matter coupled gauge theories (63 | as 
implied by effective magnetic field mappings [47| . 



S. Ising mappings of Kitaev models on general manifolds 

Here, we extend our Ising chain mappings to Kitaev 
models on general manifolds on closed oriented surfaces 
with genus or number of handles g. Kitaev's Toric code 
model on the torus (which has been the focus of our dis- 
cussion so far) corresponds to the case oi g = 1. 

To investigate what occurs on general manifolds, we 
apply the Eulcr-Lhuillier formula which states that the 
genus number g [or number of handles] satisfies 

V - E + F = 2{l- g). (120) 

Here, V is the number of vertices of the system, E the 
number of its edges, and F the number of its faces. In 
the Kitaev model, V denotes the number of sites of the 
system and as such counts the number of fields {^s}, E 
is the number of fields which live on the nearest-neighbor 
links (ij), and F is the number of plaqucttes (or number 
of terms {Bp}). For manifolds of genus number g > 1, 
the mapping to Ising chains [Eqs. PHH)) . piM|) . and pUS)) ] 
will proceed as before. The sole difference is that the 
number of edges (the number of fields {crij}) is larger 
than it was before. Here, Ne ~ Na + Nb + 2{g — 1), with 
Na^b the number of Ising spins on the two periodic chains 
[Na ^ V ^ Ns, Nb = F], and Ne is the number of fields 
{cTij}- We find that 

^^22(9-i)Zg=i. (121) 

In Eq. HHl), F = V = Ns and Zg=i the Toric code 
partition function which we found earlier [Eq. (jlOip 
for hi ~ 0]. [For general tessellation of a torus with 
F plaqucttes and V vertices, we will have ^g=i = 
( (2 cosh P)^ + i2 sinh /?) ^ ) ( (2 cosh /?) ^ (2 sinh /?) ^ ) .] Eq. 
p2ip is the Laplace transform of the density of states 
g{E). It implies that each energy eigenstates of the Ki- 
taev model on a manifold of genus g is precisely 2^^""^^ 
higher than that of two decoupled periodic Ising chains. 
In particular, the GS degeneracy is indeed [65| 

22(9-1) X 2^ = 2^9. (122) 

The precise exponential degeneracy of each level [the fac- 
tor 2^'^s~^)] holds for each level by comparison to its mul- 
tiplicity on the torus. 

B. Dimensional Reduction and Thermal Fragility 
in Wen's plaquette model 

Although appearing as different models, Wen's plaque- 
tte model of Eq. and Kitaev's Toric code model of 
Eqs. pT]|) and pT|) are one and the same. The equiv- 
alence between these two systems is demonstrated by a 
rotation which we will explicitly spell out later in this 
Section. As it must, according to this equivalence. Wen's 
model [1^ of Eq. ([T^ also exhibits a trivial dimensional 
reduction just as Kitaev's Toric code model docs. Here, 
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we have that with Ap — o'fa^a^crf the Hamiltonian of 
Eq. (HID, 

H^-KY^A, [^p,A']=0, Al = l. (123) 
p 

When periodic boundary conditions are appUcd (the sys- 
tem is on a torus), the number of plaquettes is equal 
to the number of vertices. As in the Kitaev Toric code 
model, the partition function of these decoupled plaque- 
tte Ising fields is, for an even by even lattice, given by 
the partition function to two decoupled Ising chains with 
periodic boundary conditions [ssj . In the absence of an 
external field that couples to the d = \ symmetries of 
this system, we have 

^even-even = [{2 COS\v P K)'^ ' ' + {2 sinh (3 K)^ ' ' . (124) 

For an odd by even lattice, the partition function is that 
of a single periodic Ising chain 

^even-odd = [(2 cosh /3X)^= + (2 sinh )^= ] . (125) 

The different decoupling for the even by even or even 
by odd cases follows from their constraints [l^. In both 
cases, we have that the product 

nA = l (126) 

p 

with the product over plaquettes p of the lattice. In the 
even by even lattice, we have the additional constraint 

n A = 1- (127) 

even p 

"Even p" refers here to dual lattice sites (plaquettes) in 
which the sum of the x and y coordinates is even. Here, 
the correspondence with Eqs. (|102p . (|104p . and ()105p can 
be made via 

■Ap'=2s ^ O'S, 

Ap'=2p+i ^ bp. (128) 

Similarly, as in Eq. (j95p . any GS can be expressed as 
a projection from a reference state: 

Nllil + ApW (129) 

p 

with 10) a reference state. We will soon rewrite the GS 
in an explicit form similar to that of Eq. (|92p . 

An explicit correspondence between the Kitaev Toric 
code model of Eqs. PU)) . and (|lip and Wen's model of 
Eq. is seen by simple rotations. This is illustrated 
in Fig. [TTJ To this end, we start with Kitaev's model of 
Eqs. ((To)) , and ((TT|) . We may draw the diagonal hues as 
shown in Fig. 1111 The spins which formerly lied on the 
bond of the lattice are now on the vertices of the rotated 
lattice. We may then identify every other center of the 




FIG. 11: Transformations relating Kitaev's Toric code model 
of Eqs. (fTOll . and ([11} to Wen's plaquette model of Eq. (fTi} . 
The two consecutive four spin terms correspond to a plaque- 
tte and star terms {Bp and As) in Kitaev's Toric code model 
of Eqs. (|10|) . and (lllf) . In the rotated lattice formed by the 
diagonals, Kitaev's Toric code model becomes Wen's plaque- 
tte model of Eq. (|14|) . (This occurs upon applying the spin 
rotation of Eq. (|130[) on one sublattice of the spins.) 

squares formed by these diagonals as a plaquette center 
(denoted in Fig. [TT] by a hollow circle) and the centers 
of the remaining squares by a star center (denoted by a 
full circle). Next, we rotate all of the spins (labeled in 
Fig. [TT] by the intersection points of the diagonals) on 
one sub-lattice by 180 degrees about the symmetric line 
in the x-z plane in the internal spin space, by 

dwen= n exp[z^«+a,^)]. (130) 

even sublattice 

In the rotated system formed by the diagonals, what we 
obtain is Wen's model of Eq. ([T^ [with a relabeling of the 
spin polarizations]. Due to their exact equivalence, all of 
the conclusions stated for the Kitaev model in Section 
Kill Al will hold for Wen's model as well. In particular, 
the GS for a system (of A' > in Eq. ([H])) with pe- 
riodic boundary conditions along the diagonals is given 
by Eq. ([5^ after undoing the various rotations (in both 
spin and real space), 

|^,,)=2i-^.'/2^|c'). (131) 

Here, q' = (±1,±1) denotes the transformed topological 
sectors. That is, replacing the topological numbers ±1 
for the eigenvalues of the d = 1 symmetries Zi=i,2 in Ki- 
taev's model (Eq. ([T^ ) which now refer to the direction 
relative to the diagonals in Fig. Ill[ we have 

2^ = &WenZ^O^en- (132) 
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Similarly, the states |c') denote all states in which = 
±1 for all sites a of the even sublattice and cr^ = ±1 for 
all sites b of the odd sublattice for which the product of 
these quantities around any plaquette is one. Following 
the same considerations as for Kitaev's Toric code model 
which we outlined in the previous Section, we see that the 
d = 1 symmetry operators attain a vanishing expectation 
value at any finite temperature unless a finite external 
field is applied. 



C. Finite-temperature transitions in related 
higher-dimensional analogues 

As we proved above, Kitaev's Toric code model as well 
as Wen's related plaquette model do not exhibit a finite- 
T phase transition by virtue of their mapping onto a 
system of Ising chains which exhibit a T = transition. 
It is obviously the case that other systems with TQO can 
exhibit finite-T transitions. For instance, the D = 3 Ising 
gauge theory exhibits rank-7i = 8 TQO yet, as can be 
seen by the fact that it is dual to a D ~ 3 Ising model, it 
exhibits a singularity in the free energy at a finite critical 
temperature (Tc). A related system is a D = 2 extension 
of the Kitaev Toric code model of Eqs. (fTU]) . and (fTTI) 
[6^. Here, the product in the star operator Ag of Eq. 
PT|) spans all six sites which, on the cubic lattice, are 
the nearest-neighbors of a given lattice site. Similarly, 
in the product defining the plaquette operator Bp of Eq. 
(fTTj) . we have (as in the D = 2 case), the product of cr^ 
on all four bonds which form a planar plaquette (parallel 
to either the xy, xz, or yz planes). It is easy to prove 
that the Hamiltonian that results is none other than that 
of the D = 3 Ising gauge theory (given by — J^p Bp) 
augmented by the sum of all local symmetry generators 
(— As). It is Eq. ((86|) with, in the notation of Section 
IXIIl H being the D = 3 Ising gauge theory Hamiltonian 
[i.e. Eq. ((84)) with the sum performed over all plaquettes 
that appear in the cubic lattice] and now with {Gi] the 
set of all local symmetries ({^s}) 

B H^£) Ising gauge "t" -^i;eriea; 

Ising gauge ^ ^ Bp , HyQj^iQx ^ ^ As . (133) 

p s 

Unlike the D = 2 case, in three dimensions, the pla- 
quette terms {Bp} which appear in Eq. p33p satisfy 
too many constraints that replace the single constraint 
of Eq. For instance, if we consider a single cube on 

the lattice which is of unit length and label its six faces 
by Cubep=i_2.- - ,6 then wc have that IlpeCube^P ~ 
Precisely all of these constraints (and none additional) 
appear in the D = 3 Ising gauge theory. All terms in Eq. 
()133p commute and as no constraints couple the vertex 
fields to the plaquette fields, the partition function is the 
product of the decoupled plaquette and vertex fields, 

^ — Ising gauge ^ ^ID Ising chain- (-^'^^) 



The first term in Eq. (|134p arises from the plaquette 
terms whereas the second term is that from the vertex 
(star) fields {As} which, as in the D = 2 case, adhere 
to the single constraint they satisfied in Eq. The 
counting of the degrees of freedom proceeds as follows: 

Z3D Ising gauge IS thc traCC of exp[—(3H3D Ising gauge] 

over all gauge independent degrees of freedom whereas 

ZlD Ising chain is thc traCC of C-Kp[- (3H^ertex] OVCr all of 

the gauge degrees of freedom {Gi} (which are here exem- 
plified by the star terms {A^}). In the thermodynamic 
limit, for the system with both periodic or open bound- 
ary conditions, the free energy per site is given by 

F = Fi(/3) + F3 (/?*). (135) 

Fi is the free energy per site of an Ising chain; this con- 
tribution is borne by the star fields {A^}. Similarly, F3 
is the free energy per site of a Z? = 3 Ising model at the 
dual /?* [eS]: 

sinh2/?sinh2/?* = 1. (136) 

The contribution of F3 originates from the free energy of 
the Ising gauge theory {{Bp}) which is exactly dual to 
the D = 3 Ising model. The free energy of Eq. p35p 
exhibits a T = (/3ci = 00) singularity (borne by the 
T ~ singularity of the D = 1 Ising contribution Fi) as 
well as a finitc-T singularity, 

Pc2 ^ 0.761423, (137) 

which originates from the D = 3 Ising contribution F3 . 

As in the D = 2 case, no toric operators can attain a 
finite expectation value at any finite temperature, e.g. 

{ZpJ^i n <>=0, (138) 

with Pa an entire plane. The steps of Eqs. (|112p . (|113p . 
()116|) . and pi7p may be reproduced with the former lin- 
ear contours {Ci,2} replaced by planes Pa. Due to thc de- 
coupHng of the plaquette {{Bp}) and vertex ({A^}) oper- 
ators, the expectation value on the righthand side of Eq. 
(|113|) . with the contours {Ci^2} replaced by two planes 
{Pi, 2}, is given by its value for a classical D ^ 3 Ising 
gauge theory. As seen by high and low coupling expan- 
sions [isj now applied to surface operators, the correlator 

{Zp.Zp,) = ((riyePi <^ij){UijeP2 '^u)) vanishes in the 
D = 3 Ising gauge theory as the boimding planes Pi^2 
are taken to be infinite. In both (i) the confined phase 
[P < Pc2\i a-iid in (ii) thc deconfined phase [/3 > /3c2]i the 
pair correlator (Zp-^Zp^) vanishes exponentially in the 
total area of thc planes Pi and P2 . 

For general dimension D, the partition function asso- 
ciated with a simple extension of Eq. (jl33p is given by 

Z Z jj^dijyiQjisional Ising gauge ^ Z\pi Ising chain -(139) 

In the two dimensional case, which we discussed earlier, 
the first term- thc partition function of the two dimen- 
sional Ising gauge theory- is nothing but the partition 
function of an Ising chain. 
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XIV. SHAPE AND STATE DEPENDENT 
ENTANGLEMENT ENTROPY: RESULTS AND A 
CONJECTURE 

With py the marginal density matrix obtained by trac- 
ing out all degrees of freedom exterior to a boundary sur- 
face r (of size i?), the entanglement entropy between the 
two regions which are separated by T is defined as 



S'ei 



-tr [pY Inpr]- 



(140) 



Recent work computed the entanglement entropy for a 
variety of systems of pertinence, e.g., AKLT type integer 
spin chains [6^ . The quantity Sent measures the entan- 
glement between fields in the interior of F to those lying 
outside. It was recently conjectured [g^ that a non-local 
borne deviation ^ > from an asymptotic area law scal- 
ing. 



(141) 



iTL D = 2, constitutes a precise measure of TQO. In Q we 
conjectured that in some systems, this criteria is either 
incorrect or, at best, may require an additional explicit 
average and/or limit, of a form yet to be prescribed, over 
all allowed contours as the size of the contour tends to in- 
finity. The central point of Q was that the entanglement 
entropy may depend on the shape of T [or, equivalently, 
that for a given shape we may find GSs for which the 
bare entanglement entropy criteria is inconsistent with 
the definition of TQO given by Eq. ([1])] . These deviations 
from an area law are compounded by other deviations in 
other arenas (e.g. those for tight binding fermions [70| ) 
which are unrelated to any appearance of TQO. We now 
elaborate on this. 

In Klein spin models which host TQO (by satisfying 
Eq. (P)) - e.g. arbitrary rank T > TQO of the Klein 
model on the pyrochlore lattice [s^ the entanglement 
entropy criterion is not satisfied. In a pure dimer state- 
a state which is a product of individual singlet pairs, 
all GS correlations are manifest as local singlets and, at 
T = 0, the entanglement entropy S'ent is, up to constants, 
the number of dimers partitioned by F: 



Qpure 
"^ent 



Ngr In 2. 



(142) 



Here, Ngr is the number of dimers which are partitioned 
by F: the number of singlet dimers in which one spin 
forming the singlet lies inside F and the other lies outside 
F. A random surface F intersects OdlFjl) dimers, with 
||F|| the number of dual lattice sites which lie on F. How- 
ever, by setting the contour F along the dimers present 
in any given GS, we find arbitrarily large surfaces such 
that no singlet correlations are disrupted and S'ent = 0. 
We similarly find many GSs/contours with arbitrary ^ for 
fixed contours F. The arbitrariness of the entanglement 
entropy and of the deviation ^ in systems with arbitrary 
rank - TQO suggests that the conjecture raised by [6^, 
although very insightful and correct in some cases, may 
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FIG. 12: Two singlet product states which are GSs of the 
Klein model Hamiltonian of Eq. (|34|l on a checkerboard 
lattice. The ovals denote singlet states formed by nearest- 
neighbor spins. The dashed line denotes a contour P. The 
entanglement entropy between the two regions separated by 
F is different for the two above states: it is (5 In 2) for the 
singlet product state on the left while it is In 2 for the state 
on the right. Both of these singlet product states can be dif- 
ferentiated by local observables and thus do not satisfy the 
TQO conditions of Eq. 



be inconclusive. Similar conclusions also appear for the 
more sophisticated subtractions [termed a topological en- 
tropy] which were suggested by [6^ for the extraction of 

In what follows, we elaborate on this conjecture. We 
start with (i) a rigorous demonstration that non-TQO 
states or states that, at best, can capture arbitrary rank 
TQO (yet do not satisfy Eq. (P) for all GSs) in which 
the entanglement entropy can assume any value between 
zero and a quantity which scales linearly with the area 
||F||, (ii) Show that also the topological entropy can sim- 
ilarly assume any value (and in particular) be finite in 
systems which do not host TQO, conjecture on (iii) the 
existence of specific TQO states in the Klein model whose 
TQO is manifest not necessarily in entanglement entropy 
but rather in topological invariants (line number) that 
classify sectors of the GS subspace, and finally construct 
states [not necessarily only GSs of the Klein Hamiltonian] 
which have TQO yet display vanishing entanglement en- 
tropy. 



A. Pure singlet product states 

Pure singlet product states |c) - states which are a 
product of disjoint singlets - do not support TQO (i.e. 
there the local operators - [Si ■ Sj] on nearest-neighbor 
sites ij which distinguish the states). 



1. Short-ranged SU{2) Klem models on a lattice geometry 

Spin 5=1/2 [SU[2)\ Klein models were briefly re- 
viewed in Section IVIII Here, dimer states can support 
arbitrary rank TQO yet do not support Eq. ([T]) for all 
GSs of the Klein model. Here, various contour shapes 
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r can be tailored to lead to any value of the entangle- 
ment entropy (from zero up to an upper bound which 
scales with the area ||r||). In a dimer system, nearly any 
value of the entanglement entropy bounded from above 
by ||r|| times a constant is possible (see Fig. [T2|) . By 
setting r to lie along dimers, the entanglement entropy 
strictly vanishes. Similarly, in the ferroelectric GS (see 
[36j). if we choose the boundary T (which circumscribes 
an inner domain including all sites on F itself) then the 
number of partitioned dimers for a boundary F which is 
of a square form and lying along the non-diagonal axes of 
the checkerboard lattice) is exactly half the perimeter. In 
fact, for any given dimer state, we can find a multitude 
of contours F of arbitrary size such that 



< ^e^nt 



ure dimer 



(F) < ||r|| ln2. 



(143) 



[For an elementary derivation of the factor of In 2, see 
SectionlFl] In particular, it is not a problem to construct 
dimer states for a given contour F - or equivalently con- 
struct contours F for a given dimer state - in which, any 
desired function of the linear scale of F is possible so long 
as it satisfies Eq. (|143p . To conclude, we illustrated that 
systems exist where the entanglement entropy can as- 
sume any value between zero and a quantity which scales 
linearly with the size of the contour ||F||. In the singlet 
product states under consideration, Eq. ([1]) does not hold 
and the system does not exhibit TQO. This arbitrary de- 
viation from an area law scaling does not conform with 
either a pure area law scaling anticipated for non-TQO 
systems nor to a unique topological deviation ^. 

It should be noted that here the topological entropy 
Stoim as defined by is 



Stopo — Sa + Sb + Sc ~ Sab — Sac ~ Sec 

+Sabc 



for a partition of a large contour F which spans a region 
[ABC] into three separate regions A, B, and C gives a re- 



sult which is identically zero: St 



topo 



0. The motivation 



for introducing the Stopo of Eq. (|144[) in [68| was that in 
many cases, Stopo = C of Eq. (|14ip . 



2. SU (N > 4J Klein models on a small world network 

The basic idea underlying the Klein models can be ex- 
tended to SU{N > 2) systems. This can be done for 
any lattice/graph in which we can choose N sites within 
each fundamental unit (each "plaquette" ) to belong to an 
SU{N) singlet. Obviously, for > 4, the basic "plaque- 
tte" spans more 4 sites. The covering of the lattice/graph 
with these singlets is that with hard core units (singlets) 
each of which spans A^ sites. These Klein models are 
what we will now introduce for the first time. Here, we 
will have 



Stat 



(145) 



with J > 0, and pStot=max ^j^g projection operator onto 
SU{N) spin states on a given plaquette a which belong 
to the highest irreducible representation of SU{N). Now, 
let us consider the Klein model on a small worlds net- 
work. In a small worlds network [69j , most links connect 
sites which are geometrically close to each other yet once 
in a while there is a link which connects very far sepa- 
rated sites. Such a small worlds geometry allows for the 
existence of singlet states composed of far separated sites. 

If we consider a long-ranged SU{A) singlet (see Fig. [T3|) 
which is composed of a single site in each of the 4 regions 
A, B, C, and the region exterior to ABC, then Eq. p44p 
reads 



-ilong singlet 
^topo 



ASl' 



(146) 



Here, S*™* is the entanglement entropy between one site 
and the remaining three sites of an SU{A) singlet. Sim- 
ilarly, S"!"* denotes the entanglement entropy between 
the two pairs of sites which constitute an SU{4) singlet. 
Here, we find that for an SU{A) Klein model on a small- 
worlds network [where the plaquettes can contain sites 
which are arbitrarily close or far apart] the entanglement 
entropy is equal to 



Stopo — A^/ong[4iS'™* — 2>S^*]. 



(147) 



Here, Niong denotes the number of singlets in which each 
of the four sites lies in a different region. We have that 



In 4, 



ST* = ln6. 



Inserting this in Eq. (|147p leads to 
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Stopo — Ni 



long 



In 



27 



(144) For general SU{N) systems (fc = 1, • • • ), 



ST* = In 



(148) 
(149) 

(150) 



A detailed derivation of these combinatorial forms is 
given in Appendix [F] The entanglement entropies asso- 
ciated with the completely antisymmetric SU [N) singlet 
state arc equivalent to those of an A^ fcrmion system in 
which no interactions are present. 

Returning to Eq. (|149|) we see that in general, for 
Niong 7^ in an S[/(4) Klein model on a small worlds 
network we have 



Stopo 7^ 0- 



(151) 



Nevertheless, whenever singlets are formed between 
nearest- neighbor quartets, we have a local operator 
[J2i=i S'i]'^ = 0. When no singlet is formed by these 4 
points, {[J2i=i Si]"^) > 0. We thus, in general, have a lo- 
cal operator ([X]i=i Si]'^) which attains a different expec- 
tation value in different GSs - the system does not pos- 
sess strict TQO. Nevertheless, here, in general, Stopo 7^ 0. 
Furthermore, here, as before, Stopo is quite arbitrary. 
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FIG. 13: Shown above is a singlet state formed by distant sites 
in the regions A, B, C, and the region exterior to these. The 
lines delineate the members of this 4 site clique. An S!7(4) 
Klein model Hamiltonian on a small worlds network will in- 
clude a vast number of short-range singlets yet also allow for 
singlet states between remote spins. The topological entan- 
glement entropy of Eqs. (|144|l , and (|146|) which is associated 

with the singlet state shown above is Stopo — In fy ■ 



In principle, for any scheme for Stopo similar to 
Eq. (|144p which requires partitions into (i?-l) regions, 
an SU{N > R) model on a small worlds network can 
generally give rise to an arbitrary Sfopo- In this system, 
general local order parameters can differentiate different 
pure [non-superposed] singlet product states. 



B. TQO ground states in the nearest-neighbor 
5=1/2 Klein model on the checkerboard lattice 

We can construct a set of states which satisfy Eq. (P) 
in which we conjecture that a direct topological, rather 
than entanglement entropy, marker might be the most 
transparent. As stated earlier, in individual dimer states, 
local observables can differentiate different GSs. For in- 
stance, if a dimer connects the two sites i and j then 
{Si ■ Sj) = —3/4. As such local measurements attain 
different values in different GSs, individual dimer states 
do not exhibit TQO at T = 0. Linear combinations of 
intra-plaquctte dimer states [H] exhaust all of the GSs 
of the Klein model on any lattice and may enable the 
construction of TQO states. Let us focus on the Klein 
model of Eq. p4|) on a checkerboard lattice of size Lx L. 
This Hamiltonian has GSs given by Eq. ([55)1 . Now, let 
us focus on a subset of these GSs. To this end, let us 
draw a square contour F of edge whose length is R and 



consider the states 

E 1^)- (152) 

|c) has I lines, no dimers partitioned by F 

Here, \ipi) is defined for each individual number / of lines 
in the line representation of the dimer coverings (see d 
of Section IVIII (Figs. [71 and [5] in particular)). We may 
focus on the physical regime defined by (L/2 + R) > 
I > {L/2 - R) and L R. In Eq. p52)) . the states 
\ipi) are sums over all dimer states |c) which belong to a 
very specific chosen sector. The multiplicative constant 
A/" is a normalization factor. Here, the last constraint 
amounts to the requirement that each dimer configura- 
tion summed over exhibits a dimer covering GS within 
the square contour of size R defined by F with no spins 
unpaired. A/" is a normalization factor. Clearly, there 
are no local processes linking the different As 
quasi-local {d = 0) symmetries [sums of local permuta- 
tions] link all dimer states within a given sector of / lines 
to each other, the expectation value of any quasi-local 
quantity, {ipilVltjji), cannot be used to differentiate the 
different states which are given by Eq. p52p and Eq. ([T|) 
may follow. Here, the different states and their charac- 
ter are clearly marked by the topologically conserved line 
number I appropriate to each sector. As no dimers cross 
F, the T = entanglement entropy across F may be con- 
jectured to vanish, S^^^ = for each of the states | ?/';). 
The contour F may be taken to be of an arbitrary shape. 
The states of Eq. (|152p form the analogue of the states of 
highest spin in SU{2) discussed in Section IXI Bl Regard- 
less of the form of the entanglement entropy in these and 
other cases, the topologically conserved quantities (the 
line number here) may be the simplest marker of TQO. 



C. Vanishing topological entanglement entropy in 
states with TQO 

The above conjecture can be made rigorous in a some- 
what more restricted set of states. These states will be 
constructed in a manner similar to the TQO GSs of the 
Klein model of Eq. (|152p . Nevertheless, these states need 
not be GSs of the 1/2 SU{2) Klein Hamiltonian of 
IVIII To this end, let us first note that Eq. ^ implies 
that if T = TQO appears in a set of states {</>" }a=i 
any given a = 1, 2, • • • , iVg then it appears in their prod- 
uct state. That is, let us consider a set of "fragmented" 
states {</>a}a=i ea-ch of which satisfies Eq. ^ and each of 
which has its support in a region of the lattice [or volume] 
Aa such that 

Aa n A;, = for all a ^ b, and 

AiUAaU-'-UAt = A. (153) 

In Eq. (|153p . A denotes the entire lattice [or volume] of 
the system in question. Let us next construct the product 
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states 



(154) 



which has the entire lattice [or vohime] A as its support. 
The reader can easily convince her /himself that if Eq. ^ 
holds for all states {0a}a=i of Eq. ()153|1 then the states 
of Eq. (|154p satisfy Eq. ([T]) as well. In other words, these 
product states display T = TQO. 

Next, let us consider a lattice partitioned into four re- 
gions Aa with a = 1,2, 3, 4. In the notation of Fig. (fT2|) 
these correspond to the four regions A, B, C, and the 
remainder of A which we henceforth label by D. Now, let 
us construct the following fragmented states 



\(pa;lj =M'a. 



c) has la lines in A^ 



(155) 



By la lines we allude to the number of lines piercing a 
horizontal plane as it bisects region a. Next, we examine 
the states 



|(^C;Lc) ® \(t>D:Lo 



(156) 



As in Eq. (|154p . these states display TQO. A calcula- 
tion of the topological entanglement entropy of Eq. (|144p 
reveals that in each of these states Stopo = 0. This occurs 
notwithstanding that each of these states exhibit TQO. 



D. A Summary and a new conjecture 

In conclusion, the considerations of the singlet product 
states of Section IXIV Al suggest that if the conjecture of 
[6^ indeed holds universally on lattice systems, an appro- 
priate averaging and/or limit process needs to be spelled 
out for different contours F such as those detailed here 
for, at least, some lattice systems. Otherwise, the en- 
tanglement entropy measure of TQO is inconsistent with 
the original definition of TQO [Eq. ([1])] as represent- 
ing a robustness to local perturbations. We raise here 
the conjecture that if topological entanglement entropy 
is to capture TQO then, at the very least, Stopo needs 
to be contour independent. This leads to an exponen- 
tially large - in the system size - number of conditions 
(the topological entanglement entropy for each of the ex- 
ponentially large number of partitions must be set to a 
fixed value). 



XV. GRAPH PROBLEMS: ENTROPY, 
SPECTRA, AND TOPOLOGY 

The insufficiencies of spectra and entropy in determin- 
ing if TQO is present have counterparts in the topology 
of graphs and in the Graph Equivalence Problem (GEP) 



in particular. The adjacency matrix of a graph has ele- 
ments Cij = 1 if vertices i and j are linked by an edge 
and Cij = otherwise. Vertex relabeling i p(i) leaves 
a graph invariant but changes the adjacency matrix C ac- 
cording to C — *■ C" = P^CP with P an orthogonal matrix 
which represents the permutation: P = <5j,p(i). The GEP 
is the following [tI]: "Given C and C , can we decide 
if both correspond to the same topological graph?" The 
spectra of C and C" are insufficient criteria. Entropic 
measures (72| are useful but also do not suffice. Similar 



to the current ideas in the study of TQO, both measures 
(spectra and entropy) were long ago suggested as a way 
to partially flesh out topological structures in graphs. In 
order to make this connection more rigid, we discuss in 
the Appendix new Gauge- Graph Wavefunctions - states 
defined on a graph in which TQO is manifest. 



XVI. STRING AND BRANE TYPE 
CORRELATORS 

Hamiltonians potentially capturing a few aspects of the 
properties of novel materials display non-local string or- 
ders [zl, M M M [23, M M, m MMM- Here, 

there are non-local correlators which display enhanced 
(or maximal) correlations vis-a-vis standard two-point 
correlation functions. In this way, the presence of off- 
diagonal long-range order (ODLRO) in those string cor- 
relators does not imply the existence of TQO as defined 
by Eq. ([T]). However, this ODLRO in the string/brane 
correlators signals the presence of a hidden order of a non- 
local character with important physical consequences. 

It is important to emphasize that these string correla- 
tors are, generally, different from those in gauge theories 
or more generally systems with exact low c?-GLSs. To 
make the distinction clearer, we reiterate the symmetry 
considerations of Section IIX CI A t finite temperatures, 
by Elitzur's theorem for d = [4:^ theories and by its ex- 
tension to d = 1 discrete or d < 2 GLSs [l3| , all non GLSs 
invariant quantities must vanish. Only GLSs invariant 
quantities can be finite. In the case of gauge theories, as 
well as in other theories, this only allows for string corre- 
lators (closed Wilson loops or open meson-like lines) to 
be finite. At T = 0, similar considerations often follow 
while working in the eigenbasis in which the d-GLSs are 
diagonalized: in this basis, only d-GLSs invariant quan- 
tities can attain a finite expectation value. For a d-GLSs 
U we have that {ga\U\ga) ^ Ua- If C/ is a d = 1 GLS op- 
erator then we will also have another trivial example of 
a string correlator. For higher d > 1 GLSs, we will have 
brane like correlators. These considerations are different 
from the maximal (or enhanced) string correlators of for 
example Refs. Ji, 75., 7g, 7_7, 7&, 79,, M, .81, 82, sj 
in which string correlators are not mandated by an ex- 
act symmetry. We will come back and make contact and 
comparison with these considerations towards the end 
of this Section. We will suggest an extension to higher- 
dimensional brane correlators, detailed examples of which 
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will be provided in future work. 



A. An algorithm for string/brane correlators which 
relies on ground state selection rules and non-local 
transformations 



of the total spin on sites j and (j + 1) onto a value of 



2, i.e. P 



Stot=2 
3,3 + 1 ■ 



This allows for the construction 



of GSs [8JI which are of the form 



(159) 



We now outline an algorithm for the construction of 
such non-local correlators which are not mandated by 
symmetries. (In systems with uniform global order al- 
ready present in their GS, the algorithm leads to the 
usual two-point correlators.) We seek a unitary trans- 
formation Us which rotates the GSs into a new set of 
states which have greater correlations as measured by a 
set of local operators {Vi}. These new states may have 
an appropriately defined polarization (eigenvalues {vi}) 
of either (i) more slowly decaying (algebraic or other) 
correlations, (ii) a uniform sign [partial polarization) or 
(iii) maximal expectation values Vi = Vmax for all i [max- 
imal polarization) . Cases (i) or (ii) may lead to an emer- 
gent lower-dimensional gauge-like structure for the en- 
hanced correlator. Generally, in any system with known 
(or engineered) GSs, we may explicitly construct polar- 
izing transformations Us- In some cases, these operators 
Us embody a local gauge-like structure. Our contention is 
that these unitary operators which rotate the (generally) 
entangled states onto uniformly states of high polariza- 
tion embody an operator P\o^ for which 



(157) 



for any state \g) in the GS manifold. In the AKLT spin 
chain which we will focus on shortly, there is associated 
with each pair of sites an operator for which the se- 
lection rule of Eq. (|157p is satisfied. 



1. Selection rules and polarizing transformations in the 
AKLT spin chain 

In what follows, we illustrate that in the presence of 
GS entanglement characteristic (although not exclusive) 
to TQO, maximal (or in more general instances finite) 
amplitudes of correlation functions evaluated within a 
certain GS are seen with a non-local string correlator 
borne by a non-trivial polarizing operator Us ■ To provide 
a concise known example where these concepts become 
clear, we focus on case (ii) within the well-studied S = 1 
AKLT Hamiltonian [H 



AKLT — 



Sj ■ Sj+i + - [S-j ■ Sj+iY 



n 

3 = 1 



2l 



(158) 



for a chain with Ns sites. The argument of the summand 
is, up to constants, nothing but the projection operator 



In Eq. (UnSl), each 5 = 1 is replaced by two S" = 1/2 
spins and singlets are formed between of these nearest- 
neighbor S" = 1/2 spins, finally the 5' = 1/2 basis state 
is projected onto the S = 1 basis (which is attained by 
the symmetric projection operator Piy^nm) Hi- Conse- 
quently, the total spin of any nearest-neighbor pair can 
at most be 1. Within such a state. 



(160) 



a realization of Eq. (|157p . In a system with open bound- 
ary conditions, there are 4 degenerate GSs in which there 
remain fractionalized spins at the two endpoints of the 
chain. 

In what follows, we discuss gauge-like aspects of this 
system. (Details underlying the assertions below are 
given in AppendixlGl) Here, there is a non-trivial unitary 
operator 

Us = Y[ exp[i7r5|^^], with [Us,H/^klt] + 0, (161) 

j<k 

which maps the GSs into linear superpositions of 

states in each of which the local staggered magnetization 
Vj- = [—\yS^^ is uniformly non- negative (or non-positive) 
at every site j. As all transformed states 



\Va) Us\gcy) 



(162) 



are superpositions of states with uniform sign Vj (allow- 
ing for two non-negative or non-positive Vj values at every 
site out of the three 5=1 states) , the correlator 

= {pa^y^y^Vc) = {gc.\UtV,V,Us\g^) 

= (-l)l*--'"l(.ga|Sfg,y5||5a) (163) 

can be computed to give \Gij\ ~ (2/3)^ for arbitrarily 
large separation \i — j|, i.e. Gij displays ODLRO. Here, 
Qij is a shorthand for the string product 



i<k<j 

Given that the two-point spin correlators 



G, 



{s-s- 



(164) 



(165) 



decay algebraically, we have, for \i — j\ ^ 1, 

\{S!Q^,S])\:^\{StS])\. (166) 
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The inequality of Eq. p66p was shown to hold in a 
large family of D = 1 5' = 1 Hamiltonians which ex- 
tends the AKLT point [sHi ■ In Appendix \G\ we pro- 
vide a proof of Eq. (|166p that is based on classes of 
states and not specific Hamiltonians. Similar large am- 
plitude string or non-local correlators appear in other 
arenas [zl, M, M, M, IzS, H S [H, [H. Thus far, 
systems hosting such correlations (and the form of the 
correlations themselves) were found on an example by 
example basis. In what follows, we suggest how exact 
inequalities of the form of Eq. (|166p may be derived in a 
systematic fashion. To couch our discussion in the sim- 
plest possible terms, the terminology will often follow is 
that of the spin systems (and of the AKLT Hamiltonian, 
in particular). Nevertheless, the ideas introduced below 
allow a construction of non-local string operators in any 
instance in which the GSs of a certain system are known 
and in which these GSs are entangled. Entanglement im- 
plies that the operator Ug which maps the GSs into a set 
of uniform (polarized) states {|pa)} cannot be a uniform 
product of local operators (or of sums thereof). That is, 

Us^li 0„ (167) 

with {Oj} local operators. 

This unitary transformation generalizes the exact sym- 
metry of the GSs which we discussed in earlier Sec- 
tions. Of all 3^= linearly independent basis states 
{\Sf,Si,--- ,S'w,>} only (2^»+i - 1) might have a fi- 
nite overlap with a GS: In the AKLT example, all GSs 
are superpositions of Neel-like segments of the chain 
(with one of the two sublattice parities: Sf = (—1)* or 
Sf = (—1)*+^) [m which are separated by sites at which 
5f = 0, e.g. 

I + - + -0+-0 +•••)• (168) 

The rule of thumb for determining the sublattice parity of 
the Neel orders across the zero sites is as follows: opposite 
parity Neel orders exist across an odd number of zeros (as 
in Eq. (|168p ). and similar ones exist across even length 
numbers of zeros. We find that this structure can be 
captured by the following d ~ operators within the GS 
basis, 

P^j = ~S^Q^,S^ + (1 - iS^S^j'), (169) 

for which 

P^M = \g) (170) 

for all GSs \g) and for all site pairs (?,j). That Eqs. 
P69p . and p70p are satisfied follows from the form of 
Eq. (HnH). The wavefunction of Eq. - that of 

consecutive singlets - does not allow the total S^' on any 
string to have an absolute value that exceeds 1. That is, 

|(5lE^^l5)l<l. (171) 

l=i 



It is worth to notice that Eqs. (|169p . and (|170p imply 

{g\S^Q,,S^\g) = -{g\{S^S^r\g). (172) 

That is, the hidden non-local order is equal to a local 
nematic-type correlator. ( As we will see below, this 
is true for any state that belongs to the subspace Ti-y) 
In particular, the fact that the string correlator Gij of 
Eq. (|163p is equal to (2/3)^ implies that, within the 
GS subspace, ((S'f)^) = 2/3 at all sites i. The opera- 
tors of Eq. p69p act as the identity operator within the 
GS basis: they form a realization of our general rule of 
thumb of Eq. (|157p . They should not be confused with 
obeying exact d — GLSs selection rules which appear 
throughout the entire spectrum. In Fig. 1141 we schemat- 
ically represent the set of all states {lya}} which sat- 
isfy Pijlua) = lUa)- This set includes the GS basis as a 
very special subset. The set of independent basis vectors 
IIj/q)} may be equivalently defined by finding all such 
basis states such that the target space Us\ya) degener- 
ates into partially polarized states Us\ya) = \Pa)- In the 
AKLT system, these partially polarized states are those 
with non-negative (non-positive) entries, e.g. 

I + + + +0 + +0 + •••)• (173) 

Equation p73p is the state of Eq. (|168p after it would 
be acted upon by the partially polarizing transformation 
operator of Eq. (|16ip . The polarization transformations 
are just another way of stating the selection rules of Eqs. 
(Unni), and p7D)) . The operator of Eq. pBT]) leads to a 
partially polarized state [of the form of Eq. (|173p ] for 
all states which satisfy the selection rules of Eq. (|170p 
[e.g. the basis state of Eq. p68p ]. We denote by Hy the 
Hilbert subspace spanned by {jya}}. By unitarity, Hy 
spans (2^=+^ — 1) states. This is so as the number of 
partially polarized states {jpa)} is (2^=+^ — 1). Further- 
more, let us denote by Ty the set of operations which have 
finite matrix elements only within this basis. These op- 
erators (albeit not symmetries of the Hamiltonian) con- 
stitute analogues of the (d = 0) gauge transformation 
symmetry operators of Z2 gauge theories linking degen- 
erate states. These d ~ Z2 transformations correspond 
(in the ^J^i Sj basis) to the creation (annihilation) of 
an Sj — state at any site j followed by a unit displace- 
ment of at all sites k > j. 

In the AKLT Hamiltonian, the entire set of basis vec- 
tors {Ij/a)} exhibits an invariance under local Ising gauge- 
like operations (Z2 (g) Z2 (S) • • • (81 Z2 with a Z2 appearing 
for each of the j = 1,2, ■ ■ ■ , sites) which keep the non- 
negative (non-positive) entries of {—lySj non-negative 
(non-positive) . Gauge- like operators enumerate the num- 
ber of independent rays {|?/a)} (the size of the Hilbert 
space in which the GSs exist in the ^^^i \Sj) basis): 
(2^-'+^ — 1) with, asymptotically, a factor of 2 given by 
an Ising (Z2) degeneracy present at each of the Ng sites. 

The structure which emerges is not that of the GSs 
(which is that of a global Z2 (8) Z2 symmetry) but rather 
that of the S'^ basis in which the GSs reside. With 
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{f^GSs} the set of unitary operators which have their 
support only within the GS basis, we have the obvious 
{Ty} D {?7gss} with the subscript referring to the set 
of symmetries which leave the larger GSs basis support 
Ij/q) {local gauge in the AKLT case) or the smaller set of 
GSs only un-admixed with other states outside that basis 
or GS manifold respectively. With generalized polarized 
states giving rise to non-local string correlators, we gen- 
eralize the notion of unitary transformations associated 
with the manifold of the GSs to the larger set of trans- 
formations living on the manifold of states which satisfy 
the GS selection rule of Eq. (|170p (of which the GSs are, 
by definition, a subset). This is schematically illustrated 
in Fig. [H 

The gauge structure that emerges is manifest in the 
string correlator. The string correlator is precisely the 
gauge invariant correlator involving charges (e^ and Cj 
which are here portrayed by the spin components Sf and 
Sj) in the presence of a local photon gauge field A asso- 
ciated with the local gauge structure of the large set of 
states which become partially polarized under Us- The 
group element associated with this gauge structure in a 
path linking site i to site j is 



Qij = exp[i ^ Aki], 



(174) 



Enhanced, albeit decaying, string correlators with an un- 
derlying gauge-like structure also appear in doped spin 
chains [731. 1761. IttII . 

We may also generate maximal (Gij = 1) string corre- 
lators (case (iii) above); here, the number of states with 
maximal polarization is finite and no local gauge-like 
structure emerges. If we consider unitary transforma- 
tions {f/f ^'^l of the four GSs of the AKLT chain to the 
two orthogonal states given, in the basis by 



which form the endpoints of the chain. Here, the spin- 
spin correlation is maximal and equal to one (it is not 4/9 
as for the partially polarized case where only all non-zero 
entries attained a uniform sign). 

Having derived the string-type correlators of the AKLT 
5 = 1 chains, we briefly review and notice links with the 
FQHE problem. Just as in the spin chain problem, where 
a polarization operator led to long range order of string- 
like objects, there exists a similar transformation in the 
FQHE (see [H). 

Although conceptually clear, we wish to re-emphasize 
a simple yet important point. In the presence of entan- 
gled GSs, many transformations may map the original 
GS basis to new states in which the correlations, albeit 
not being maximal nor tending to a smaller finite con- 
stant for arbitrarily large separations between their end- 
points, are larger than in the original basis. In the spin 
language employed in this Section, these transformations 
do not lead to uniformly polarized states. Such transfor- 
mations are afforded and may be studied by, for exam- 
ple, Jordan- Wigner and other non-local transformations 



2. Non-local correlators in general gapped systems 

We next construct non-maximal string and brane type 
correlators in general systems with a spectral gap (of size 
A£') . To this end, we rely on the work of Ref . [93| which 
showed that the GSs of a gapped system can be made 
arbitrarily close to a generalized matrix product state of 
finite-size representation. Such a general product state 
is similar to what we wish to have in order to allow for a 
transformation to a highly polarized state which we just 
discussed. Let us consider Hamiltonians of the type 



(178) 



\P- 



1 ) 

max/ 

' ) 

max/ 



,1), 

1-1,-1,-1,- 



-1) 



(175) 



and to two other independent states of high polarization, 
e.g. 



\ymaxi 

^ ) 

max/ 



|l,-l,l,•••,(-l)^•••), 
|-i,i,-i,.-- ,(-i)'-+\ 



(176) 



with an inverted staggered spin at the j'th site of the chain 
then, following the same sequence of steps as undertaken 
in Eqs. (fT62)) . and ^63]) . 

\{Pma.\S^S]\praa.)\ = | (.g„ | " ^ 5| [7,"° ISo) I 

(177) 



= 1. 



The transformation Uf-^iax involves all sites of the lattice. 
Unlike the transformation of (|16ip . Eq. ()177p does not 
only contain fields between (and including) sites i and j. 
Eq. p77p is a string operator for the two sites i and j 



For example, in the AKLT Hamiltonian of Eq. (|158p . 
Hi = Si ■ Si+i + ^{Si ■ Si+i)'^. To construct the approxi- 
mate local projective states, we set {q is a constant) 



H,it) 



AE 



dtH,{t). 



\r2TTq J- 

We next define the quasi-local operators 
AE 



(179) 



i7f""=(t) = e 

Hloc — 



E 

j:\l-j\<lpr 



H, 



(180) 



In Eq. (|180p , the quasi- local Hamiltonian iJ/oc is the sum 
of {Hj} with the distance between sites i and j satisfying 
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d{i,j) < [Iproj — R)- They key point is that the quasi- 
local Hamiltonians Mi approximate the original Hi and 
enable the construction of a product state. Here, if we 
write all expressions explicitly, we find that the projec- 
tion operators which take us to the low-energy sector of 
{Mi} and their sums and with Lg*-"-*) low-energy normal- 
ized projected states defined itcratively as in [ost lead 
precisely to size n string or brane operators. These pro- 
jection operators play the role of the selection rules em- 
bodied in Eqs. (|157p . and (|169p . These expectation val- 
ues are those of a string or brane like operator Q within 
the GS manifold, 

(^(") I^H) ^ (5|Q(») |5) > 1 - 6„, (181) 

with bounded 6„ as defined in For general gapped 
Hamiltonians for which we do not know the GSs and 
the unitary transformations which will maximally polar- 
ize them, the string/brane correlators that we find in 
Eq. (|181[) . are the best that one can do. 



3. Low-energy selection rules in general systems 

We now return to make connections between (i) the po- 
larization/'pToiect'rve algorithm which we outlined above 
for generating string correlators in systems where no 
GLSs are present to (ii) systems with exact GLSs for 
the special case of the AKLT chain in which all GSs sat- 
isfy the d = selection rules of Eqs. ^M)) . and p70)l . In 
systems with no GLSs, general two-point and other cor- 
relators may be non-zero. This forces us to think of only 
string-like correlators. Obviously, unlike the situation in 
gauge and other theories (e.g. Kitaev's Toric code model 
of Eqs. (fTO|) . and (fTT|) to which we will return to mo- 
mentarily), general spin-spin and other correlators may 
be non-zero in the AKLT and many other systems with 
low-range string order, there are no GLSs symmetries 
which force these correlators to vanish by an application 
of Elitzur's theorem and its generalization. Nevertheless, 
the form of the GSs as adhering to low-d selection rules 
(as in Eqs. and pTD)) for the AKLT chain) enables 

maximal correlators of a form similar to that encountered 
in systems with exact local symmetries. Such an exam- 
ple is furnished by Kitaev's model of Eqs. (fTU]) . and ITTI) 
for which the GSs attain the form of Eq. ((95|) . Here, the 
general d — projection operator 



P. 



71-. 771-. 



B„) (182) 



is not only a symmetry but also acts as the identity oper- 
ator within the GS basis (similar to Eqs. (|169p . and (|170[) 
for the AKLT chain). The expectation value 



J5) = l, 



(183) 



for any GS \g). Here, the set {si,--- ,s„;pi, ,^ 
contains an arbitrary number of sites (n) and plaque- 
ttes (m) out of those in the entire lattice. In general 



the expectation value of Eq. (|183p leads to a product of 
string-like objects. At finite temperatures, this expecta- 
tion value will be degraded and we have 



.1 + tanh/3. 



7l + 7n 



(184) 



These brane-type correlators are closely related to (and 
are sums of) Wilson loop-type correlators [3 which are 
formed by the boundaries of connected clusters of pla- 
qucttc and vertex terms on the direct {Bp) and dual {Ag) 
lattices. 

Similar to the discussion prior to Eq. (|18ip . for more 
general gapped systems, we can construct and evaluate 
string correlators at finite temperatures. A more non- 
trivial example is furnished by &rane-type correlators re- 
cently found [oil for the Kitaev model on the hexagonal 
lattice. In that system, the spin basis GS [sl] is not of 
the form of Eq. (|76p with a uniform / = 1. 

Obviously, Kitaev's Toric model and related systems 
arc very special: most systems cannot have GSs which 
can be written in the simple form of Eq. (j95|) . Never- 
theless, in most entangled systems, there is a non-local 
unitary transformation which maps the GSs to states of 
maximal correlation. When written out in terms of local 
fields, this unitary transformation, which will generally 
involve all sites in the system, may lead to string- (or 
brane-) type correlators. In practical terms, for gapped 
systems, there is generally an effective matrix product 
form which allows to us construct string correlators as 
shown in this Section. 



B. The non-equivalence of string/brane orders and 

TQO 

Our general polarization (or selection rule) algorithm 
relies on the existence of a general non-local transfor- 
mations which maps the GSs onto uniformly factoriz- 
able states. The very generality of this algorithm sug- 
gests that the definition of TQO as given by Eq. ([1]) and 
the presence of string (or brane) orders are two indepen- 
dent issues. This distinction of the two orders is made 
clear by a simple example. The local expectation value 
{9a\Si^i\ga) depends on \ga), violating Eq. ([T]) and sug- 
gesting that the AKLT chain is not topologically ordered 
(see Appendix [G]) . 



XVII. ON GROUND STATE DEGENERACY, 
TOPOLOGY, AND ROBUSTNESS 

One of the earliest definitions of TQO was that as clas- 
sifying systems in which the GS degeneracy depends on 
the top olog y of the manifold on which the system is em- 
bedded [20| . In that way, the degeneracy depends upon 
a global property. In what follows, we will note that this 
earliest of definitions also encompasses systems which do 
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diniHy = dimT^. 



FIG. 14: A cartoon of the manifold of states {{\ya}}) where 
a local gauge structure exists. The double arrow above de- 
scribes the action of the non-local unitary operator Us (and 
U^^). Upon applying Ub, the states \ya) become partially po- 
larized, Ua\ya} = \Pa} = ||/q), aud have enhanced correlations 
relative to a chosen local observable (see text). By its very 
definition, the space Tiy — {\ya)} includes the GS manifold. 
Generally, it is far larger. In the AKLT problem, TCy is of 
size 2^^^+^ — 1, while the net number of GSs is no larger than 
four. The unitary transformation Ua preserves the dimension 
^2^^3+1 _ "i^-j The total number of possible linearly inde- 

pendent states in the 5 = 1 system is 3^° . The subspace Ti.y 
is an infinitesimally small fraction of all possible states (7i). 
States in TLy are linked to each other by local gauge-like oper- 
ators. [See Eq. (|169|l and the description of these gauge-like 
operators in the text.] 

not satisfy the robustness condition of Eq. ([T]) . Wc will 
illustrate this point by three examples. 

(i) Let us first consider the Klein model of Ref. [9^. 
By the theorem proved in [sl] all GSs of this model on 
the square lattice are superpositions of singlet dimer cov- 
erings. In these dimer coverings, at least one dimer must 
appear in each plaquette. An explicit calculation reveals 
that all of these special dimer coverings on the square 
lattice are orthogonal to one another. Thus, the degen- 
eracy is exactly equal to the number of dimer coverings in 
which each plaquette houses at least one dimer. Armed 
with this result, we may next compute the GS degen- 
eracy of this square lattice system when it is embedded 
on manifolds of different topology. These are easily tab- 
ulated: a system with open boundary conditions has a 
degeneracy equal to 

degeneracYKiein^open-bc = 2^"^+^ + (185) 

Here, Ns denotes the number of sites while Nd denotes 
the number of diagonals. 

A system with periodic boundary conditions (a system 
on a torus) has fewer states in the GS. The 3Ns defect 
states appearing in Eq. (jl85p are incompatible with peri- 
odic boundary conditions. For different sorts of periodic 
boundary conditions, the degeneracy can vary. Systems 
with periodic boundary conditions along the diagonals 
support more GSs than systems with periodic boundary 
conditions along the horizontal/vertical axes. 



Thus. Klein models of this type have a GS degeneracy 
which depends on the topology of the embedding man- 
ifold. Nevertheless, as we discussed in earlier Sections, 
the bulk of such Klein model GSs do not obey the ro- 
bustness condition of Eq. ([T]): In any given dimer state, 
we may compute the quasi-local expectation value {Sj-Sj) 
to see whether or not the spins at nearest (or, at worst, 
next-nearest) sites form a singlet dimer or not. 

Although in the above system, the GS degeneracy de- 
pends on the topology, it also depends on other local 
transformations. We next construct a system whose de- 
generacy depends only on the topology of the manifold 
on which it is embedded. 

(ii) Let us consider a classical D = 2 Ising model on a 
closed oriented manifold of genus g. Let us mark all of 
the spins by vertices of a polyhedra on the manifold. 
Let us next consider the Hamiltonian 

Here, V is the number of vertices - the total number 
of Si ~ +1 spins, a is an arbitrary positive constant, E 
is the number of edges (of "H — h" pairs on neighboring 
sites on the polyhedra) on the polyhedra formed by the 
spins, and F is the total number of faces. As by the 
Euler-Lhuillier formula for closed oriented surfaces with 
genus or number of handles g [Eq. (|120p ] , we will have 
that the GSs of H will correspond to all polyhedra (or 
Ising states) with a value of (X](ij) ^i^jY which is given 
in terms of the genus number g. For a given size (number 
of spins) of the Ising system, the degeneracy will depend 
only on g and nothing else. 

(iii) Let us finally consider the classical D ^ 2 Ising 
gauge theory of Eq. ([M]) . Similar to Eq. (|122p . and 
the discussion of Section IXIII A 3[ the degeneracy on a 
manifold of genus 5 > 1 is [86i] 2^-2'^s-i^ The large pref- 
actor of 2^= can be avoided in a legitimate local gauge 
fix which reduces the number of degrees freedom (cfj) 
yet does not introduce additional constraints on the pla- 
quette variables {Bp} beyond those seen in Eq. ((99|) . To 
this end, let us draw a closed non-intersecting path on 
the lattice which passes through the centers of all of the 
lattice plaquettes. This path is of length F - the num- 
ber of plaquettes. Let us set afj = 1 on all of the bonds 
which this path does not intersect [13 ■ By Eq. (|120p 
with V = Ng the number of vertices, our path leads 
to = {Ns + 2g — 2) uncrossed bonds. There remain 
two unfixed bonds (the bonds which are crossed by the 
path) in each plaquette - each bond being shared by two 
plaquettes. Following this procedure, the classical Ising 
gauge Hamiltonian of Eq. ((84|) reduces to that of a cir- 
cular Ising chain of length F which has two GSs. Such 
a procedure sets [Ns + 2{g — 1)] bonds to fixed values 
whereas there are only [Ng — 1] gauge fix degrees of free- 
dom: {2g— 1) additional Ising degrees of freedom appear 
in the Ising gauge system by comparison to that of the 
Ising chain. A legitimate gauge fix leads to a few plaque- 
ttes which host more than two unfixed bonds - (2g — 1) 
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formerly pinned bonds by our procedure need to be set 
active. This leads to a degeneracy 

degeneracy|3i,g ^3,^, = (187) 

in any legitimate gauge fix. This result is, as it must, 
identical to that of Kitaev's model [Eq. (|122p ] which, 
in its GS sector amounts to a particular (non-classical) 
gauge fix: Ag = 1. Local measurements can differentiate 
between the various GSs (thereby violating Eq. ([1])) of 
this system whose GS degeneracy depends on topology. 
Examples such as this can be avoided if we only consider 
rank-n TQO [see the discussion of Section IIII A| : That 
is, if we consider n GSs which satisfy Eq. ([T]). However, 
considering GSs which satisfy Eq. ([1]) defeats the purpose 
of an independent conjectured check on TQO - that in 
which the GS degeneracy depends only on the topology. 

Putting all of the pieces together, we note that the 
link between (1) the dependence of the GS degeneracy 
on topology and (2) the robustness condition of Eq. ^ 
is tenuous. 



XVIII. CONCLUSIONS 

To conclude, we analyzed Topological Quantum Order 
(TQO) from the perspective of generahzed Gauge-Like 
Symmetries (GLSs) and gauge structures. This approach 
allows us to identify and construct new physical systems 
with TQO. This work details, in great length, the abbre- 
viated results summarized in an earlier Letter [J . Below 
we list the new results which are provided in our work: 

(1) We provide an extension of the zero-temperature 
definition of TQO (Kitaev's original definition which re- 
lies on robustness to local perturbations) to finite tem- 
peratures. 

(2) We relate GLSs to general (non-Abelian) 
Aharonov-Bohm type transformations. We discuss how 
low-dimensional GLSs lead (i) to conserved non-local 
quantities, (ii) the absence of quasi-particle excitations in 
these systems, (iii) how these symmetries enable a natu- 
ral identification of charge fractionalization. and to (iv) 
the appearance, in high dimensions, of topological terms 
which are usually only associated with low-dimensional 
theories (dimensional reduction). 

(3) We show how finite-T TQO is mandated by low- 
dimensional GLSs. Our results are particularly strong 
for gapped systems with d <2 continuous GLSs. 

(4) We discuss the selection rules associated with these 
symmetries and show how these enable the construction 
of states with TQO. In one form or another, these se- 
lection rules account for all of the known examples of 
TQO to date. We stress that GLS selection rules alone 
can enforce a vanishing effect of local perturbations - no 
spectral gaps are invoked. 

(5) We further show how GLSs can enforce the ap- 
pearance of degenerate states in the large system-size 
(thermodynamic) limit. This degeneracy is borne by van- 
ishing of off-diagonal matrix elements (in the thermody- 



namic limit) of quasi-local operators in the basis which 
is spanned by eigenstates of the GLSs. 

(6) We identify the GLSs associated with all of the 
prominent examples of TQO and identify several new 
systems with TQO. It is important to stress that it is 
not necessary to have Hamiltonians with exact GLSs. 
The same occurs also in systems with standard (local) 
order parameters: although no real system displays per- 
fect symmetries, the fixed points associated with certain 
symmetry may have a finite regime of validity even when 
perturbations are added. What is important for TQO, 
much as in standard phase transitions, is that theories 
may be adiabatically connected with systems in which 
the GLSs are exact. The adiabatic link between Frac- 
tional Quantum Hall systems on a two-torus and the (m- 
rized)Peierls chain (in which d = 1 Im symmetries are 
exact) is an example of this principle. 

(7) We show that in several well-known examples of 
TQO (Kitaev's model and Wen's model) an exact di- 
mensional reduction occurs. As the partition functions 
of these systems arc equivalent to those of Ising chains, 
no finitc-T phase transition occurs. No less important, 
in spite of the existence of a spectral gap, some zero- 
temperature topological quantities are eradicated at any 
finite temperature if all sites in the system may freely 
change. A finite spectral gap might not protect the ro- 
bustness of topological quantities even at infinitesimal 
temperatures. 

(8) We show that in certain systems, the entanglement 
entropy might not always reflect TQO. To this end, states 
are constructed in which local perturbations lead to a null 
effect in the GS basis (and thus these states are topolog- 
ically ordered). It is shown that these states can have 
a vanishing entanglement entropy in spite of their TQO 
character. 

(9) We provide a general algorithm for the construc- 
tion of non-local string orders in systems where the form 
of the GSs (i) can be shown to adhere to certain selec- 
tion rules and (ii) specifically show how string or higher- 
dimensional brane orders appear in general gapped sys- 
tems in arbitrary dimensions. We identify an underlying 
gauge structure associated with some of these string or- 
ders. 

(10) We show that, on its own, the spectrum of a 
theory is not sufficient to determine whether or not TQO 
appears: The information about TQO is encoded in the 
state itself. In order to prove this, we show that the 
spectra of several theories with TQO (e.g. Kitaev's and 
Wen's models, gauge theories, etc.) arc equivalent to 
theories with standard (local) orders (Ising chains, Ising 
models in D~3 dimensions, etc.). 

(11) We demonstrate that the expectation value of the 
non-local string correlator of the AKLT problem (and 
other D = 1 models whose GSs belong to a particular 
class) is related to the expectation value of a local corre- 
lator of the nematic type. 

(12) We illustrate that systems can have a GS degen- 
eracy which depends on the topology of the surface on 
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which they are embedded yet not display TQO in the 
sense of robustness to local perturbations. 

(13) We remark on TQO on graphs and construct 
Gauge-Graph Wavefunctions with TQO. 

(14) We provide (see Appendix IC 2p an explicit Dirac- 
form expression for the symmetry operators of the (m- 
rized) Pcierls chain. 

(15) We introduce general SU{N) Klein models on 
small world networks and remark on their topological en- 
tanglement entropy content. 

Note added in proof. We were delighted to learn that 
some time after we reported the results concerning the 
spectra and thermal fragihty in [3l and C. Castel- 
novo and C. Chamon nicely confirmed our conclusions. 
Castelnovo and Chamon [90| found an entanglement en- 
tropy signature of the T = transition of the Kitaev 
model that we first found in Q and Q . 
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APPENDIX A: DEGENERATE PERTURBATION 
THEORY (DPT) IN A NUTSHELL 

We now present DPT in the Rayleigh-Schrodingcr ver- 
sion in a very simple and pedagogical fashion. This will 
allow us to establish the conditions (on the matrix ele- 
ments of the local perturbation V , between states in the 
degenerate eigenspace) for removal of such degeneracy in 
a certain order of DPT. We want to illustrate the differ- 
ence between these and the conditions of Eq. ([1]). 

The equation we want to solve perturbatively is 

(i/o + Af)|«') , (Al) 

where the two series 

n n 

I*) = ^^l*^''^) , and e = ^ X^e^P^ , (A2) 

p=0 p=0 

are assumed to be continuous and analytic functions of 
A, for A € [0,1]. 

Let us assume that we want to compute the correc- 
tions to a degenerate eigenspace Mo = {\ga)} {a € Sq = 
[1, Ng]) due to a given perturbation V. The unperturbed 
state ((^-Wl^-W) = 1) 

I*^°^)-E7"l5") (A3) 

a=l 



is a linear combination of the degenerate eigenstates \ga) 
satisfying (eq = 

Ho\ga) = eo\ga) , {9a\gp) = Sa(3 (A4) 

with coefficients that need to be determined using V 
for the perturbation expansion to be analytic. The solu- 
tion is subject to the normalization condition = 
1, which implies 

e = eo + («'^''^|AF|*) . (A5) 

The corresponding Rayleigh-Schrodinger 7T,-th order so- 
lution can be obtained recursively 

|<I»(»)) = |<i»{0)) + Coit - |«'("-i))(*W|)Ay|«'("-i)) 

(A6) 

with Go — (eo ^ ^o) ^ and II the identity operator. 

Let us establish the conditions on V to remove the 
degeneracy of Mq. Let's start assuming that the degen- 
eracy is removed to first order in V. That means that 
{ga\V\'Sf'-°^) = e(i)(.g„|«'(o)) = e(i)7„, or equivalently 

(^£^-e^'^)7o+ E 7/5K!^'=0, (A7) 

Va, with V^^J = {ga\V\gp). We need to determine e^^^ 
by solving det[W^^j3^]=0, with wi]} = V^i^ - e^^) and 
wj^p = V^^p , otherwise. It is clear that if 

Vl'J = Vpp, Va,/3 ,and 

V^}J = 0, V«,/3 («^/?) , (A8) 

then e*^^' = Va^ , and the degeneracy is not removed. If 
any of the conditions above is violated, one can compute 
e*^^) and the coefficients 7q,. At this point, at least part 
of, the degeneracy has been removed. 

Let's assume now that Eqs. (jASp are satis- 
fied and investigate whether degeneracies are removed 
to second order. That means that {ga\VGo{l — 
|5f{o))(^(o)|)y|^(o)^) ^ e(2)-y„, or (as can be shown as- 
suming Eqs. (K8^ ) 

{ga.\VGoV\¥'>^) = e'-'h^ , (A9) 

with Go = GqP±, and P± = X^r^es^ a projector 

onto the subspace orthogonal to Mo {Ho\v,^) = e^\v,^), 
{vri\v^) = Sri^). Therefore, the following conditions 

Vll^ = VfflK Va,/3 ,and 

V^'J = 0, Va,/3 {a^(3) , (AlO) 

with V^^p — {ga\VGoV\gp) , guarantee that the degen- 
eracy is not removed to second order. By induction, if 
the degeneracy is removed to order ^ -f 1 in perturbation 
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theory, then the following set of conditions need to be 
satisfied (Vn e [1,£] ) 



where 



V^;\ Va,/3 ,and 



with 



V^;> = 0, \fa,(3 (a^P) 

V^f = {ga\ yGoVyGoV \9p) , 
n factors V 



(All) 
(A12) 



and one of the following conditions 

V^'c^'^ ^ ypp^^\ for some a, /?, or 

^S^^' 7^ 0. for some a, /? (a 7^ /?) . (A13) 

For all these expressions to be analytic we have to assume 
that there is no level crossing, i.e. eq 7^ e,,, Vt; G 5(f = 
iV„ + !,•••. 



APPENDIX B: GENERALIZED ELITZUR'S 
THEOREM 

We review the results of Ref. |l7| needed for the pur- 
poses of this paper. Reference showed that "The ab- 
solute mean value of any local quantity (involving only 
a finite number of fields) which is not invariant under a 
d-dimensional symmetry group Qd of the i'-dimensional 
Hamiltonian H is bounded from above (as well as be- 
low for quantities of fixed sign) by the absolute mean 
value of the same quantity computed for a d-dimcnsional 
Hamiltonian H which is globally invariant under Qd and 
preserves the range of the interactions" . Non invariant 
means that the quantity under consideration, f{(t>i), has 
no invariant component: 



E/feik(0i)] = o. 



(Bl) 



For a continuous group, this condition is replaced by 
/ f[Si{'Pi)]dS = 0. To determine if SSB occurs, wc com- 
pute 



(/((/);)) = lim lim (/(0i)),,^, 



(B2) 



where {f{4'i))h.Ns is the mean value of f{(t>i) computed 
on a finite lattice of Ng sites, and in the presence of a 
symmetry breaking field h. Since A ~ Ui^i' the site i 
belongs at least to one set Cj . It is convenient to rename 
the fields in the following way: = i/^i if i ^ Cj and 
01 = 77i if i e Cj . The mean value {f{4'i))h,Ns is given by: 



E 



-/J(ff({V'i,rJi}) + hEieC- ')i> 



Z{,p} 



-/3(ff({Vi,r;i}) + 'iEiec. 'Ji) 



(B4) 



FromEq. (|B3|) : 



< 



Hi'} 



(B5) 

where {ipi} is the particular configuration of fields "01 that 
maximizes the expression between brackets in Eq. (|B3j. 
H({ipi,rii}) is a d-dimensional Hamiltonian for the field 
variables t/, which is invariant under the global symmetry 
group of transformations fjk over the field 771. We can 
define H{{r]i}) = H{{tpi,r]i}). The range of the interac- 
tions between the 77-fields in H{{r]i}) is clearly the same 
as the range of the interactions between the ^-fields in 
H{{(j>i}). This completes the demonstration of the main 
theorem in • Note that the frozen variables ^pi act like 
external fields in H{{(j)i}) which do not break the global 
symmetry group of transformations JTjk- 

Corollary I: Elitzur's theorem [4^. Any local quantity 
(i.e. involving only a finite number of fields) which is not 
invariant under a local (or c? = 0) symmetry group has 
a vanishing mean value at any finite temperature. This 
is a direct consequence of Eq. (|B3|1 and the fact that 
H{{i]{\) is a zero-dimensional Hamiltonian. 

Corollary II [13] ■ A local quantity which is not gauge 
invariant under a d = 1 GLS group has a vanishing mean 
value at any finite temperature for systems with finite 
range interactions. This is a consequence of Eq. (jBSP and 
the absence of SSB in onc-dimensional Hamiltonians such 
as H[{rii\) with interactions of finite range and strength. 
Here, /(jji) is a non-invariant quantity under the global 
symmetry group Qd [see Eq. (jBll) ]. 

Corollary III [l^. In finite-range systems, local quan- 
tities not invariant under continuous d = 2 symmetries 
have a vanishing mean value at any finite temperature. 
This results from [Eq. (|B5[) ] together with the applica- 
tion of the Mermin- Wagner theorem [3^ 



lim lim 



-/3(ff({v;i,T,i})-H/iEiec. Vi) 



= 0. 



Hi'} 



(B6) 

We invoked that Qd is a continuous symmetry group of 
is 3- non-invariant quantity for Qd [see Eq. 
(|B1[) ]. and il{{r]i}) is a d = 2 Hamiltonian that only 
contains finite range interactions. 

The generalization of this theorem to the quantum case 
is straightforward if we choose a basis of eigenvectors 
of the local operators linearly coupled to the symmetry 
breaking field h. Here, the states can be written as a 
direct product = lip) (g) \ri) . Eq. (jB5|) is re-obtained 
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with the sums replaced by traces over the states I77) 

-/3(ff({V;i,r,i})+/i2:ig(,. ,,0 

l(/(0i))/.,^J < ■'■ 



tr 



-0{H{{i:i^Vi}) + hT.;ec, ^0 



(B7) 

In this case, 1-0) corresponds to one particular state of 
the basis 1-0) that maximizes the right-hand side of Eq. 
(|B7p . Generalizing standard proofs, e.g. 9l|, one finds a 
T = quantum extension of Corollary III in the presence 
of a gap: 

Corollary IV . If a gap exists in a system possessing 
a d < 2-dimensional continuous symmetry in its low en- 
ergy sector then the expectation value of any local quan- 
tity not invariant under this symmetry, strictly vanishes 
at T = 0. As shown in the present paper, though lo- 
cal order cannot appear, multi-particle (including TQO) 
order can exist. The same holds for emergent discrete 
d < 1 symmetries for T = 0+. (Here, although T = 
SSB occurs wherein in some GSs local observables attain 
different values, SSB is prohibited for T = 0"*".) 

Corollary V |5l| . The absolute values of non-symmetry 
invariant correlators \G\ = KOiea- '^^^^ *^ Q ^^"^ 
bounded from above by absolute values of the same corre- 
lators \G\ in a d-dimensional system defined by Cj in the 
presence of transverse non-symmetry breaking fields. If 
no resonant terms appear in the lower-dimensional spec- 
tral functions (due to fractionalization), this allows for 
fractionalization of non-symmetry invariant quantities in 
the higher-dimensional system. 

Note. 

The existence of a spectral gap in a high-dimensional 
system may follow from an exponential bound in gapped 
low-dimensional systems in conjunction with Corollary 
V for all non symmetry invariant correlation functions, 
e.g. the connected G'^( I r|) < ^ exp(— |r|/^). The relation 
between spectral gaps and exponential decay of correla- 
tions was investigated in (9^ where it was proven that 
exponentially decaying correlations may imply a spectral 
gap. 



APPENDIX C: PEIERLS PROBLEM: 
POLYACETYLENE STORY 



THE 



In this Appendix, we first review the trimerized (m = 
3) polyacetylene chain (Appendix IC ip and then derive 
new universal symmetry operators for the general m- 
rized chain f Appendix IC 2[) . These symmetry operators, 
forming a group, endow the system with a primitive 
charge quantized in units of e/m. Thus, the Appendix il- 
lustrates the relation between symmetry, degeneracy, and 
fractionalization. 



1. The trimerized Peierls chain 

We start with the electron-phonon Hamiltonian for a 
chain of Ng sites introduced by Su, Schricffer, and Heeger 



(SSH) M 

H = -^tjj + l[4^C^_^.i^ + c]_^i„Cja] + E^iastic, (CI) 

where EeiasUc = Kf^ represents the (clas- 

sical) elastic energy with Uj denoting the displacement 
from equilibrium of the j-th atom. The operators Cj^ 

(Cj^) are the electron annihilation (creation) operators 
at site j (periodic boundary conditions are assumed, i.e. 

~ (^la)- By the SSH arguments, the kinetic hop- 
ping term ijj+i is modulated by the lattice displace- 
ments — for small separation between the atoms, the 
hopping element is enhanced whereas for large separa- 



tion it decreases. Thus, tj.j+i = to + a{uj — Uj+i) with 
a > 0. Next, we consider the trimerized case via the 
Born-Oppenheimer approximation and set 



Uj = ucos(— - 6) 

which leads to an elastic energy 
_ 3Ku^ 

^elastic — ^ 

and to hopping amplitudes 



(C2) 



(C3) 



ti 



to + a u sin( — 
o 

to + Vs a. u sin 6 



to + a/S a u sin( 



57r 



(C4) 



With this in hand, note that 9 I 
ti — * ^3; ^2 ti, tz 
Similarly, 6 6 — ^ leads to 

ti t2, t2 ^3, t3 



^ effects the shift 



t2. 



(C5) 



(C6) 



In what follows, a shorthand will be employed for the 
reciprocal lattice vector Q = ^ (the lattice constant a 
is set to one). Fourier transforming. 



1 



Cka 



(C7) 



the electronic portion of the Hamiltonian reads 

E 



keRBZ.cr 




,(C8) 



where RBZ is the reduced Brillouin zone with boundaries 
=F(5/2, and 




-ie 



Vi e 



ie 



V2 e 

£1 Vo e-^' 
e^' £2 



(C9) 
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with £n= — 2^0 cos(fc + nQ), and Vn = —i^/iaucos{k + 
nQ). The eigenenergies are solutions of the cubic equa- 
tion 

E^{k) - [itl + '^a^u'^]E{k) - 

cos3fc[i(%/3 aufsuiZe+^a^u^ to - 2t^^. (CIO) 

If the chain is 2/3 filled (i.e. only two thirds of the 
Ns sites are occupied by electrons) then only the lower 
most band is filled (out of the three bands) by spin up 
and down electrons. The net kinetic energy per electron 
is given by 

Tk^n = \{Ei{k = 0)+Ei{k = Q/2)) (Cll) 

with Eiik) the lowest of the three energy bands given by 
Eq. (|C10|1 . The net electronic energy is the number of 
particles (2A^s)/3 times the kinetic energy per electron, 
i.e. [Ekin = iTkinNs]- The total energy is 

Etot = Ekin + Eelastic- (C'12) 

Minimizing Etot with respect to 9 we find the three min- 
ima 9,nin = f I By Eq. (|C4|) . these correspond to 
the three possibilities 

I : ti = t2 = to H — ^ au, t^—to — Vs au 

W : ti = to — Vs au, t2=t3 = to H — ^ ctu 
\/3 

III : ti = ~ to ^ — — au, t2—to — au. (CIS) 

The electronic degrees of freedom ({Cj^, cj^}) are slaved 
to the distortions {uj} caused by the phonons. (This is 
indeed the essence of the Born-Oppenheimer approxima- 
tion.) 

2. New universal symmetry operators in m^rized 
Peierls problems 

To couch our results in a simple familiar setting, we 
will in the spirit of the earlier Appendix focus on the 
trimerized (m = 3) chain. Later on, we will illustrate how 
to extend these results for arbitrary m-rized chains and 
derive a universal m independent form of the symmetry 
operators. Before delving into the specifics, let us make a 
few general remarks. First, it is clear that the electronic 
only kinetic portion of each of the Hamiltonians 

Hk,n = j + l [c]aCj + la + '^j + la^ja] (CM) 

for the three distortion patters (a = 1,11,111) is non- 
degenerate. Henceforth, we will label Hkin for the hop- 
pings {tjj+i} dictated by \uj(ff)] by Ha. The two thirds 



full electronic state GS for each of the corresponding 
electron-only Hamiltonians is unique and given by 

I'Ac.)^ n /^UalO), (C15) 
fcGRBZ.cr 

where, in each case, we occupy all of the states of the low- 
est band (£'i(fc)). Note, however, that each single particle 
eigenstate /3|.j.g..Q|0) is different for each of the Hamiltoni- 
ans Ha. Notwithstanding, the single and multiple parti- 
cle spectra of Ha are all the same. This assertion follows 
from a relabeling of all indices (i.e. in going from II — > I 
we may define Cj+io- — *■ dj„ with a fermionic problem in 
dj„ whose Hamiltonian is given precisely by iJ] [d^^., rfj^] 

in lieu of 7J|[Cj^, cj^^]. As both sets of fermionic opera- 
tors obey the same algebra, the spectra of H\\ and H\ 
are identically the same.) The transformations between 
H\ — > H\\ — > H\\\ amount to discrete gauge transforma- 
tions. The basis vectors are permuted or relabeled (e.g. 
cjg. rf^g.) yet this change of variable is trivially innocu- 
ous in the spectra. A moment's reflection reveals that 
this gauge transformation is not a true symmetry of the 
Hamiltonian (e.g. \H\,H\^ ^ 0). 

Let us consider the full Hamiltonian (Eq. (jCip ) con- 
taining both the lattice degrees of freedom {%} and the 

electronic degrees of freedom {cj^ic]^}. It is clear that 
as the electronic states are slaved to the elastic deforma- 
tions, the existence of three GSs of the elastic deforma- 
tions leads to three GSs of the entire Hamiltonian. 

Next, we construct operators linking the different GSs. 
Consider the unitary operator 

Uei = n C^fetT , where 

fceRBZ,(T 
2 

Uka = exp[z y^(fc + nQ) nfc+„Q^], (C16) 

ri=0 

with Uka ~ c\a'^k<j- It readily verified that 

UeiHkuMuli ^ Hk^n{e + Q). (C17) 

In particular, it is clear that U^iH\U\i = H\\, 
U^iH\\Uli — H\\\, and U^iH\\\Ul, — H\ (whence we remind 
the reader that Omin = f 7 correspond to I, II, III). 

The full degeneracy between the different lattice dis- 
tortions and the electronic configurations slaved to them 
is captured by an operator relating 9 ^ 9 + Q for any 
of the pertinent values of 9 (and ensuing lattice distor- 
tions Uj given by Eq. (jC2[) ) and at the same time shift- 
ing the electronic operators accordingly. With the hop- 
ping amplitudes left explicitly as functions of 9 (i.e. 
^j-j+i ~ ^j(^) = ^j) the unitary operator 

Ug = exp[^QP0] = exp[g-^] (C18) 

shifts 9 by Q. [When the Baker-Campbell-Hausdorff 
identity is applied to Eq. (jC9[) . we attain UgHkinUe = 
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Q)]. With Eqs. (|C5p in mind, the operation 



UgHkin [tl, t2,tz\Ug — Hkin [^87 ^li ^2]- 



(C19) 



Thus, it is clear (by a similar relabeling as before with 
the fermionic c — > d) that, for any 0, the spectra of 
UgHkin[ti,t2,t3]Ul is identical to that of i?fci„ [^i , t2 , ^s] • 
By compounding both lattice distortions and opera- 
tions on the electronic operators we generate the full 
symmetry of the Hamiltonian 



U = UlU^i = exp[-lQPg] Y[ Uka 



(C20) 
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The operation of [/^ undoes the shift ^ + Q ef- 
fected by Uei. Trivially, UH{ti,t2,h)W = H{ti,t2,h) 
or, stated otherwise, [H, U] = 0. We immediately see 
that {1, ?7, U"^} form a Z3 group. As a byproduct at the 
special values of Omin = f 1 ^ the operator U links 
the degenerate Hamiltonians -fffei„ between states I II, 
II — > III, and III ^ I. [/ is a unitary operator living within 
each triplet of momentum values (fc. A: -I- Q, /c + 2Q) and 
corresponding 9 values {9,9 + Q,6 + 2Q). J7 is a direct 
product of 5[/(3) operators for each value of < fc < Q 
and < 6* < Q. Following the same prescription for other 
periods of the Peicrls distortion, we find a universal sym- 
metry which holds regardless of the value of m (m ~ 2 
in the dimerized case, m = 3 in the trimerized chain and 
so on). Very generally, the d = 1 connecting operators 



[/fco- , where Uka = ey^vVknu^]. (C21) 
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Here, the product spans the entire Brillouin zone (BZ) 
with k £ (— TT, tt] regardless of the size of the RBZ (deter- 
mined by m). In Eq. (jC21[) . the operator displacing the 
electrons by one lattice constant is Uei = exp[iPtot] with 
Ptot the total momentum. 

We just listed the symmetry operators for the elec- 
tronic degrees of freedom via the fc-space representation 
of Eq. (|C2ip . To make the string correlators slightly 
more transparent, we now transform Uei to real space. 
On a chain of length Ng (with Ns even). 



fceBZ 



(C22) 



with kernel {Ku = , y = Xj — Xi, and e = j^) 



R'li ^ — k 



^ Ns/2 
* n=-Nj2+l 



— —id,, — \ : 1 

e , TV^ey 2i e^'^ 

cos ■ sm 



which at eNg = 2tt reduces to 

e 



e'^y - 1 



("If, 



(C23) 
2 ^' 

(C24) 



and, in the thermodynamic limit (for y ^ Ns), simplifies 
to 



K{y) 

Putting all of the pieces together. 



with (for y < 



J|exp[i^/\, 



(C25) 



(C26) 



(C27) 



Taking into account the overall parity modulation 
(— 1)^*"^^ on the two sublattices, we have a universal 
(m-independent) form for the symmetry operators of the 
m-rized Peierls chain. 



1/^ 7!^ 



(C28) 



Here, Ura = Ci=2rcr,Vra = Ci=2r+la, "so- = Cj=2so-, and 

= Cj=2s+i(T are the electronic annihilation operators 
on the even/odd sites. The coordinates t/i = 2(r — s), 
2/2 = 2(7- — s) -I- 1, and j/3 = 2(r — s) — 1 which may 
further be recast in terms of the familiar Pauli matrices 
(or two-dimensional Dirac matrices encountered in field 
theoretic formulations of the SSH problem — similar to 
those appearing in massless chiral problems). 

The off-diagonal character of Eq. (jC28[) coupling even 
and odd sites is reminiscent to the coupling between such 
two chiral field in the field theoretical treatment of the 
SSH problem where the anomalous behavior and frac- 
tional Fermi number are partially brought about by off 
diagonal Dirac operators. 



APPENDIX D: A LIGHTNING REVIEW OF 
FRACTIONAL QUANTUM HALL SYSTEMS 

1. Symmetries 

Assume that the electronic system resides on a surface 
with periodic boundary conditions, i.e. a flat torus. The 
group symmetry operators that leave the Hamiltonian of 
Eq. (jl6p invariant are the magnetic translations which 
can be written as 



Tid) = Y[t,ia), 



(Dl) 



with a ~ a^ix + o-yBy a vector in the plane, ti{d) — 
cxp[jrki ■ a] , and ki = H^ + |i3 A the generators of mag- 
netic translations. With the condition that the magnetic 
flux is an integer number of flux quanta 00, i.e. N^(I)q, 



{T{Lyey/N^)T\ (D2) 
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with m, and n integers {L^Ly/N^ = 27r^^, where i = 
y^hc/{\e\B) is the magnetic length), are translations 
which respect the boundary conditions and satisfy 



T1T2 = exp[-i2nh']T2Ti, 



(D3) 



where the filling fraction v = N^lNip = p/q with p and q 
mutually prime integers. Most notably, for < 1, one set 
of these operators {Tfln^o,!,- or {T^}m=o^,■■■ ,q-i 
suffices to link all GSs to each other, i.e. the GS subspace 
is at least (7-fold degenerate. Since continuous d = D = 2 
symmetries of magnetic translations along a given space 
direction (which are a subset) suffice to link all GSs, TQO 
follows from our theorem. Therefore, there is no (local) 
order parameter that characterizes a FQH liquid state. 



To understand the topological degeneracy of real FQH 
liquids (i.e, including, for example, random potentials) 
Wen and Niu [23| assumed well-defined quasiparticle- 
quasihole elementary excitations with fractional statis- 
tics, and finite energy gaps. The symmetry opera- 
tions they considered correspond to d = 1 quasiparticle- 
quasiholc propagation T^y an integer (n) number of cy- 
cles around the toric axis (which are further related to 
integer flux insertions along each of the toric directions). 
It is assumed that such processes map GSs to GSs and, 
by definition, satisfy an equation similar to Eq. (|D3[) . 
The operators T linking the various states lead, for finite 
Lx.y, to an exponential splitting between the various GSs 
of the infinite L^.y system. The considerations discussed 
in [20j mirror those which we elaborated on in the Ap- 
pendix on the Quantum Dimer Model, although in the 
latter case we could write down the linking operators in 
terms of the original degrees of freedom. 



Clearly, real FQH systems are not exactly represented 
by the H of Eq. (fTB)) . The presence of impurities, for 
example, spoils the invariance of the Hamiltonian under 
magnetic translations. But this is not a problem as long 
as the perturbations added to H do not close the excita- 
tion gap. Of course, those perturbations (which are not 
invariant under magnetic translations) will remove the 
GS degeneracy; however, as we have shown, the split- 
ting will be exponentially small in the system size and 
will eventually vanish in the thermodynamic limit. This 
is indeed the essence of TQO, and the H of Eq. pB]) 
is a minimal fixed-point Hamiltonian capturing the TQO 
of FQH systems. We emphasize that the key point to 
prove existence of TQO in FQH liquids is the realization 
of d = 2 continuous symmetries connecting all possible 
GSs which, by our theorem, cannot be broken. . 



2. Reduction to an effective one-particle problem 

In terms of the effective field theory of [11], the La- 
grangian density is 



(p*{ido - ao - eAo)4i - -^(l)*{idi - ai ~ eAi)'^(j) 
2m 



inq 



(D4) 



Under the decomposition into global and local excita- 
tions, ai + eAi = + Sai{x, t), this leads to 

Z = /i..exp[./dt^(.4..-.4M 



This is invariant under the combined transformations |2C 







exp 
f 27rpi 



2TTi 



PlXl 

Li 
2np2), 



P2X2 \ 
L2 ) 



(D5) 



with integer pi and p2 such that 0' is single valued. The 
first of Eqs. (|D5[) is, in the thermodynamic limit, a con- 
tinuous d = 2 U{1) symmetry while the second is a dis- 
crete d = transformation. The two global coordinates 
(6*1,2) account for non-energetic excitations. Effectively, 
all short-range fluctuations 5ai can be thrown out. The 
resulting problem has only two fields 61 and 02- In the 
low-energy sector the Hamiltonian = H^ -I- H^ , where 
= de^ — iAe^ and the gauge field corresponds to m 
flux quanta threading the torus parameterized by (^i ,^2)7 
B ~ m/(27r). The Hamiltonian is now a single particle 
operator which is invariant under the magnetic trans- 
lation group symmetry. The magnetic translation group 
cannot be broken for this d — Q (single particle) problem. 
The degeneracy is clearly invariant under the magnetic 
translations linking the different GSs. This is a central 
result found by Wen and Niu [2^. Local deformations 
of the a field cost finite energy and can therefore be re- 
moved in the low energy sector. We stress that if we do 
not restrict ourselves to the lowest-energy sector then the 
d = 1 (or, more precisely, the d = 1 -I- 1) quasi-particle 
evolution around a toric cycle is the symmetry operator 
of the original many-body problem. 



Reduction to a one-dimensional problem with 
discrete d = 1 symmetries 



When placed on a thin torus, the D = 2 Quantum Hall 
system reduces, in the limit in which the thickness of the 
torus tends to zero, to a D = 1 problem [2l|, [l^. This in- 
finitely thin limit of the torus is adiabatically connected 
to that of the original D = 2 systems where both of its 
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sides {Lj. and Ly) are comparable. As recently noted 
in [U, [22] , anyons in Quantum Hall systems may be re- 
garded as domain walls in Z? = 1 systems for which we 
derived for the d = 1 symmetry operators of the Peierls 
chain in the previous Appendix. 



APPENDIX E: THE QUANTUM DIMER MODEL 
AND ITS RELATION TO KITAEV'S TORIC 
CODE MODEL 

Other systems have GSs identical in form to that of 
the the Kitaev model (Eq. albeit only at iso- 

lated points. These systems do not have Hamiltonians 
which conform to those of the generalized Kitaev type of 
Eq. ([86|l . The consequence of TQO generally follows for 
all states of the form of Eq. (17^ . One of these systems is 
the Quantum dimcr model [56| which recently has been 
looked at anew in the context of quantum computing in 
a sp irit very similar to Kitaev's Toric code model, e.g. 

M- _ 

On a triangular lattice this model takes the form [95| . 



+ «E(l ")( + • (El) 



Here, the sum on i runs over all of the Np plaquettes, 
and the sum on a over the three different orientations of 
the dimer plaquettes, namely V\ rotated by and ±60°. 
We refer to the plaquettes with a parallel pair of dimers 
as flippablc plaquettes. As a complete orthonormal ba- 
sis set we use {|c) | c = 1, • • • , iVc}, where |c) stands for 
one of the Nc possible hard-core dimcr coverings of the 
triangular lattice. V is diagonal in this basis, with 
V\c) = nfi{c)\c) measuring the number, nfi{c), of flip- 
pable plaquettes in configuration c. In this system, the 
d = I invariants correspond to parity conservation in the 
following way (sgI . [95l | : draw two non-contractible loops 
Ci^2 passing through the bonds and wind the system in 
the horizontal and vertical directions respectively. The 
parities of the number of dimer crossings (—1) '''•'''='-.^ 
reflect two d ~ I symmetries. On general lattices, when- 
ever t = V (the so-called "Rokhsar-Kivelson point" ) [5^ , 
the GS in each sector of {d = 1) topological numbers is a 
superposition of all dimer configurations- each appearing 
with an equal weight- this GS is precisely of the form 
of Eq. (j92p with the sum now over all dimcr coverings 
of the lattice belonging to a specific sector q defined by 
these d = 1 invariants q. The GSs of the Kitaev model 
are formally the same as those of the Quantum Dimer 
Model at the Rokhsar-Kivelson point of Eq. ([5^ with 
the states |c) one coding for all dimer coverings which 
belonging to a given topological sector- the GSs are all 
eigenstates of the topological d = I operators. Unlike 
the square lattice sibling of Eq. (jEip . on the triangular 



lattice there is a finite range of parameters for which no 
local order appears. [9^ Formally, in all these cases, the 
GSs can be expressed (up to an innocuous scale factor 
M) as a projection of reference states {\<j))) onto a given 
topological sector q, \g) = M x P^l^). 

We should nevertheless emphasize the difference in the 
form of the Hamiltonian giving rise to these GSs here and 
in the Kitaev model. We cannot view the RK model as 
a realization of Eq. ([86]) . Formally, we may set vV to be 
H and Gi to be the local symmetries (in fact, additional 
four dimer moves are necessary to link all dimcr con- 
figurations). Here, however, the maximal eigenvalues of 
{Gi} occur when we have a maximal number of flippablc 
plaquettes. However, a maximal number of flippablc pla- 
quettes elevates the energy V to a maximum. Thus, we 
cannot follow the same procedure as before. 

Ordered valence bond solid crystal GSs cannot be 
linked to one another by the exclusive use of d = 1 sym- 
metries and as a consequence SSB can occur and TQO 
need not appear. The RK model maps onto a fcrmionic 
model [sB]. The different topological sectors amount 
to different boundary conditions (symmetric or antisym- 
metric along the two directions) on the free fcrmion prob- 
lem leading to the partition functions Z±^±. The spec- 
trum is gapped in the case of the triangular lattice and 
is gapless for the square lattice [IBl. The split in energy 
between the different topological sectors is bounded from 
above by the value of the split for a one-dimensional ray 
[see, e.g. the ray along R where R maximizes Eq. (18) in 
Ref. [9^ {R = 0{L))]. For a finite size system, tunnel- 
ing between the various GSs of the infinite size system 
lifts the GS degeneracy. The bound is that corresponding 
to tunneling events in a d = 1 system which are medi- 
ated by a soliton. The relations of (see Appendix 
IB|) bound this split from above by the energy split gen- 
erated by such events in a d-dimensional system. In our 
context, it is this small value of the split A which leads 
to an exponential decay (in the system size) of all local 
quantities. In the thermodynamic limit, these tunneling 
events lead to zero modes. This exponentially small split 
between the system with the two different boundary con- 
ditions is precisely the same as that for a d = 1 gapped 
system. This forms a. d = 1 analogue of the bounds of 
17[ employed throughout our work. The event which 
bounds the expectation value of local observables is the 
propagation of a soliton along a d = 1 trajectory. On a 
torus, such solitons (visons) must appear in pairs. Simi- 
lar to the FQHE problem discussed in Appendix ID 31 the 
trajectory of a vison pair around a d = 1 cycle of the 
torus is the operation which links different GSs in the 
thermodynamic limit. 



APPENDIX F: VON NEUMANN ENTROPY OF 
SU{N) SINGLETS 

In what follows, we provide all of the details under- 
lying the result of Eq. (|150p and related combinatorial 
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relations in Section IXIVI 

To build some of the intuition, let us first consider an 
SU{2) singlet, 



with 



101) - |10) 



(Fl) 



Here \ab) denotes a product states of two particles: \ab) - 
\a)i (g) 16)2 (a, b = {0, 1}). The reduced density matrix 



Pi = tr 2\ip){M 



E («i^)(V'i«> 

= {04}2 



(F2) 



Substituting (|FT1) into Eq. (|F2| and employing 2(O|-0) 
and 2(11^) = ^|0)i we have 



Here, 



Pi = ^(|0>(0| + |l)(l|)= (^^ ? 



^r* = -tri[pi lnpi]=ln2. 



(F3) 



(F4) 



This value is what was employed in the SU{2) analysis 
of Section EVZl 

Next, let us turn to SU{3) singlets (a, b,c = {0, 1, 2}). 
Here, 



|012) + |201) + |120) - |102) - |210) - |021) 
V6 



(F5) 



The basis generated by particles two and three spans 3^ = 
9 dimensions. The resulting single particle density matrix 



1 



Pi 



mm 



|2>(2|), 



(F6) 



leads to S"™* = — tr i[pi Inpi] = In 3. Similar analysis for 
general SU{N) singlets results in S'f"* = \nN. 

Let us next turn to the computation of the entangle- 
ment entropy between two particles and the remaining 
{N — 2) particles in a general SU (N) singlet. For N = 2, 
we have the obvious result that S'|"* ~ 0. This is so as 



P2 = \'4'){''P\ and consequently S"! 



-tr i,2[p2lnp2] = 



since p2 is a density matrix for a pure state. Alterna- 
tively, this can also be seen by the reciprocity relation 
— tr [pAlnp^i] = — tr [pslnps] if the sets A and B arc 
complementary to each other: An B = and AU B = 1 
- the complete space. The complementary set to the two 
particle set in an SU{2) singlet is empty and consequently 
SI"* = = 0. 

When applied to SU{3), the reciprocity relation shows 
that S"!"* = SI1\ = Sf"' = In 3. For purposes of a gen- 
eralization that will later follow, let us now forgo this 
shortcut in this particular case and compute S*™* long- 



(F7) 



hand. Here, we have from Eq. (|F5P that 
P2 = tr 3|'0)(V'I, 



3(0|V) 
3(1|V'> 



|12)-|21) 

V6 
|20) - |02) 

V6 
101) ^ |io) 

^/6 ■ 



(F8) 



This then leads to 
1. 



P2= g[|01)(01| 

+ |21>(21|- 

-^[|oi)(io| 

+ |21)(12| + 



H10)(10| 
-|02)(02|- 
H10)(01| 
|02)(20| 4 



H12)(12| 
-|20)(20|] 
H12)(21| 
|20)(02|]. 



(F9) 



In the basis spanned by the states 
(|01), |10), |02), |20), |12), |21)) (in that order) this 
leads to the following block diagonal reduced density 
matrix. 



P2 





/I 


-1 















-1 


1 














1 








1 


-1 








6 








-1 


1 






















1 


-1 




Vo 











-1 


1/ 



(FIO) 



The basic block diagonal structure is that of the matrix 

-1 "1') • (Fll) 

which has the two eigenvalues e_(- = and e_ = 2. The 
non-zero eigenvalues of p2 arc thus £i — £2 = (3 = 1/3. 
This leads to S^* = - Ci hie; = In 3. 

In general, an SU{N) singlet lip) will belong to the 
totally antisymmetric representation of the permutation 
group Sn- Thus, it will contain A^! terms in its expan- 
sion. Of these, A!/2 terms will appears with a sign 
and Nl/2 terms will appear with a (— ) sign. All of the co- 
efficients are equal to 1/ VAT . The matrix p2 has ^ 

eigenvalues {ei} different from zero. All of these eigen- 
values are equal to one another and 



E In = - 



lint 



(F12) 



Here, e = me^/N\ = (2m)/Af! where m = (A^ - 2)! 

A 
^ 2 

A^ ' 

2 I . i iie tactor ot m origi- 
nates from the following considerations. Suppose that we 



Putting all of the pieces together, we have e = 1/ 
which leads to S*™* = In ( ) . The factor of 
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isolate the first two particles (1 and 2) and trace over the 
rest (3, 4, • • • , N). For each pair of states, e.g. (0,1), we 
will have the two kets jOlafec • • • ) and |10a6c • • • ). Here, 



, c, • • • 



{2, 3, 4, • • • , (iV - 1)} with a ^ 6 ^ c • 



The number of distinct states jOlafcc • • • ) that we have is 
(7V-2)! = m. 

Let us next consider the evaluation of S*™* in an 
SU{N) singlet with general k. S'^"* is given by 

5r* = -tri,2,-,fe[pfclnpfc] (F13) 

where 

Pk=tT k+i.k+2.-,N\^)i^\- (F14) 
We partition each state in the N particle basis, e.g. 



■N-1) . 

{N-ky. 

This leads to the following structure for 



(F15) 



Pk 



kl 



■ [k\ terms] [k\ terms] H ^ (F16) 



terms 



Here, we employed the shorthand 

[k\ terms] = \0l2---k- 1) 
- \102---k- 1) 



all other even perm, 
all other odd perm. 



The submatrix [fc! terms] [fc! terms] is of the form 



A 



One —One 
—One One 



(F17) 



where the square matrix One is given by 



One 





1 


1 • 




(I 


1 


1 • 










■ 1 








• 1 


Vi 


1 


1 • 


• 1/ 



(F18) 



The spectrum of the matrix A is given by 
Spec{A} = {0,0,0,- • • ,0,A:!}. 



k\-l 



On the other hand, pk is block diagonal with 



(F19) 



blocks of the type A. Therefore, the number of non- 

' N ' 



zero eigenvalues of pfc is ( ^ ) . All of these eigenvalues 



are equal to one another: ei 



(2 



This leads to 



N 
k 



^ ^ elne = lnf ^ ) ,(F20) 



as stated in Eq. pSOp . 



APPENDIX G: NON-LOCAL STRING 
CORRELATORS: THE CASE OF 5 = 1 SPIN 
CHAINS 

1. Symmetries 

Given an S = 1 chain of length Ng, we define the 
following global su{2) spin operators 



N, 



S'^ = Sj , with p. ^ x,y, z. 



(Gl) 



The global symmetry operators, exp[i7rS'^] and 
exp[i7rS'^] can be related to global symmetries of non- 
local order parameters. The latter non-local order pa- 
rameters have longer range correlations (the string cor- 
relators of the Haldane chain). 

With Us ^ n,<feexp[z7r5|5f], {Us = Ul = U'^) 



^ UsS]U-^ = S-JexpiiTT SI], 



k=j+i 



k=l 
k=l 



(G2) 



with global transformed operators S*^ = 

Let us now examine the operators Pz = exp[i7rS'^] and 
Px = exp[z7rS'^]. We claim that these operators are none 
other than exp[z7rS'^] and exp[i7rS'^], respectively. The 
calculation is elementary. First consider 

exp[iaS^] = 1 4- ^K^sina) + {S^f{cosa- 1). (G3) 

For a — TT, 

exp[i7rS'|] = 1 -2(51)2. (G4) 
By rotational symmetry {Sj = {S^ , Sj, S^)) 

exp[i7rS^_,- • n] = 1 - 2{S^ ■ nf . (G5) 

Let us turn to the evaluation of the various terms in 
Px.z- From Eqs. (jG2p . the exponential transforms as 



eMiTTS'^]^UsCy.p[iTTSm-^^Us{l - 2{S^^f)U-^{G&) 
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The only non-trivial term is 

k=] + l k=3 + l 

= {S^reM^2nY,S^,] = iS^r, (G7) 
k=j+i 

where we invoked that for any integer spin (including the 
current 5 = 1 of interest, see also Eq. (|G3P ). 



Thus, 



exp[i27rS'^] = 1. 



Px = exp[i7rS"^] = exp[i7rS"^]. 



(G8) 



(G9) 



By the same token, = exp[i7rS'^]. It is obvious that 
{Hx, Px,^z, Pz} satisfy a Z2 x Z2 group algebra. First, 
note that P^ z = %,z- Next, we immediately see that 
[Px,Pz] = 0. More explicitly. 



[Po^,P^] = [[1(1 -2(5D'), 11(1 -2(5D')]-(G10) 



fe=i 



fc=i 



As 



[{S^)MS^,f)]=Q for any J, fc. 



(Gil) 



Px and Pz commute and the set {tx,Px, ^z,Pz} form a 
Z2 X Z2 group. The physical meaning of the operators 
Px.z is very simple — a rotation by tt about the S'^'^ 
axis. It is evident that this is a symmetry for any global 
rotation-invariant spin chain. In particular all Hamil- 
tonians discussed in Ref. [8^ . which includes i^AKLT, 
display a global Z2 x Z2 symmetry. A rotation by tt 
about the internal axis may be viewed as the prod- 
uct of a rotation by tt about (the operator P^) and 
a rotation by tt about (Pz)- Thus, the rotation by tt 
about S"^ does not constitute a new group element and 
we merely have a global Z2 x Z2 group. Furthermore, 
PxPz ~ PzPx = Py = cxp[i7rS'^], with commutators 



[Us,Px] = [Us,Pz]^0. 



(G12) 



Moreover, each element in the product forming the non- 
local transformation Us commutes with the GS connect- 
ing operators Px.z, i-e. 

[eMi^S]S^],Px] = [exp[iTTS]S^],Pz] = 0. (G13) 

For the AKLT Hamiltonian (Eq. p55)) ) reveals that 

[i?AKLT,f/s] ^0. 



2. Local operators distinguishing the AKLT 
ground states and connecting operators 

In general, the degeneracy of the GS subspace depends 
on boundary conditions. For Haklt, the GS is four- fold 



degenerate in the case of open, and non-degenerate in the 
case of periodic, boundary conditions. We now construct 
local operators which lift the four-fold degeneracy of the 
AKLT GSs. We further comment on the relation between 
the role of the non-trivial string operator Us and the 
operators which links pairs of GSs. 

In the string transformed basis, the Hamiltonian 



is given by 



Haklt = UsHakltUs 



2l 



(G14) 



(G15) 



where 

h, = -SJSJ+^ + ^Je^^^^+^^'+i'S-J+i - 5|5|+i.(G16) 

The two-sites Hamiltonian [hj + ^h^] is easily diagonal- 
ized nil . The GS subspace is 4-dimensional and spanned 

by the states \4>a)j ^ \(t>a)j+i, 01 = 1,2,3,4, where 

I'/'i). = ^[|0>,+^/2|+),], 

\<j^2), = -i=[|0>,-V2|+),], 

|(/.3), = -L[|0),+V2h),], 

\<Pi), = -L[|0)^.-V2 h),], (G17) 



and 



1. 



1) 



[hj + M\(i^o.)j®\^o), 



+1 



H0a>j® |0a)j + l.(G18) 



In this way the (non-orthogonal) AKLT GSs can be com- 
pactly written as 



(G19) 



We may ask the question whether a local (or quasi- 
local) operator distinguishes two orthonormal AKLT 
GSs. To this end we need to consider, for example, the 
two orthonormal GSs 



\9i) 

\92) 



Us\~9i), 

Us 



V9 



1 



(G20) 



and measure the z-component of the spin in the first lat- 
tice site with the result 

l = {9i\S'i\9i)^{92\St\92) = ~l, (G21) 

meaning that we can certainly distinguish these 2 states 
through local measurements. This suggests that the GSs 
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of this S* = 1 spin chain arc not topologically ordered 
according to the definition of Eq. (jT]) . 

To see how these GSs are connected, we apply the sym- 
metry generators to them. In the (tilde) basis, the global 
Z2 symmetry generator is exp[i7rS'^] = Yiji^ ^ 2(5^)^). 
An explicit evaluation reveals that 

(l-2(5j=f)|0i), = -|03>,, 
(l-2(5;f)|02>, = -104),, 

(l-2(^p2)|^3)^. ^ 

(1-2(^^)2)104),. = -102),. (G22) 

Similarly, 



(1 




= 102),, 




(1 




= I01>,, 




(1 




= 104),, 




(1 


-2(^1)^)104), 


= 103),. 


(G23) 



Thus, the global Z2 x Z2 symmetry operators {Px,Pz} 
link the GSs to each other and are thus named connecting 
operators. 



3. String correlators and string operators as partial 
polarizers 

The objective of this Appendix is to flesh out the de- 
tails of the reason why the string correlator is always 
larger or equal to the usual Necl correlator. Wc use gen- 
eral classes of states to prove it instead of Hamiltonians; 
we find this procedure more transparent and also more 
easily generalizablc to other problems. In this way, we 
extend and highlight physical aspects of a previous the- 
orem by Kennedy and Tasaki [8a |. 

Consider the parent states represented by a string of 
length (Ns — AI) of ± followed by a string of length M 
of zeros, 



ImL) 



Wli) = |--f--f-...+00---0). (G24) 



states |i/^) with a = 1,2,- •• ,Nm, with the property 
ivU^Ai) = 5c,f3Sv^, and M\7jfj) = M\r^fj) (77(0 = M, ^, 



and the operator M = J2j=ii^ ~ i^j)'^) counts the num- 
ber of zeros). Moreover, they satisfy 

PzKl) = (-1) 



(G26) 



with, much as in earlier Sections, Pz = exp[z7rS'^], and 
Px = exp[z7rS"^] {exp[iTTSJ]\mj) ~ — 1(— w,)) for z-axis 
polarization rrij = 0, ±1). Thus, the effect of Px is to 
conserve the number of zeros and to flip -f to — and vice 
versa. In this way, one can construct states which are 
also eigenstates of Px and Pz 

Ivtii) = ^(|MM>±kM», (G27) 

Any state with a well-defined number M of zeros and 
can, most generally, be written as ((0a/|0a/') = 

5 MM' ) 



Nm 

hi) = ^^aalrili). 

^Nm |„ 12 



(G28) 



This state satisfies {J2a=i = 1) 

A^I0A/), 



M|0m) 
Pz\M 
510m) 



(-1)^-^10m), 
Cm I 0m), with 



(G29) 



^ I + (-)! if Ns-M is odd and j] = 
'^^ y if TVs - M is even ^ ' 

Lemma: For a general state |0a/), 

|(0M|5fQ,,5'||0M)| > |(0M|5f5||0M)|, where 

Q,, = n exp[z^5|]. (G31) 

i<k<j 

The proof of this assertion is immediate. By explicit 
evaluation. 



Ns-M 



M 



(0M|5fO.,5'||0A/) 



Nm 



(G32) 



In what follows, the appearance of the zeros in these 
states allows for partial polarization (vis a vis full po- 
larization if the zeros were absent). Similarly, consider 
states with all possible rcshuflling of the zeros. For ex- 
ample, the states 

Imm) = I+0- + - + 0---0), 

= I-0 + - + - + --- + 0---0), (G25) 

and numerous others. Thus, for a fixed number M < Ng 
fN-\ 

of zeros there are ( ^ M = Nm states \fJ.%j) and Nm 



where Mfj = {VMl^iQijSjl'ilM) ■ Here, i and j are ar- 
bitrary site labels, i, j = 1, • • • , iVg. If either Sf\r]1j) or 
Sj\ri'^j) vanish the expectation value above vanishes (as 
docs ('7"/|<5'j^<S'||?7"/)). There are four remaining cases: 

St\rili) = 5|k?M) = ^ Q^J = -1 , 

Therefore, 



j^o.^\ if St\iit,) = 0, or ^IhJ,) = ^ ^^33^ 
1—1 otherwise 
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while (5« = ivUSfS^lvt,)) 

itSf\7jfj)^Q,orS^\rj%j)^0 
Sr,^{+1 if S^I^m) ^ S^lvt,) . (G34) 

As the string correlator has a uniform sign, 

\{MS!Q^jS]\M\ > \{Ms^s]\M\- (G35) 



Notice that this proof is easily extended to the case 
of states which are also eigenstates of {\4>m±) = 

J2a=i^a\'r]'M±): thcy are not eigenstates of unless 
TVs — M is even). 

Let us now generalize these statements and examine 
what happens when the number of zeros (M) is not a 
good quantum number but remains being a good one. 
We define a class of states via = 1) 



with the property {J2m l/^J^/P = 1) 



(G36) 



(G37) 



We have the simple 

Corollary: For the class of states defined in Eq. (|G36p 

mstQ.,s]m>ms^s]m (G38) 

The proof of this assertion proceeds as for the Lemma, 
{^\S^Q,,S]\^) = \PMf{HI\S:Q^JS;\c^M), (G39) 

M 

and we apply the Lemma, leading to the inequality 
(IG38|. 

If we seek states j^*) which are also eigenstates of Pz 
(see text just after Eq. (jG9|) ) 



P;,!^') = ±1|*), (G40) 
we then have the following constraints 

P,|*) = +11^-) ^ 1^-) = f^Mlhi), (G41) 

M,(Ns-M)& even 

P,|*> = -1|1'>^|«')= ^/3m|0a/>. (G42) 

AI,{N^~M)e odd 

In order to understand the reason why one gets larger 
string (non-local) correlations than local ones we will 
investigate the effect of Us = Y\j<k^^^\^''^'^j^k\ = 
n,<fe(l - 2(5|)2(5-)2) on e.g. l/xj^). In what follows, 
we will use a shorthand notation for the states |?7f/) and 



consider only those modes involved in the action of the 
operators. We find that {\mjmk)) 



1 _ of qz cx\2 

1 ^{JjOf^j 



+> 



1 - 2{S]Sr' 



f \--) 



1 - 2(5| 5,^)2 



1 — 01 qx\2 

1 ^\DjDi^j 



l^2{S'^Slf\\Qmk) = lOmfc), 
1 - 2(5|5^)2] |m,0) =-K0) (to, ^0).(G43) 

Let us next consider three sites (iTOim^mfc)) with 
Uijk ~ exp[i7r SfS'J] exp[i7r S'f 5)?] cxp[i7rS'|S'^] 

t/,jfc| + -+> = - 
U^,k\-+-) = - 
C/.jfe| + -0) = -I + +0), 

f/*,fe|-+o) = -I--0), 

C/»jfc|0 + -) = -|o + +), 

C/»,fc|+0-) = I+0+). (G44) 

For a generic chain of length Ns , 

UsWli) = (-i)'^M + + + + + --- + oo---o), 



Ns-M 



UsWl,) = (-i)^^i 



M 

^00_^),(G45) 

M 



where the phase is ipi = Ns{Ns - l)/2 - M{M - l)/2. 
Thus, the string operator Us acts as a polarizer. 
For generic states {|?7f/)}, 



(G46) 



where the state (Wm)) denotes the fully polarized 



state of (Ns — M) spins having m 
M intercalated zeros. For example, 



(-l)'^^l +0 + + + + + 
(-1)^^1-0 



-1 (to = —1) and 

0---0), 
- • •0),(G47) 



and, in general, the phase is given by 
NsiNs - 1) 



(G48) 



where {io} spans the position of the zeros. 

GSs of known Hamiltonians represent particular exam- 
ples of the classes |^m) or For example, the GS of 
the t-Jz model [tI, [tI] is a particular |0m) state, while 
the GSs of the Hamiltonians considered by Kennedy and 
Tasaki [85| belong to the class I^P). For example, con- 
sider the AKLT GSs of ^aklt, Eq. ((GT5|) . It is quite 
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clear that they represent linear combinations of polarized 
states with different number of zeros M 

Iffc) = E/^mI0m>, with, 

M 

Nm Nm 

l4f^> = EI'"^m), and|0i^/)) = ^|z>2,). (G49) 

7—1 7—1 



From what we have just shown 



= J2 (^mUSIi) = E /3m I0m) , with, 

M M 
Nm 

= E(-1)^1mX/>, and 

7 = 1 
Nm 



7 = 1 



which unambiguously shows that they belong to the class 
of states \^). 

Similar constructs for operators Us polarizing given 
GSs and bringing them to uniform reference states having 
high correlations of one sort or another appear in doped 
Hubbard chains [t^, III, IzS] • its t-J^ approximant 
the system may be mapped onto a 5' = 1 chain similar to 
the one above by a Jordan- Wigner transformations fl^ : 
there the local state S"? = corresponds to the presence 
of a hole. Similar correlations also appear in _D = 2 t-J^ 
models [t^ , and in Quantum Hall systems [t^ . 



APPENDIX H: GAUGE-GRAPH 
WAVEFUNCTIONS 

In the text we made some analogies between systems 
with TQO and graphs. To complete this circle of ideas, 



we now illustrate how simple wavefunctions on a graph 
may be concocted to have TQO. 

Let us consider a graph which on any link has a gauge 
field. In what follows, we will focus on a graph composed 
of {2Ns) points — we label Ns of these site as "+" sites 
and the other Ns sites as "— " sites. We now consider a 
state 10) which has M = NsZ dimers (crf^ j_ = —1 if 
and i- are linked by a dimer; af^ = 1 otherwise) con- 
necting sites to "— " sites. Here, z is the coordination 
number of the graph. In what follows, we examine the 
case of z = 0{Ns). Given any such state {cf)) we define 
the graph wavefunctions 



1^ 



1,2/ 



1 



E 



(-1) 



(HI) 



The permutations P+,- permute only the Ns or 
Ns " sites amongst themselves — they do not inter- 
permute "+" sites with "— " sites. The phase factor 
(—1)^+^" denotes the parity of the combined permuta- 
tions. Inspection reveals that (ipi\il.'2) = 0. We claim 
that |'0i,2) exhibit rank-n = 2 TQO as we now show. It 
is apparent that {'fpi\(7pg''\ipi) = {ip2\crpl,'' \i'2) where {p,q) 
are any two points in the graph {p £ +,q £ —). It is also 
clear that for any product of {(Jpq} we will have equiva- 
lent expectation values in and | ■;/'-)■ If z = 0{Ns) 
then no local operation links to \il'2)- Any oper- 
ation which links the two states involves 0{Ns) graph 
links. With <5+^_ being being any permutation of the 
or "— " sites respectively, we have Q±\ipi) = \ipi) 
and Q±\ip2) = (— l)*^**^" |'02)- Any such permutation 
involves, at least, 2z gauge links. As z = ©(TV^), the vi- 
able order parameters Q± are not quasi-local and involve 
0{Ns) fields. 
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non-local. In some cases, see (i) below, local fields alone 
may also chose a given GS or a sector of topologically 



ordered GSs. Case (ii) - that of non-local fields - is more 
common. The core analysis of our work centers precisely 
on the general non-local d > 1-dimensional GLS opera- 
tors which enforce such boundary conditions. 
We list examples of fields for two prominent models: 
(i) Local fields: 

(a) Let us consider the simple case of a lattice gauge the- 
ory such as a standard D = 2 Z2 gauge theory (Eq. (|13|) 
with h — 0)) on a, square lattice with closed boundary 
conditions. 
Its Hamiltonian is: 

□ 

In this case we may consider 
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with As the local star operator of Kitaev's Toric code 
model [see Eq. (|lip]. For e = 0+, He favors states which 
lie in the four-dimensional GS sector of the Kitaev model: 
{ga}t^i - To distinguish between the four different states 
a non-local (d = 1) field is required. 

Conversely, if the gauge theory is on a square lattice with 
open boundary conditions then for e = 0"'' clearly 
favors the equal amplitude completely symmetric GS - 
the equal amplitude sum over all spin configurations in 
which the product of a" around any plaquette is one. 
This gauge system given by // is a canonical example of 
rank-n — 4 finite-T TQO. For any non-gauge-invariant 
operator V, {V} = (because of Elitzur's theorem). 
However, what is coupled to the infinitesimal symmetry 
breaking field is a local operator. Elitzur's original proof 
invoked precisely a similar type of argument, 
(b) Next, let us consider a system for which Eq. ^ holds 
for aU GSs. The Kitaev model (Eqs. and ^) on 

a system with amended open boundary conditions has 
a unique GS - the equal superposition of all of the 4 
GSs of the Eq. (|92|l . Here, by amended open boundary 
conditions [as opposed to the open boundary condition 
Kitaev system analyzed in Section IXIII Aj we allude to 
the system of Eq. (|10p which further includes all star 
operators As of Eq. (Illf) which have only two or three 
operator products when s is near the boundary: two spin 
product when s is a corner site and three spin product on 
all other boundary sites. Written equivalently, the GS in 
this case is the equal amplitude superposition of all states 
in the basis which have Bp — 1 for all plaquettes p. 
[The system of (Eqs. (|lUp . and pip is defined on a plane 
with open geometry by the inclusion, for all boundary 
sites s, of all bonds (whose number for boundary sites is 
three or two (for corner sites)) which have one these sites 
as their endpoint.] 

According to the classification of Section [III a system 
with short-range interactions which is subject to closed 
boundary conditions can be made into one with open 
boundary conditions (and vice versa) by the addition of 
quasi-local operators; these operators live on the bound- 
ary of the system with open boundary conditions or along 
an arbitrary cut of the system with periodic boundary 
conditions. Here, we see how quasi-local terms which are 
added to Kitaev's model on a torus (that with closed 
boundary conditions) favor order in a particular GS - the 
unique GS of the system with open boundary conditions, 
(ii) Non-local fields: In the Kitaev model (Eqs. (|10p . 
and pip ) and many other systems, a non-local field fa- 
vors ordering in a given state. This non-local (d — 1- 
dimensional) field lives on the D = 2 boundary and 
directly emulates the effect on applying boundary con- 
ditions which favor order in one or another topological 
sector (characterized by d = 1 invariants). In this work, 
we precisely analyze the characteristics of such d = 1 
symmetry operators. 

As shown in the text, Kitaev's model has GSs given by 
Eqs. ((92]). Here, Ha = i/-(/iiZi -1-/12^2) with Zi given by 
Eq. p2p favors ordering in the GS of H whose eigenvalues 
of are given by sign(fti). Here, the operators Zi are 

non-local. On an L x L lattice, they span 0{L) fields. 
In general, the fields in Ha favoring the order in the state 
\ga) can be gleaned from the order parameter which for 
TQOs spans a volume in Dtopo < D dimensions. For the 
Kitaev model, this has Dtopo = 1. The procedure for de- 
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